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1 . Introduction 


Nearly  every  machining  operation  has  a potential  for  Improvement  in 
productivity.  This  potential  arises  from  the  fact  that  the  optimum  machining 
conditions  vary  with  the  Job,  the  machine  tool,  the  cutting  tool,  and  the 
operator. 

Machining  conditions  for  a particular  job  are  often  selected  by 
Taylor's  tool  life  equation,  through  the  use  of  handbooks  such  as  the 
Machining  Data  Handbook  [ 14]  , or  on  the  basis  of  the  engineer's  or 
operator's  experience.  These  methods  in  conjunction  with  trial  runs  allow 
the  determination  of  machining  conditions  which  apparently  produce  satis- 
factory results.  mTthese  machining  conditions,  iiowever,  are  often  "ball  park" 
estimates  of  the  optimal  conditions.  They  are  in  a sense  similar  to  the 
estimates  of  the  optimal  operating  conditions  for  a chemical  facility 
obtained  by  pilot  plant  operations.  A process  of  "tuning"  still  remains 
to  be  done.  However,  the  production  personnel  may  leave  the  process 
"untuned"  in  order  to  concentrate  their  efforts  on  more  pressing  problems. 

The  objective  of  this  study  is  to  develop  and  to  implement  a system 
for  optimizing  machining  conditions  for  single-operation  and  multiple- 
operation  (numerically  controlled)  machine  tools.  The  measures  of 
productivity  considered  are: 

1.  Production  cost  per  piece,  and 

2 . Production  rate . 

This  system  is  designed  to  use  data  collected  by  production  personnel 
and  to  provide  feedback  so  that  the  machining  parameters  can  be  adjusted 
accordingly.  This  system  is  based  on  the  concept  that  the  machining  process 
not  only  generates  the  configuration  of  a part,  but  also  generates  the  neces- 
sary metal-cutting  information  on  "how  to  machine  the  part.  " A primary 
consideration  in  the  development  of  this  system  is  that  the  process  under  study 
should  not  initially  be  drastically  changed  or  upset.  Thus  changes  in  operating 
conditions  should  be  minor  rather  than  major.  This  is  accomplished  through 
a set  of  rules  for  normal  operation,  so  that,  without  serious  danger  of  loss 
through  the  manufacture  of  unsatisfactory  parts,  an  evolutionary  influence 
is  at  work  which  steadily  and  automatically  moves  the  process  toward  its 
optimal  operating  conditions. 

This  planned  program  of  perturbations  of  the  process  variables  of  a 
machining  process  differs  from  programs  of  planned  statistical  experimenta- 
tion in  the  following  two  major  ways  [ 9 ] : 
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1.  The  program  is  conducted  on  the  machining  process  during  actual 
production  of  a product  which  is  expected  to  be  shipped  to  the 
customer.  In  a planned  experiment,  the  testing  is  usually 
conducted  in  a laboratory  or  pilot  plant  during  product  develop- 
ment. Often  a scaled-down  version  of  the  final  product  line  is 
involved.  Since  actual  production  is  being  carried  out,  the 
amplitude  of  perturbations  introduced  during  the  program  is 
generally  small,  and  sometimes  effects  may  only  be  determined 
statistically.  In  a standard  experimental  program,  the  amplitude 
of  perturbations  are  often  maximized  to  determine  the  effects 
resulting  more  expediently. 

2 . The  program  is  frequently  conducted  over  an  extended  time 
period  rather  than  on  a one-time  basis  as  in  the  case  of  most 
standard  experimental  programs. 

Certain  process  characteristics  are  favorable  to  this  evolutionary 
manner  of  operation,  while  others  are  unfavorable.  Table  1.1  gives  a 
list  of  some  of  these  characteristics.  It  is  particularly  important  to  note 
that  low  volume,  non-repeating  job  shop  orders  are  not  suitable  candidates 
for  this  type  of  optimization.  These  jobs  are  usually  completed  before 
sufficient  information  is  gathered  to  suggest  optimum  operating  conditions. 

The  optimization  system  includes  computer  programs  which  analyze  the 
data.  The  computer  output  is  then  technically  evaluated  by  a committee. 

The  major  task  of  this  committee  is  to  discuss  the  implications  of  current 
results  and  to  suggest  potential  changes.  This  committee  should  be  composed 
of  production  personnel,  such  as  the  foreman  and  an  operator,  as  well  as 
staff  personnel,  such  as  an  industrial  engineer  and  a tool  engineer.  This 
committee  should  meet  on  a regular  basis  to  review  the  current  operations 
and  to  suggest  future  studies. 

The  technical  exchange  which  occurs  in  these  meetings  can  be  more 
important  than  the  information  provided  by  the  computer  programs.  For 
example,  the  programs  can  suggest  new  speeds  and  feeds  at  which  to 
operate  and  can  be  used  to  compare  different  tool  materials  and/or  different 
types  of  tool  inserts.  However,  they  can  not  suggest  a new  tool  material. 
These  types  of  suggestions  must  come  from  personnel  familiar  with  the 
particular  problem.  The  programs,  however,  help  to  motivate  these 
discussions,  and  the  improvements  are  accomplished  through  the  informa- 
tion exchange  which  takes  place  at  these  meetings.  The  programs  can  also 
be  used  to  evaluate  and  to  document  the  performance  of  different  cutting 
tools.  This  documentation  could  serve  as  a justification  to  procure 
proven  tools  of  particular  brand s^at  possibly  higher  costs  for  similar  future 
applications. 
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Table  1 . 1 


Process  Characteristics  Which  are  Favorable 
or  Unfavorable  to  the  Use  of  EVOP[  9] 


• Characteristics  Which  Are  Favorable 

Characteristics  Which  Are  Unfavorable 

1.  The  process  involves  high 
volume  production  over  a 
t reasonably  extensive  time 

period . 

1 . The  process  is  a job  shop  with 
few  or  no  repeats  of  units  with 
identical  specifications. 

1 . The  potential  benefits  of 
process  improvements  are 
large  (the  process  is  an 
important  one  and  is  not 
already  operating  at  optimum 
conditions) . 

2.  The  cost  for  process  improvement 
exceeds  the  potential  benfefits. 

3.  The  process  variables  can 
be  perturbed  readily. 

3.  The  process  variables  cannot  be 
perturbed  readily. 

4.  The  process  stabilizes 
rapidly  after  a process  change. 

4.  The  process  requires  a long  time 
to  stabilize  after  a process  change. 

5.  The  process  response  can  be 
obtained  rapidly. 

5.  The  process  response  is  not 
obtained  rapidly  (for  example,  if  the 
response  variable  is  time  until  failure 
on  a life  test) . 
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Two  computer  programs  are  given  in  this  report  for  optimizing  machining 
parameters,  the  Performance  Index  Method  (PIM)  program  and  the  Machining 
Optimization  (MACHOP)  program. 

The  PIM  program,  referred  to  as  the  "on-hand"  program  in  the  contract, 
was  available  from  AWC^  as  a computer  listing.  It  was  designed  for  opti- 
mizing machining  conditions  of  single-operation  machine  tools.  During  the 
initial  phase  of  this  project,  the  PIM  program  was  modified,  debugged,  and 
tested  on  simulated  data.  During  this  same  period,  production  data  was 
collected  at  the  Rock  Island  Arsenal  (RIA)  Operations  Division  shops.  Section 
2 of  this  report  contains  a discussion  of  the  problems  encountered  with  the 
PIM  program  and  the  recommendations  concerning  its  use.^ 

Because  of  the  difficulty  in  collecting  data  for  the  PIM  program  and  in 
using  the  program  with  shop  data,  the  MACHOP  program  was  developed.  The 
latter  program  can  be  used  to  optimize  single-operation  machine  tools  as  well 
as  multiple  operation  (numerically  controlled)  machine  tools. 

The  MACHOP  approach  is  described  in  Section  3.  Section  4 describes  the 
use  of  MACHOP,  including  the  data  handling  system.  Section  5 contains  the 
MACHOP  analysis  of  the  data  collected  at  the  RIA  Operations  Division  shops. 

Appendix  A contains  a sample  output  of  a regression  program  which  was 
used  to  verify  the  regression  modules  of  the  PIM  and  MACHOP  optimization 
programs.  Appendix  B contains  a simulation  program  listing  and  sample  out- 
put. This  program  was  developed  for  preliminary  analysis  of  the  logic  in  the 
PIM  and  the  MACHOP  programs.  Appendix  C contains  the  data  collection  forms 
for  the  MACHOP  program.  An  analysis  of  the  PIM  design  module  is  given  in 
Appendix  D.  Appendix  E contains  sample  output  from  the  Carboloy  Systems 
Computerized  Machinability  Program.  A program  listing  for  the  MACHOP  program 
is  given  in  Appendix  F.  Program  documentation  necessary  for  the  implementation 
and  maintenance  of  the  MACHOP  program  is  included  as  Appendix  G. 


2.  The  Performance  Index  Method  (PIM)  Program 

The  first  phase  of  this  contract  concerned  the  adaptation  of  the  "on- 
hand"  program  (herein  called  the  PIM  program)  for  use  with  single-operation 
machine  tools.  This  phase  of  the  contract  required  approximately  3 months. 

During  this  period,  production  operations  in  the  Rock  Island  Arsenal 
(RIA)  Operations  Division  shops  were  surveyed  and  data  collection  schemes 
were  determined.  Simultaneously  the  PIM  program  was  converted  from  a 


*U.S.Army  Weapons  Command,  now  U.S.Army  Armament  Command. 

? 

Section  2 can  be  bypassed  without  loss  of  continuity  if  the  MACHOP 
program  is  of  primary  interest. 
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program  listing  to  punch  cards.  The  program  logic  was  checked  and 
corrected,  and  the  program  was  tested  on  simulated  data.  The  results 
of  this  testing  indicated  some  additional  programming  errors.  Following 
the  correction  of  these  errors,  the  program  performed  satisfactorily  on 
the  simulated  data.  Details  are  given  in  Sections  2.1  and  2.2. 

Major  problems  were  encountered  in  the  collection  and  analysis 
of  production  data  with  the  PIM  program.  A discussion  of  these  problems 
is  given  in  Section  2.4. 

2.1  Recommendations  on  the  Use  of  the  PIM  Program 


The  problems  which  may  be  encountered  in  using  the  PIM  program 
restrict  its  applicability.  On  the  basis  of  the  experience  obtained  in 
using  the  PIM  program  on  RIA  data,  the  following  recommendations  are 
made: 

1.  The  PIM  program  should  be  used  only  for  those  machining 
processes  or  operations  where  the  production  personnel  are 
readily  able  to  find  nine  feed-speed  combinations  (three  levels 
of  feed  for  each  of  three  levels  of  speed),  which  can  be  run 
without  disturbing  the  production  process  or  risking  the 
production  of  scrap  parts. 

2.  The  validity  of  the  performance  index  prediction  outside  of 
the  region  where  data  has  been  taken  is  highly  questionable. 
Hence  the  range  of  prediction,  which  is  specified  in  the 
program  by  the  usable  speeds  and  feeds,  should  be  restricted. 
In  particular  the  usable  feeds  (speeds)  should  be  limited  to 
one  feed  (speed)  level  above  and  below  the  feeds  (speeds)  at 
which  data  are  to  be  collected. 

3.  The  data  collection  forms  developed  for  the  MAC  HOP  program 
should  be  used  for  collection  of  data  for  the  PIM  program. 

4.  The  PIM  program  can  be  used  for  multiple-operation  tools 
(numerically  controlled),  replacing  the  speed  and  feed  variables 
by  increments  of  the  speed  and  feed  overrides,  respectively. 

2.2  Program  Logic  Check 

The  PIM  program  was  available  as  a program  listing.  This  was 
converted  to  punch  cards  and  verified. 

Although  the  PIM  program  had  been  previously  "debugged,  " a few 
errors  were  noted.  A list  of  the  changes  made  in  the  program  follows: 
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1 . Definition  of  Performance  Index  (PI) 

The  PI  was  defined  in  the  final  Report  DAAF01-70-C-106a  [10]  as: 

PI  = Q.Pr  + (1  - Q)/Cu. 

In  the  computer  program  (subroutine  PERIND)  the  following  defini- 
tion was  used: 


PI  = (1  - Q).Pr  + Q/Cu. 


This  was  corrected  to  agree  with  the  report. 


2 . Programming  Enrors  in  the  PICK  Subroutine 


Following  the  determination  of  the  optimum  feed  and  speed  based 
on  a given  set  of  experimental  observations , the  program  selects 
the  next  set  of  feeds  and  speeds  by  one  of  two  methods  depending  on 
whether  the  previous  optimal  point  was  on  a boundary.  If  the  previous 
optimal  point  was  not  on  a boundary.  Ham  [10]  indicates  that  nine 
points  are  to  be  picked  around  the  previous  optimal  point,  subject 
to  the  maximum  feed,  maximum  speed,  and  maximum  horsepower 
constraints  and  tests  made  for  uniqueness.  See  Appendix  D for  an 
analysis  of  the  design  strategy. 


The  program  (PICK  subroutine),  however,  incorrectly  selected 
at  most  eight  points,  one  of  which  was  the  first  point  in  the 
previously  analyzed  set.  The  other  seven  points  (2-8)  are  shown  in 
Figure  2.1. 

• • • 

7 3 5 


Feed 

(Ipr) 


2 (previous  optimum) 


• © • 

8 4 6 


Speed  (rpm) 


Figure  2.1  Points  Selected  by  Design  Module 
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To  be  consistent  with  the  report  and  to  achieve  a reasonable 
design,  line  #15  of  the  program  was  changed  to  NPT=0  and  line 
#25  was  changed  to  IF(J.EQ.l)  GO  TO  75. 

3 . Format  Statements 

Some  format  statements  in  the  program  were  changed  to  allow 
costs  of  over  $10  per  piece  to  be  printed  out. 

2.3  Program  Testing  with  Simulated  Data 

Prior  to  obtaining  plant  data,  a simulation  model,  which  is  discussed 
below,  was  constructed  for  testing  the  PIM  program.  Some  of  the  logic 
errors  reported  in  Section  2 . 1 were  first  noted  when  the  data  from  these 
simulations  were  analyzed  by  the  PIM  program. 

2.3.1  The  Simulation  Model 


The  simulation  model  was  based  on  the  turning  operation  on  the  Monarch 
lathe  (RIA,  ID  #30303)  for  a recoil  cylinder  (part  #10895646).  This  model 
was  constructed  for  the  preliminary  evaluation  of  the  PIM  program.  This 
simulation  is  different  from  that  discussed  by  Ham  [ 10]  . He  used  Taylor’s 
tool  life  equation  to  determine  the  tool  life  and  thus  the  cost.  Then  a 
random  error  term  was  introduced  into  the  cost  function.  In  the  simula- 
tion model  given  here,  individual  tool  life  data  are  generated. * On  the 
basis  of  these  tool  lives,  the  number  of  pieces  produced  and  the  number 
of  edges  used  during  one  shift  of  operation  are  determined.  Thus,  the 
data  are  in  the  same  format  as  the  production  data. 

The  simulation  model  was  developed  according  to  the  following 
assumptions: 

1 , The  equation  used  for  determining  the  tool  life  is 

VTaF0  = C, 

where  a. , fj  , and  C are  constants , 


D = Diameter  of  work  piece  in  inches , 
N = spindle  speed  in  rpm. 


*The  tool  life  model  does  not  include  the  possibility  of  tool  breakage 
and  the  accompanying  loss  of  unused  edges. 
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F = feed  in  ipr,  and 
T = tool  life  in  minutes. 

If  one  takes  logarithms,  the  above  tool  life  equation  can  be 
written  as 


In  T = — [In  C-8  In  F - In  V ] . 
a 

A random  error  term,  e,  is  added  to  the  above  equation  to 
introduce  variability  into  the  tool  life.  The  error  term,  e,  is 
assumed  to  have  a normal  distribution  with  mean  zero  and 
standard  deviation  of 

- [In  C - 8 In  F - In  V ] *ERR, 
a 

where  ERR  is  the  percentage  error  appropriate  for  the  process. 

2.  An  eight  hour  day,  allowing  for  the  operator's  personal  time, 
is  assumed  to  consist  of  420  minutes  of  production  time. 

3.  If  at  the  end  of  a shift,  a part  is  over  75%  finished,  it  is 
assumed  to  be  completed. 

A flow  chart  based  on  the  above  assumptions  is  given  in  Figure  2.2. 

The  computer  program,  written  in  FORTRAN  IV  and  a sample  output  are  given 
in  Appendix  B. 

2.3.2  Analysis  of  Simulation  Data 

The  center  of  the  initial  set  of  points  selected  was  V = 220  rpm  and 
F = 0.0187  ipr.  Although  the  PIM  program  will  operate  with  as  few  as 
six  design  points.  Ham's  [10]  recommendation  of  nine  points  in  a geometric 
pattern  of  a 3^  factorial  design  was  followed.  A large  range  of  feeds  and 
speeds  were  available  on  the  Monarch  lathe.  However,  the  usable  speed 
range  was  limited  to  95  rpm  - 330  rpm  and  the  usable  feed  range  was 
limited  to  0.007  ipr  - 0.0337  ipr. 

The  data  generated  by  the  simulation  program  for  the  first  experiment 
are  given  in  Table  2.1. 

Based  on  these  data,  the  on-hand  program  selected  the  optimal  speed 
and  feed  as  V = 220  rpm  and  F = 0.0337  ipr.  Note  that  this  is  on  the 
usable  feed  boundary.  Hence,  feed  is  no  longer  adjustable.  The  points 
suggested  for  the  next  experiment  and  the  simulated  production  data  for 
these  points  are  given  in  Table  2.2. 
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Figure  2 . 2 Simulation  Flow  Chart 
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Table  2 . 1 First  Data  Set 


Speed 

192 

Feed 

0.0168 

No.  of  parts 
14.0 

Time 

420.0 

No.  of  edges 
8.0 

192 

0.0187 

14.0 

420.0 

11.0 

192 

0.0210 

15.0 

420.0 

14.0 

220 

0.0168 

14.0 

420.0 

16.0 

220 

0.0187 

15.0 

420.0 

20.0 

220 

0.0210 

16.0 

420.0 

21.0 

255 

0.0168 

14.0 

420.0 

36.0 

255 

0.0187 

15.0 

420.0 

32.0 

255 

0.0210 

16.0 

420.0 

31.0 

Table  2.2  Second  Data  Set 

Speed 

166 

Feed 

0.0337 

No.  of  parts 
17.0 

Time 

420.0 

No.  of  edaes 
9.0 

220 

0.0337 

19.0 

420.0 

17.0 

290 

0.0337 

19.0 

420.0 

45.0 

330 

0.0337 

18.0 

420.0 

64.0 

Based  on  these  data,  the  PIM  program  selected  the  optimal  speed  and 
feed  as  220  rpm  and  0.0337  ipr.  Since  the  two  previous  optimal  points 
were  not  identical,  the  PIM  program  selected  a new  set  of  five  points  on 
the  feed  boundary.  The  feed  and  speed  settings  and  the  simulated  data 
for  these  5 points  are  given  in  Table  2.3. 

The  point  220  rpm,  0.0337  ipr  was  again  chosen  as  the  optimal  point. 
The  analysis  was  terminated  since  the  same  point  was  chosen  on  two 
consecutive  runs. 
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Table  2 . 3 Third  Data  Set 


Speed 

126 

Feed 

0.0337 

No.  of  parts 
16 

Time 

420.0 

No.  of edges 
6.0 

166 

0.0337 

18 

420.0 

2.0 

220 

0.0337 

18 

420.0 

22.0 

290 

0.0337 

19 

420.0 

46.0 

330 

0.0337 

18 

420.0 

63.0 

A graphical  presentation  of  these  results  is  given  in  Figures  2.3 
and  2.4. 

The  PIM  program  required  18  observations  to  obtain  the  optimal  feed 
and  speed.  Table  2.4  gives  the  cost  per  piece  in  dollars  for  some  of  the 
usable  feeds  and  speeds  for  the  case  where  ERR  is  0. 


Table  2.4  Cost/Piece  ($)  as  Given  by  Taylor's  Equation 


\ Speed 
\(rpm) 

145 

166 

192 

220 

255 

290 

330 

FeedN. 
(ipr)  \ 

0.0153 

11.76 

10.89 

10.24 

10.47 

10.05 

10.44 

10.95 

0.0168 

10.78 

10.05 

10.24 

9.72 

9.38 

9.74 

10.22 

0.0187 

10.78 

10.08 

9.51 

9.07 

9.38 

9.74 

10.22 

0.02 10 

9.95 

9.36 

8.88 

9.10 

8.79 

9.14 

9.58 

0.0240 

9.24 

8.74 

8.35 

8.53 

8.30 

8.60 

9.02 

0.0337 

8.11 

7.71 

7.42 

7.58 

7.41 

7.67 

8.47 

11 


Feed 

(ipr) 


Figure  2.3  Design  Points  for  the  PIM  Method 


*the  circled  numbers  indicate  the  experiment  number, 

i.e. , i indicates  that  this  point  is  in  the  i* **1  experiment. 
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Feed 

(lpr) 


Figure  2 . 4 Cost/Piece  Estimates  ($)  at 
Observed  Feeds  and  Speeds 


*The  superscripts  refer  to  the  experiment  number. 
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2.4  Program  Testing  with  Rock  Island  Arsenal  (RIA)  Data 


This  phase  of  the  study  required  a major  portion  of  time  and  dollar 
expenditures  of  the  research  project.  All  previous  experimentation  using 
the  PIM  program  was  based  on  simulated  data  and/or  a very  limited  amount 
of  laboratory  experimental  data. 

A turning  operation  on  a recoil  cylinder  using  the  Monarch  lathe  was 
selected  for  the  first  study.  Other  information  concerning  this  operation 
is  given  in  the  job  description  in  Table  2.5. 


Table  2 . 5 Job  Description  for  First  Study 
Operation:  - Turning 
Part:  - Variable  recoil  cylinder 
Material:  - Steel  tube  4140. 

Dimensions:  47.5"  long,  8.5"  dia. , depth  of  cut  1/8" 

Job  order  #:  0016011 


Part  #:  10895646 


Cutting  tool:  - Titanum  coated  carbide  insert,  multi-edged 


Tool  cost/edge:  - $.42 


Machine  tool  specifications:  Monarch  stepped  lathe. 

50  HP.  RIA  ID  #30303 

Available  speeds  (in  rpm):  84,  95,  110,  126,  145,  166,  192,  220,  255,  290, 

330.  380.  435 


Available  feeds  (in  ipr);  .0032,  .0035, 

.0060,  .0065, 
.0105,  .0120, 
.0187,  .0210, 
.0374,  .0421, 
.0673,  .0748, 
.1224,  .1346, 

Labor  and  overhead  rate:  - $ 18/hr. 


0037, 

.0038, 

.0042, 

.0047, 

.0013, 

0070, 

.0073, 

.0076, 

.0084, 

.0093, 

0129, 

.0140, 

.0146, 

.0153, 

.0168, 

0240, 

.0259, 

.0293, 

.0306, 

.0337, 

0451, 

.0518, 

.0561, 

.0585, 

.0612, 

0841, 

.0962, 

.1036, 

.1122, 

.1171, 

1496, 

.1683. 
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A 3^  factorial  experiment!  was  conducted  with  V = 220  rpm  and 
F = 0.0187  ipr  as  a center  point.  Each  of  the  nine  feed-speed  combina- 
tions was  used  for  one  shift  of  operation.  The  same  operator  was  used 
for  all  of  the  points.  The  data  in  Table  2.6  were  recorded  by  the  machine 
operator. 


Based  on  the  above  production  data,  the  cost  per  piece  in  dollars  is 
calculated  in  each  case  and  is  given  in  the  last  column  of  Table  2.6.  In 
this  case  the  performance  index  is  taken  as  1/Cu.^  The  PIM^ program  fits 
a regression  equation  of  the  form  PI  = bf 


+ bjV 


+ b2F 


fb3V' 


+ b4I* 2 


+ b5v. 


to  the  data.  (See  Ham  [10]  fora  more  complete  discussion.)  The’ performance 
index  is  evaluated  at  all  usable  feed-speed  combinations. 


Table  2.6  Data  Summary  for  Experiment  1 

Results  for  Recoil  Cylinder  - First  Study 


Speed  (rpm) 
V 

Feed  (ipr) 
F 

No.  of  parts 

Production 

time 

No.  of  tool 
edqes 

Cu 

$/piece 

192 

0.0168 

14 

434 

17 

9.51 

192 

0.0187 

18 

410 

18 

7.28 

192 

0.0210 

14 

377 

18 

8.63 

220 

0.0168 

13 

396 

23 

9.90 

220 

0.0187 

15 

393 

15 

8.28 

220 

0.0210 

20 

440 

44 

7.52 

255 

0.0168 

10 

217 

17 

7.22 

255 

0.0187 

11 

267 

29 

8.39 

255 

0.0210 

14 

357 

31 

8.58 

The  optimal  feed-speed  combination  is  selected  as  the  point  which 
maximizes  the  predicted  performance  index  PI  and  thus  minimizes  the 
predicted  cost,  Cu. 

*An  experiment  conducted  at  3 levels  on  each  of  2 factors,  resulting  in 
a total  of  9 experimental  points. 

2 

The  reciprocal  of  Cu,  the  cost  per  piece,  in  $/piece. 
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The  feed-speed  combination  selected  as  optimal  was  V = 95  rpm 
and  F = 0.0187  ipr.  At  this  combination  the  performance  index  predic- 
tion is  PI  = 0.212  and  hence  the  cost  estimate  is  Cu  = $4. 72/piece. 
However,  a few  calculations  illustrate  that  the  predicted  cost  cannot 
be  achieved  at  this  feed  and  speed. 

The  machining  time  per  piece  is  47.5/(95  x .0240)  = 20.8  min/piece. 
Allowing  nine  min/piece  handling  and  no  time  for  tool  change  we  arrive  at 
a total  time  of  29.8  min/piece.  Using  the  labor  plus  overhead  rate  of 
$ 18/hr  and  neglecting  the  tool  costs,  we  obtain  a conservative  estimate 
of  cost  per  piece  of  $9. 00/piece.  This  is  nearly  double  the  cost  predicted 
by  the  PIM  program.  It  is  also  higher  than  the  cost  per  piece  for  many  of 
the  original  nine  feed-speed  combinations. 

The  cost  estimates  Obtained  by  the  PIM  program  for  a number  of  feed- 
speed  combinations  are  given  in  Table  2.7.  (These  are  not  printed  by  the 
PIM  program  but  are  easily  obtained.) 

The  regression  calculations  in  the  PIM  program  were  also  verified 
by  analyzing  the  data  on  a separate  regression  program.  (See  Appendix  A.) 
No  substantial  differences  were  noted.  The  prediction  equation  explained 
37.3%  of  the  variability  in  the  data. 


Table  2.7  Cost/Piece  Predictions  ($)  Using  the  PIM  Program 


Speed 
N.  (rpm) 

95 

145 

166 

192 

220 

255 

290 

Feed  Ny 
(ipr)  \ 

0.0153 

14.59 

14.04 

12.75 

10.97 

8.75 

6.88 

0.0168 

8.67 

9.81 

9.90 

9.66 

9.01 

7.87 

6.61 

0.0187 

6.24 

7.46 

7.84 

8.10 

8.09 

7.65 

6.85 

0.0210 

5.10 

6.49 

7.08 

7.76 

8.32 

8.59 

8.19 

0.0240 

4.721 

6.69 

7.82 

9.53 

11.71 

14.98 

16.62 

I 


Indicates  the  optimal  point. 
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The  major  reason  for  this  failure  of  the  PIM  program  is  the  large 
amount  of  variability  present  in  the  data.  Due  to  this  high  variability 
in  the  shop  data,  the  prediction  equation  used  in  the  PIM  program  was 
valid  only  over  a small  region  of  feed-speed  combinations  where  the 
data  were  collected. 

One  method  of  reducing  the  variability  of  the  prediction  equation  is 
to  take  data  for  at  least  two  shifts  at  each  feed-speed  combination.!  How- 
ever, this  would  require  data  from  eighteen  shifts  of  operation.  This  does 
not  appear  to  be  a viable  alternative  because  of  the  excessive  length  of 
time  required  to  collect  the  data  prior  to  feedback. 

It  should  also  be  noted  that  the- specification  of  constraints  such  as 
the  horsepower  constraint  did  not  seem  to  be  of  any  value.  This  is 
probably  attributable  to  the  fact  that  the  mathematical  expressions  for 
the  constraints  are  only  approximations.  Determination  of  feasible 
operating  feed-speed  combinations  would  be  better  left  to  shop  personnel. 
These  points  can  be  determined  on  the  basis  of  experience  and  trial  runs. 

Another  problem  was  encountered  in  collecting  shop  data  for  the  PIM 
program.  The  program  requires  a minimum  of  six  feed-speed  combinations. 
(Actually  nine  combinations  are  recommended  in  the  users  instructions.) 

On  the  operations  observed  in  the  RIA  shop  it  was  difficult  to  select  three 
feeds  and  three  speeds  (yielding  nine  feed-speed  combinations)  without 
taking  a high  risk  of  producing  scrap  parts  at  one  of  these  combinations. 

Since  the  PIM  program  functioned  satisfactorily  with  the  simulated 
data,  some  commentary  seems  to  be  in  order  concerning  the  differences 
between  the  simulated  data  and  the  production  data. 

Whereas  the  simulation  model  assumed  that  the  handling  time  and  the 
total  production  time  were  constant,  namely  15  min/part  and  420  min/day, 
respectively,  the  quantities  were  highly  variable  in  the  actual  production 
situation.  The  handling  time  varied  from  9 min/part  to  15.8  min/part, 
while  the  production  time  varied  from  2 17 min/day  to  440  min/day. 

The  variation  in  the  handling  time  leads  to  much  higher  variability  in  the 
estimates  of  the  performance  index.  The  production  time  variability 
illustrates  the  need  for  recording  down  time,  rather  than  using  the  eight 
hour  work  shift  as  a time  base.  Another  major  difference  between  the 
simulation  model  and  the  real  production  situation  is  that  the  tool  lives 
for  the  simulation  model  were  determined  by  the  tool  life  equation  with  a 
20%  random  error  term  added.  In  the  shop  the  variability  was  even  higher 
due  to  tool  breakage.  Since  "triangular  throw  away"  inserts  were  being  used, 
breakage  of  the  tool  could  result  in  a loss  of  one  to  six  edges  at  one  time. 


The  complete  experiment  should  be  performed  using  the  same  operator, 
if  possible,  so  that  an  additional  source  of  variability  is  not  introduced. 
Data  from  long  runs  by  multiple  operators  may  be  averaged. 
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3,  The  Machine  Optimization  (MACHOP)  Program 


The  problems  encountered  in  the  collection  of  data  and  in  the  applica- 
tion of  the  PIM  program  indicate  the  importance  of  the  following  objectives 
for  development  of  a uaable  optimization  program: 

1.  A limited  number  of  feeds  and  speeds  should  be  attempted  and 
analyzed  initially  in  order  to  provide  early  feedback  and  to 
minimize  disruption  to  the  shop  operation. 

2.  Since  changes  in  feeds  and  speeds  must  be  gradual,  it  is 
initially  sufficient  for  the  optimization  program  to  determine 
the  direction  of  the  optimal  operating  conditions  rather  than 
selecting  a particular  feed-speed  combination  as  optimal. 

This  procedure  will  greatly  reduce  the  chances  of  incorrectly 
determining  the  optimum  feed-speed  combination.  This  is 
particularly  important  because  of  the  variability  in  shop  data. 

3.  At  most  two  variables  should  be  varied  simultaneously  to 
facilitate  the  usage  and  understanding  by  production  personnel. 

The  MACHOP  program  is  designed  to  accomplish  the  above  objectives 
in  addition  to  those  specified  in  the  contract.  No  constraints  other  than 
the  definitions  of  the  speed  and  feed  environment  are  considered  in  this 
process.  Although  certain  constraints  were  directly  incorporated  into  the 
PIM  program,  they  were  not  incorporated  into  the  MACHOP  package.  The 
data  collected  in  this  project  indicates  that  the  mathematical  expressions 
for  these  constraints  (i.e. , the  horsepower  or  the  surface  finish)  were  only 
crude  approximations.  Hence  it  seemed  better  to  rely  on  the  shop  personnel's 
experience  than  to  eliminate  a particular  operating  point  because  of  a 
constraint  equation.  In  addition,  these  constraints  are  not  as  critical  to 
the  MACHOP  program  as  they  were  to  the  PIM  program,  since  the  MACHOP 
program  moves  in  small  steps  toward  the  optimum  operating  conditions. 

To  simplify  the  collection  of  data  in  the  shop,  observations  are  always 
taken  in  sets  of  four  at  two  adjacent  feeds  and  two  adjacent  speeds.  For 
numerically-controlled  multiple-operation  machine  tools  the  feeds  and 
speeds  are  controlled  by  feed  and  speed  overrides.  This  procedure  avoids 
the  cost  of  preparing  new  NC  tapes  for  each  run. 

For  each  feed-speed  combination,  the  following  data  are  collected 
during  one  shift  of  operation: 

1.  the  number  of  parts  produced, 

2 . the  number  of  tool  edges  used  for  each  tool  type,  and 
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3.  the  production  time.^ 

These  results  are  submitted  to  the  MACHOP  program. 

On  the  basis  of  these  data  the  MACHOP  program  calculates  two 
performance  indices: 

1.  the  cost/piece  (Cu)  in  dollars/piece,  and 

2.  the  production  rate  (Pr)  in  pieces/minute. 

The  MACHOP  program  then  performs  two  types  of  analyses  for  each  of  the 
responses  or  performance  indices: 

1.  an  evolutionary  operation  analysis,  and 

2.  a response  surface,  regression  analysis. 

The  evolutionary  operation  analysis  is  based  on  the  work  of  Box  and 
Draper  [ 3]  . An  automatic  feedback  provision  for  systematic  optimization 
is  given  in  this  portion  of  the  program.  The  program  recommends  a set  of 
four  feed-speed  combinations  for  the  next  phase  or  cycle  of  operation. 

(This  set  of  feed-speed  combinations  may  be  the  same  as  the  set  just 
completed.)  The  evolutionary  operation  analysis  makes  no  assumptions 
concerning  the  form  or  shape  of  the  response  surfaces.  Its  purpose  is  to 
evaluate  the  differences  in  the  observed  responses  at  the  four  feed-speed 
combinations  as  compared  to  the  variability  of  the  process  and  to  suggest 
the  direction  of  movement  toward  the  optimal  machining  conditions . 

The  response  surface,  regression  analysis  uses  the  same  data  as  the 
evolutionary  operation  analysis.  The  parameters  of  two  regression  equations 
are  estimated,  and  these  equations  are  used  to  predict  the  responses  at  the 
feed-speed  combinations  in  the  region  in  which  data  have  been  collected. 
When  sufficient  data  are  available,  second  order  equations  are  fitted  to  the 
natural  logarithms  of  the  feeds  and  speeds.  Each  of  these  second  order 
equations  requires  the  estimation  of  six  parameters,  and  hence,  at  least 
six  different  feed- speed  combinations  are  necessary  to  estimate  all  the 
parameters. 

When  data  has  been  collected  at  fewer  than  six  feed- speed  combina- 
tions, lower  order  equations  are  used.  For  example,  for  the  initial  set  of 
four  feed-speed  combinations,  first  order  models  are  fitted  to  the  data,  and 
used  to  make  the  predictions.  In  addition  to  the  predictions,  analysis  of 
variance  tables  are  also  provided/  An  understanding  of  regression  analysis 
is  not  necessary  for  program  use.  Production  personnel  will  generally  be 
interested  only  in  the  predictions. 


^Actual  cutting,  work  piece  handling  and  tool  changing  time. 

2 

A complete  discussion  of  regression  and  analysis  of  variance  tables  is 
given  in  references  [ 4 ] and  [ 5 ] . 
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The  output  of  these  two  portions  of  the  program  serves  as  the  basis  for 
the  MACHOP  committee’s  discussion  concerning  the  next  phase  or  cycle  of 
experimentation.  In  general,  the  committee  should  follow  the  MACHOP 
suggestions,  unless  their  collective  experience  indicates  otherwise.  The 
prediction  output  gives  them  an  indication  of  both  the  unit  costs  and  the 
production  rates  at  adjacent  feeds  and  speeds. 

On  the  basis  of  this  information  the  committee  can  also  make  decisions 
such  as  the  following: 

1.  try  a new  cutting  tool, 

2.  try  a different  machining  process,  or 

3.  prepare  a new  tape  for  an  N/C  job. 

A new  cutting  tool  might  be  suggested  upon  collaboration  by  the  tool 
engineer  and  the  operator.  It  might  be  determined  that  a more  expensive 
cutting  tool  is  warranted  by  a corresponding  increase  in  production  rate 
and/or  reduction  in  total  cost.  A different  machining  process  might  be 
proposed  in  order  to  decrease  high  machining  costs  for  particular  parts. 

The  suggestion  to  prepare  a new  tape  might  result  from  the  fact  that  the 
MACHOP  program  suggests  that  an  override  be  increased  by  5%,  although 
the  upper  limit  has  already  been  attained.  Under  these  circumstances  it 
may  be  economical  to  prepare  a new  tape.  This  situation  could  also  occur 
if  the  program  suggested  increasing  an  override,  but'the  operator  felt  that 
this  change  would  have  a negative  effect  on  one  or  more  of  the  machining 
operations.  In  this  case  a new  tape  could  be  prepared  which  would  change 
the  feeds  and/or  speeds  of  a particular  set  of  operations  relative  to  the 
feeds  and/or  speeds  of  the  other  operations.  Following  any  of  these  changes, 
the  usual  data  should  be  collected  and  submitted  to  the  program  in  order  to 
investigate  and  to  document  the  effect  of  the  change. 

The  technical  aspects  of  the  MACHOP  program  are  discussed  in  the 
subsequent  sections.  Section  3.1  gives  information  on  the  performance 
index  calculations.  The  details  of  the  evolutionary  operation  portion  of 
the  program  are  given  in  Section  3.2,  and  the  response  surface-regression 
analysis  portion  of  the  program  is  discussed  in  Section  3.3. 

3.1  Computation  of  Performance  Indices 

The  two  performance  indices  used  in  the  MACHOP  program  are 
Cu  = Cost  per  piece  ($/piece) , and 
Pr  = Production  rate  (pieces/minute) . 


20 


These  quantities  are  calculated  for  each  of  the  test  points. 

For  a given  test  point  the  following  data  are  collected: 

Np  = Number  of  parts  produced, 

Nti  = Number  of  tool  edges  of  tool  i used,  and* 

T = Time  period  of  test  (in  minutes). 

The  test  will  usually  consist  of  running  at  a given  feed  and  speed  for 
one  shift.  Forms  for  collecting  these  data  are  provided  in  Appendix  C. 
Note  that  the  length  of  the  test  will  usually  be  less  than  the  shift  time 
because  of  interruptions  for  other  activities  such  as  safety  meetings  and 
personal  time. 

The  production  rate  computation  for  single-operation,  multiple- 
operation,  and/or  numerically  controlled  processes  is 

Pr  = Np/T. 

For  a single-operation  machining  process  the  cost  per  piece  is 
Cu  = [ (RLO)T  + (TLC)Nt  ]/Np 


where 


RLO  = Labor  and  overhead  rate  in  dollars  per  minute,  and 
TLC  = Tool  cost  in  dollars  per  edge. 

If  the  objective  is  to  minimize  tool  cost  per  piece  then  RLO  is  set  equal 
to  zero.  Alternatively,  the  overhead  may  be  removed  from  the  problem  by 
setting  RLO  equal  to  the  labor  rate. 


For  a multiple-operation  machining  process  the  cost  per  piece  is 


where 


n 

Cu  = [ (RLO)T  + 2 (TLC1)Nti]/N 
i=l 


n = Number  of  different  tools,  and 
TLCj  = Cost  of  tool  i in  dollars  per  edge. 

*Ntj  is  shortened  to  Nt. 
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3.2  Evolutionary  Operation  (EVOP) 


Evolutionary  operation  (EVOP)  is  a method  of  process  improvement 

that: 

1.  is  readily  conducted  under  actual  processing  conditions  by 
production  personnel, 

2.  operates  despite  the  presence  of  large  experimental  error, 

3.  provides  an  efficient  basis  for  scientific  or  technological 
feedback,  and 

4.  does  not  assume  knowledge  of  the  functional  form  of  the 
response  surface  or  any  explicit  knowledge  of  the  response 
function  except  that  it  is  smooth. 

EVOP  computations  and  procedures  for  analyzing  the  results  of  the  2^ 
factorial  experiment  are  based  on  the  work  of  Box  and  Draper  [ 3]  In 
this  analysis,  the  two  performance  indices,  Cu  and  Pr,  are  examined 
as  a function  of  the  feed  and  speed  of  the  machining  operation.  The  effect 
that  each  of  these  factors  exerts  on  the  performance  indices  is  evaluated 
and  two  (possibly  identical)  sets  of  operating  conditions  are  suggested  for 
future  operation. 

3.2.1  EVOP  Calculations 


A cycle  is  a set  of  four  observations  in  a rectangular  pattern  at  two 
adjacent  feeds  and  two  adjacent  speeds.  This  set  of  four  observations 
allows  the  determination  of  the  effect  of 

1.  feed, 

2.  speed,  and 

3.  feed-speed  interaction 

on  both  the  production  rate  and  the  cost  per  piece. 


*Box  and  Draper  also  consider  a 2^  factorial  design  with  an  additional 
reference  condition.  This  approach,  if  used  here,  would  require  taking  data 
at  three  feed  and  speed  levels.  Experience  with  the  machining  operations 
at  the  RIA  shops  indicated  that  this  was  too  wide  a range  and  could  lead  to 
production  of  scrap.  Hence,  the  22  factorial  design  without  reference 
condition  was  chosen. 
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A phase  consists  of  repetitions  of  the  same  cycle.  In  other  words, 
phase  indicates  which  set  of  four  points  is  being  observed  and  cycle 
indicates  the  number  of  observations  which  has  been  taken  at  each  feed- 
speed  combination  during  this  phase.  Cycles  are  repeated  within  a phase 
until  sufficient  information  is  gathered  on  the  effects  of  feed,  speed  and 
their  interaction  to  suggest  a new  set  of  feed-speed  combinations,  i.e. , 
a new  phase  of  operation.  The  new  phase  begins  when  observations  are 
taken  at  a new  set  of  feed-speed  combinations.  The  concepts  of  phase  and 
cycle  are  important  in  understanding  the  EVOP  analysis . A minimum  of  two 
cycles  of  observations  are  required  in  Phase  1 to  estimate  the  variability 
of  the  data.  Prior  estimates  of  the  variability  can  also  be  used  as  will 
be  discussed  later.  The  estimate  of  variability  is  necessary  to  determine 
the  significance  of  the  effects  of  feed  and  speed  on  the  responses.  Moving 
to  a new  phase  constitutes  the  movement  across  the  response  surface. 

Taking  another  cycle  means  additional  information  is  being  gathered  to 
determine  the  direction  of  movement. 

The  EVOP  computations  of  Box  and  Draper  [ 3 ] have  been  modified  to 
accomodate  machining  data  and  to  facilitate  computerization.  The  computa- 
tional procedure  is  summarized  in  Figure  3.1.  A complete  discussion  of 
these  calculations  is  given  for  phase  M and  cycle  n. 

The  first  step  in  the  EVOP  analysis  is  the  calculation  of  the  phase 
averages  and  phase  ranges.  These  computations  are  the  same  for  each  of 
the  performance  indices.  The  performance  index  will  be  denoted  by  y 
in  the  following  discussion.  Each  of  the  performance  indices  are  calculated 
for  each  feed-speed  combination,  for  which  data  are  available,  according 
to  the  equations  given  in  Section  3.1.  For  cycle  n the  new  observations 
at  each  of  the  four  feed-speed  combinations  are  denoted  by  y^n  (i=l  ,2,3,4). 
The  previous  cycle  sum  (PCSj)  at  each  of  these  four  points  is  given  by 

n-1 

PCSj  = S yjj . 


Thus,  the  previous  cycle  average  is 

yj  = PCS. /(n-1). 

The  four  differences  (d^)  are  calculated  by  subtracting  the  new  observa 
tions  from  the  previous  cycle  averages,  i.e. , 


d 


i 
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These  differences  are  used  to  estimate  the  variability  of  the  data.  The 
new  sums  (NSj)  are  computed  as  the  sum  of  the  previous  cycle  sums  and 
the  new  observations 

NS^  = PCSj  + yin. 

Finally  the  new  averages  are  calculated  as 

yA  = NSj/n 

for  each  of  the  four  feed-speed  combinations. 

The  new  averages  are  the  values  used  to  determine  the  effects  of  the 
factors  on  the  performance  index.  The  new  cycle  sums  and  the  new  cycle 
averages  become  the  previous  cycle  sums  and  previous  cycle  averages  for 
the  next  cycle  (i.e. , cycle  n + 1) , if  the  current  phase  is  continued. 

For  the  first  cycle  of  any  phase,  the  previous  cycle  sums  and  the 
previous  cycle  averages  are  zero  by  definition.  The  new  sums  and  new 
averages  (NS^  and  y^)  are  just  the  new  observations.  The  range  is  not 
meaningful  and  hence  is  not  calculated  for  this  cycle.  If  this  is  the  first 
phase,  there  is  no  measure  of  error  with  which  to  compare  the  effects  of 
the  factors.  This  is  the  reason  that  a minimum  of  two  cycles  in  the  first 
phase  is  required  to  produce  an  initial  estimate  of  the  standard  deviation 
of  the  error  of  the  response.  In  later  phases,  the  estimate  obtained  from 
previous  phases  is  used. 

The  effects  of  the  factors,  feed  and  speed,  are  calculated  using  the 
new  averages  y..  Consider  the  2^  factorial  design  configuration  as  in 
Figure  3.1.  Eacn  yj  is  associated  with  a point  i in  Figure  3.1.  Speed 
increases  from  left  to  right  as  depicted  by  the  arrow.  The  effect  of  speed 
is  given  by 

SPEED  EFFECT  = J ((  V2  + Yj  “ ( ^ + 74»  = j ( \ + Y 3 “ V 1 “ 74) . 

Feed  increases  as  indicated  by  the  labeled  arrow  and  the  effect  of  feed  is 
given  by 


FEED  EFFECT  = j ((  7g  + 74)  - ( 7j  + 73))  = | ( 72  + 74  - 7X  - 73> . 
The  effect  of  the  interaction  of  speed  and  feed  is  given  by 

SPEED-FEED  (INTERACTION)  EFFECT  = 

J «yi  +y2)  ■ (^3  +Va)]  =7{*1  +72  " ^3  ' ‘ 
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Figure  3 . 1 Evolutionary  Operation 
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The  variability  of  the  data  is  measured  by  the  standard  deviation  (S) 
of  the  observation  errors,  and  the  computational  procedure  to  obtain  it 
for  phase  M cycle  n is  described  subsequently.  The  new  S for  cycle 
n is  computed  as  the  range  multiplied  by  the  factor  f^  n (given  in 
Figure  3.1)  as  follows:* 

New  S = Range  * f^  n.  , 

where  Range  = [ MAX  (d)  - MIN(dJ  ] . The  new  sum  S is  then  computed  as 

i i 1 

New  sum  S = (S'  + S) , 

where  S'  is  the  previous  sum  S.  The  new  average  S (S)  is  equal  to  the 
new  sum  S divided  by  2: 

S = (New  Sum  S)/2 . 

As  indicated,  the  standard  deviation  is  a weighted  average  of  the  standard 
deviations  of  all  previous  cycles. 

When  a prior  estimate  of  the  standard  deviation  of  the  response  variable 
is  known,  it  can  be  used  as  a substitute  for  the  computed  standard  devia- 
tion during  the  initial  phase.  Having  reached  the  second  cycle  in  any 
phase,  however,  the  prior  estimate  of  the  standard  deviation  is  disregarded, 
and  the  computed  standard  deviation  is  used. 

For  any  cycle  (n  > 1),  the  standard  deviation  is  recomputed  as  indicated. 
For  the  first  cycle  of  a new  phase  (i.e. , n = 1),  the  standard  deviation  from 
the  last  phase  is  used.  For  n>  1 the  standard  deviation  is  updated  as  above. 

The  standard  error  (SE)  is  calculated  as  follows: 

SE  = S/Vn, 

where  n is  the  cycle  number,  and  S is  the  estimate  of  the  standard  deviation. 

3.2.2  Direction  of  Movement  on  the  Response  Surface 

Based  on  the  results  of  the  EVOP  calculations,  the  MACHOP  routine 
selects  a set  of  operating  conditions.  This  may  be  another  cycle  in  the 
current  phase  (i.e. , new  observations  at  the  same  four  feed-speed  combina- 
tions) or  a new  phase  may  be  initiated  (i.e. , observations  at  a different 
set  of  four  feed-speed  combinations) . In  either  case,  MACHOP  recommends 
the  four  points  (i.e. , two  adjacent  feeds  and  two  adjacent  speeds)  where 
new  observations  should  be  taken. 


The  factor  f^  is  used  to  convert  a range  to  an  estimate  of  the 
standard  deviation. 
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The  recommendations  are  derived  by  comparing  the  calculated  effects 
with  two  standard  errors  of  the  effect: 

1.  If  the  absolute  value  of  an  effect  is  greater  than  or  equal  to  the 
two  standard  errors,  t&e  effect  is  judged  significant,  i.e.. 

If  | EFFECT  | 5 2SE 

then  the  EFFECT  is  judged  significant. 

2.  If  the  absolute  value  of  an  effect  is  less  than  two  standard 
errors,  the  effect  is  judged  not  significant,  i.e. , 

If  | EFFECT  | < 2SE, 

then  the  EFFECT  is  judged  not  significant. 

When  one  or  more  effects  are  judged  significant,  a new  set  of  operating 
conditions  is  recommended  by  the  program  (i.e. , a new  phase)  in  accordance 
with  a decision  table.  The  decision  table  for  minimizing  Cu  is  given  in 
Table  3.1.  To  maximize  the  production  rate  (Pr)-,f  the  negative  of  Pr  is 
minimized.  Hence  the  same  table  is  used  by  the  program  with  the  production 
rates  replaced  by  their  negatives. 


3.3  Response  Surface  and  Regression  Analysis 


Regression  analysis  is  a technique  for  estimating  the  parameters  of  an 
equation  relating  a response  variable  to  a set  of  independent  variables.  The 
resulting  equation  is  called  a regression  equation.  In  the  MACHOP  program 
two  response  variables,  the  cost  per  piece  and  the  production  rate,  are 
considered  and  regression  equations  are  developed  for  each.  The  two 
independent  variables  are  the  feed  (F)  and  the  cutting  speed  (V). 


The  regression  equation  is  used  to  predict  the  response  (y)  for  the 
feeds  and  speeds  in  the  region  where  the  data  were  collected.  Five 
different  forms  of  the  prediction  equation  are  considered: 

(3.1)  y = bQ  + bj  In  V + b2  In  F 

(3.2)  y = bQ  + bj  In  V + b2  In  F + bg  In  V • In  F 

(3.3)  y = bQ  + bj  In  V + b2  In  F + b3  In  V • In  F + b4(ln  V)2 

(3.4)  y = bg  + b j In  V + b3  In  F + b3  In  V * In  F + b4 (In  F )2 

(3.5)  y = bQ  + bj  In  V + b2  In  F + b3  In  V • In  F + b4(ln  V)2  + bg (In  F) 2 , 
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Table  3.1  Decision  Table  for  Minimization  of  Cost* 


Significance  and 

Direction  of  Effects Action 

Reobserve 

same  Inc.  Deer.  Inc.  Deer. 

SpeedfV)  Feed(F)  Interaction  (VxF)  settings  Speed  Speed  Feed  Feed 
0 0 0 

P 0 0 

NO  0 

OP  0 

ON  0 

P P 0 

N P 0 

P N 0 

N N 0 


0 0 P 

0 0 N 


If  yj  is  min. , deer,  speed  and  feed; 
if  y2  is  min. , incr.  speed  and  feed; 
if  yg  is  min.,  incr.  speed  and  deer, 
feed;  if  is  min.,  deer,  speed 
and  incr.  feed. 


P 

P 

N 

N 

0 

0 

0 

0 

p 

p 

p 

p 

N 

N 

N 

N 


0 

0 

0 

0 

P 

P 

N 

N 

P 

P 

N 

N 

P 

P 

N 

N 


P 

N 

P 

N 

P 

N 

P 

N 

P 

N 

P 

N 

P 

N 

P 

N 


* 

* 

* 

* 

* 

* 

★ 

★ 

★ 

* 

★ 

★ 

* 

★ 

* 


★ 

★ 

★ 

★ 

* 

★ ★ 
★ 

★ 

★ 

★ 

* 

* 

★ 

★ 

* 

★ 


*P(N)  indicates  a significant  positive  (negative)  effect;  A 0 indicates 
the  effect  is  not  significant. 
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The  natural  logarithmic  transformation  of  feed  and  speed  was  deemed 
appropriate  on  the  basis  of  the  work  of  Ermer  and  Wu  [8]  and  the  experience 
gained  in  analyzing  RIA  shop  data.  The  program  can  be  easily  modified  to 
use  other  transformations.  One  of  the  equations  3.1  - 3.5  is  selected 
depending  on  the  number  of  data  points  available  and  the  number  of  levels 
of  the  feeds  and  speeds  through  the  logic  depicted  in  Figure  3.2. 

The  general  methodology  of  the  regression  program  is  briefly  explained 
in  Section  3 . 3 . 1 for  m independent  variables  and  n observations.  An 
understanding  of  this  section  is  not  necessary  for  use  of  the  program. 

Section  3.3.2  contains  a discussion  of  the  regression  calculations  for 
phase  1,  cycle  1. 

3,3.1  General  Methodology 

A multiple  regression  equation  may  be  expressed  in  the  following  form: 
y - y =0J  (x2  - Xj)  +j32(x2  - x2)  +...+  /3m(xm  “ *m)  + e' 
where  y and  x are  sample  averages. 

Given  the  set  of  observations: 


X = 
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xA  = E xij/n. 

j=l 
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Selection  of  the  Regression  Equation 
Figure  3.2 


The  b vector  contains  the  unknown  coefficients 


b = 


b 


2 


and  the  ^ vector  is  calculated  from  the  dependent  variable  observations  by 


where 


(yx  - y) 
<y2  - y) 

(yn  1 y) 


/ 


n 

y = E y/n. 
j=l 


The  least  squares  estimates  of  the  unknown  coefficients  are 

b = (X'XT^'y., 

and  the  resulting  prediction  equation  is 

y = y + bj(Xj  JTj)  +. . .+bm(xm  . 

The  analysis  of  variance  table  is  given  in  Table  3.2. 


Table  3.2  Analysis  of  Variance  (ANOVA) 


Degrees  of  Sum  of 
Source Freedom Squares 

Total  n - 1 y.'JL 

Regression  m b'(X'y) 

Residual  n - 1 - m y|y-b'  (X'y) 


Mean 

Square F 


b'(X'y)  b'(X'x)(n-l-m) 

01  mly'^-b' (X'y)  ] 

X'y-b'QC'y) 

n - 1 - m 
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The  significance  of  the  resultant  regression  may  be  tested  by  an  F test. 
Draper  and  Smith  [5]  state  that  the  equation  should  not  be  considered  a 
satisfactory  predictor  unless  the  F value  is  at  least  4 times  greater  than 
the  test  statistic  F(m  n_i_m  i_a\,  where  a is  the  significance  level  for 
the  test. 


The  ratio  (regression  sum  of  squares)/(total  sum  of  squares) , denoted 


, is  a measure  of  the  "goodness*1  of  fit  of  the  regression  equation.  An 
R2  = 1 indicates  a perfect  fit  of  the  data  to  the  function,  while  R2  = 0 
indicates  bj  = b£  =. . .=bm  = 0.  The  variance,  o2,  is  estimated  by  the 
mean  square  of  the  residuals. 

3.3.2  Illustration  of  Regression  Calculations  for  Phase  1,  Cycle  1 

To  illustrate  the  regression  methodology  we  consider  the  cost  per  piece 

prediction  for  phase  1,  cycle  1.  For  phase  1,  cycle  1,  the  following  informa- 

tion  is  available 

(m  = 

2 , n = 4) : 

V1 

F1 

NPTj  NTCj  NTj 

V2 

F2 

npt2  ntc2  nt2 

W 

F_ 

NPT  NTC  NT 

3 

3 

3 3 3 

V A 

F* 

NPT.  NTC  . NT. 

4 

4 

4 4 4 

where 

V 

= 

speed  (rpm) , 

F 

= 

feed  (ipr) , 

NPT 

= 

no.  of  parts  produced. 

NTC 

- 

no.  of  tool  changes,  and 

NT 

=: 

total  production  time  in  minutes. 

The  logarithm  of  the  feed  (In  F)  and  of  the  speed  (In  V)  are  the  independent 
variables.  The  X matrix  is 


X = 


(In  V,  - In  V)  (In  F,  - In  F) 

• • 

• • 

(In  V4-IhV)  (In  F4  - in  F) 


/ 


where  In  denotes  the  averages  of  the  logarithms. 
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The  b vector  is 


The  y.  vector  is  determined  from  the  equations  in  Section  3. 1 by 


The  b vector  is  given  by 


Yl  - y 


y2  -y 
y3  -y 
y4-y 


b = (X'X) 


The  prediction  equation  is  given  in  the  form 

y = y + bj  (In  V - In  V)  + b2  (In  F - In  F) . 

The  ANOVA  table  is  calculated  according  to  Table  3.2  using  the  expressions 
given  above. 

The  MACHOP  output  contains  an  analysis  of  variance  (ANOVA)  table, 
the  coefficient  of  multiple  regression  (R^) , the  regression  coefficients, 
their  t values  and  the  standard  deviations  for  both  the  cost  per  piece  and 
the  production  rate.  Predictions  of  the  cost  per  piece  and  the  production 
rate  are  made  for  all  feeds  and  speeds  within  one  level  of  the  extremes  of 
the  feeds  and  speeds  where  observations  have  been  taken. 
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4.  Use  of  MAC  HOP 


The  MACHOP  program  incorporates  evolutionary  operation  and  response 
surface-regression  analysis  to  determine  the  optimal  machining  conditions 
for  single-operation  and  multiple-operation  (numerically  controlled)  processes. 
It  is  intended  for  use  by  managers  or  foremen  directly  concerned  with  the 
operation  of  individual  machining  processes. 


4.1  Collecting  Data  for  MACHOP 

To  simplify  the  collection  of  data  in  the  shop,  observations  are  always 
taken  in  cycles — a set  of  four  observations  at  two  adjacent  feeds  and  two 
adjacent  speeds.  For  example,  if  the  feeds  and  speeds  selected  are 
(0.0210  ipr,  0.0240  ipr)  and  (192  rpm,  220  rpm) , respectively,  then  data 
is  collected  at  the  following  feed-speed  combinations: 


Feed 

(ipr) 


0.0240  4-- 

I 

0.0210 f - - 

I 

192 


I 


220 


Speed 

(rpm) 


The  MACHOP  committee,  composed  of  the  foreman,  the  appropriate 
technical  personnel  and  the  machine  operator,  selects  an  initial  set  of 
operating  conditions  (two  adjacent  speeds  and  two  adjacent  feeds)  for  the 
machining  operation  under  study.  Before  collecting  data  for  the  MACHOP 
program,  trial  runs  should  be  made  at  each  of  the  four  feed-speed  combina- 
tions to  insure  machinability  under  these  conditions.  If  problems  are 
encountered  during  trial  runs  at  any  of  these  four  feed-speed  combinations, 
a new  set  of  four  points  should  be  selected,  omitting  the  feed-speed 
combinations  where  problems  were  encountered.  Trial  runs  should  be  made 
at  these  feed-speed  combinations.  This  process  is  continued  until  four 
feed-speed  combinations  (at  two  adjacent  feeds  and  two  adjacent  speeds) 
have  been  selected.  Each  feed-speed  combination  is  then  run  for  one 
shift  of  operation,  and  the  following  information  is  collected: 


1.  the  number  of  parts  produced, 

2.  the  number  of  tool  edges  used  for  each  tool  type,  and 

3.  the  production  time.l 


1 


Actual  machining,  work  piece  handling  and  tool  handling  time. 
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The  quality  of  the  parts  being  machined  should  always  be  observed,  and  if 
at  any  time  the  operator  feels  the  produced  parts  are  in  danger  of  not 
meeting  standards,  he  should  stop  taking  data  at  this  feed  and  speed 
combination. 

If  data  collection  is  discontinued  at  any  feed-speed  combination  during 
a production  run,  a decision  must  be  made  concerning  the  use  of  the  data 
already  collected  at  this  point.  If  the  number  of  parts  machined  is  sufficient 
to  give  accurate  measures  of  the  performance  indices,  the  results  can  be  used 
for  input  to  the  MACHOP  program,  even  though  a full  shift  of  data  has  not  been 
collected.  If  problems  of  workpiece  accuracy  are  encountered  or  if  it  becomes 
apparent  that  scrap  or  tool  breakage  costs  will  exceed  other  potential  gains 
during  a production  run  prior  to  obtaining  a sufficient  number  of  parts,  a new  set 
of  four  feed-speed  combinations  should  be  selected.  After  collecting  data  at 
each  of  the  four  combinations  and  recording  this  information  on  the  forms 
provided,  the  data  is  summarized  and  submitted  to  the  MACHOP  program. 

(Sample  data  collection  forms  are  given  in  Appendix  C,  and  a sample  summary 
form  is  given  in  Figure  G.  1 of  Appendix  G.) 

4,2  Analysis  of  the  MACHOP  Output 

Sample  MACHOP  outputs  appear  in  Figures  5.1  - 5.3.  Each  output 
gives  the  following  information  on  separate  pages: 

a.  Input  information  such  as  feed-speed  limits,  labor  and  overhead 
rates,  tool  costs,  and  prior  estimates  of  standard  deviations 

of  the  response  variables  (if  any) . 

b.  The  specified  feed-speed  environment  represented  "graphically." 

o.  The  cumulative  input  data  and  the  computed  responses  (cost 
per  piece  and  production  rate) . 

d.  A table  of  the  computed  responses  for  this  phase  and  cycle. 

e.  Evolutionary  operation  calculations  for  both  the  cost  per  piece 
and  production  rate  responses. 

f . A recommended  set  of  operating  conditions  for  the  next  experiment 
(which  may  be  Identical  to  the  current  set) , based  on  the  cost 
per  piece  analysis. 

g.  A graphical  representation  of  the  operating  conditions  recommended 
in  i.  . 

h.  The  recommended  set  of  operating  conditions  for  the  next 
experiment  (which  may  be  identical  to  the  current  set) , based 
on  the  production  rate  analysis. 
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i.  A graphical  representation  of  the  operating  conditions  recommended 
in  h. 

j.  The  results  of  the  regression  analysis  and  the  accompanying 
analysis  of  variance  table  for  the  cost  per  piece  analysis. 

k.  The  results  of  the  regression  analysis  and  the  accompanying 
analysis  of  variance  table  for  the  production  rate  analysis. 

t.  The  predicted  response  surfaces  for  cost/piece  and  production 
rate  calculated  from  the  respective  prediction  equations. 

The  output  information  in  a,  b,  c and  d should  be  reviewed  for  detection 
of  any  input  errors.  If  any  recording  and/or  transcription  errors  are  noted, 
the  input  should  be  corrected  and  resubmitted  to  the  MAC  HOP  program.  If 
it  is  determined  that  any  of  the  data  were  collected  under  abnormal  operating 
conditions,  new  data  should  be  collected  at  these  points.  The  results 
should  then  be  resubmitted  to  the  MAC  HOP  program. 

If  the  output  information  in  a,  b,  c and  d appears  acceptable,  the 
remainder  of  the  output  e-l  should  be  examined.  Parts  f,  g,  h,  i,  and  £ 
are  of  particular  interest  to  production  personnel.  Analysis  of  these  results 
should  allow  a decision  to  be  made  concerning  operating  conditions  for  the 
next  set  of  observations.  (The  other  parts  (e,  j,and  k)  give  additional  details 
concerning  the  computations.)  If  the  results  and  recommendations  seem  to 
be  reasonable,  the  data  should  be  collected  at  the  set  of  points  suggested 
by  the  MACHOP  program.  After  the  data  have  been  collected,  the  above 
procedure  is  again  followed.  This  process  continues  until  MACHOP  committee 
decides  that  sufficient  data  have  been  collected  in  order  to  determine  optimal 
operating  conditions. 

The  MACHOP  recommendations  should  be  treated  as  possible  courses 
of  action.  The  final  decision  rests  with  the  MACHOP  committee.  If  the 
committee  decides  to  collect  additional  data,  the  following  rules  must  be 
observed.  The  new  set  of  observations  must  be  taken  either  at  the  same 
four  points  (another  cycle  in  the  current  phase)  or  at  another  rectangular 
pattern  of  adjacent  points  (a  new  phase).  (A  more  complete  discussion  of 
phases  and  cycles  was  given  in  Section  3.2.1.)  We  suggest  the  movement 
be  restricted  to  feed-speed  combinations  immediately  adjacent  to  the  previous 
feed-speed  combinations.  Either  speed  or  feed  or  both  may  be  varied  at  one 
time. 

4 . 3 Scope  of  Application 

The  MACHOP  routines  can  be  used  for  single  operation  as  well  as 
multiple-operation  (numerically  controlled)  processes. 
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For  multiple-operation  (non-numerically  controlled)  machine  tools,  data 
must  be  collected  on  all  the  cutting  tools.  In  addition,  one  operation  of 
the  machine  must  be  selected  for  the  initial  study,  since  only  one  set  of  feeds 
and  speeds  can  be  varied  at  a time.  The  selectionoof  this  particular  opera- 
tion is  important.  In  general,  it  should  be  the  operation  which  has  the 
greatest  potential  for  improvement.  This  will  often  be  the  operation  requiring 
the  longest  time.  After  this  operation  has  been  optimized,  another  operation 
can  be  selected  for  study.  Sample  data  collection  forms  are  given  in 
Appendix  C . 

For  numerically  controlled  processes  a different  problem  is  encountered. 

In  this  case  the  feeds  and  speeds  for  the  operations  are  pre-programmed. 
However,  the  feeds  and  speeds  can  be  adjusted  by  overrides.  The  optimiza- 
tion program  is  conducted  through  the  use  of  these  overrides  in  order  to 
avoid  reprogramming.  Again  a rectangular  pattern  is  used,  now  in  terms 
of  percentage  feed  and  speed  overrides.  Increments  of  five  percent  seem 
reasonable.  If  analysis  of  the  data  collected  indicates  that  either  the 
speed  6r  feed  increment  should  be  adjusted  beyond  its  presently  programmed 
limits,  consideration  should  be  given  to  reprogramming  the  tape.  At  this 
point  the  relationship  of  the  various  speeds  and  feeds  should  also  be 
discussed.  Whenever  the  tape  is  changed,  the  MACHOP  program  must  be 
restarted.  That  is,  all  previous  data  should  be  removed,  and  Phase  1 
should  begin  again.  The  only  information  which  can  be  used  from  the 
previous  output  are  the  estimates  of  the  standard  deviations  and  the  costs. 

The  estimates  of  the  standard  deviations  from  the  last  output  should  be 
submitted  as  prior  estimates  of  the  standard  deviations.  For  numerically 
controlled  processes,  data  is  recorded  on  the  same  form  as  for  multiple- 
operation  processes.  However  instead  of  recording  only  feeds  and  speeds, 
the  feed  and  speed  overrides  are,  also  recorded. 

4.4  Data  Handling  System 

During  the  initial  set-up  of  the  problem  (phase  1,  cycle  1),  the 
necessary  data  such  as  speed  and  feed  environments  must  be  recorded. 

Once  this  is  accomplished,  all  succeeding  cycles  use  the  same  input 
with  the  addition  of  a new  set  of  observations  and  punched  output  from 
the  previous  cycle. 

Once  the  initial  information  is  given,  less  than  five  minutes  per 
cycle  will  be  required  for  data  recording  and  preparation.  Consequently, 
no  elaborate  data  handling  system  is  considered  necessary  for  this  process. 

A volume  of  four  new  cards  per  cycle  is  not  sufficient  to  justify  elaborate 
equipment.  The  most  effective  method  for  accomplishing  this  is  for  the 
individual  directly  interested  with  the  process  (ultimately  the  machinist) 
to  be  responsible  for  preparing  the  data  sheets  for  forwarding  to  the  data 
processing  unit  for  key  punching. 
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5 . MACHOP  Analysis  of  RIA  Data 


In  this  section  the  results  of  the  MACHOP  analysis  of  the  turning 
operation  (recoil  cylinder,  Part  No.  10895646)  on  a Monarch  lathe  are  given. 
The  details  of  this  particular  operation  are  given  in  Table  2.5.  A summary 
of  the  data  collected  is  given  in  Table  5.1.  For  this  machining  operation 
the  responsesvariables  (cost  per  piece  and  production  rate)  are  calculated 
from  the  production  data  as  explained  in  Section  3.1.  In  addition  to  the 
input  information  required  by  the  MACHOP  program,  a prior  estimate  of  the 
standard  deviation  of  the  cost  per  piece  of  0.60  dollars  was  specified,  based 
on  observations  made  during  the  initial  phases  of  the  study .1  No  estimate 
was  made  of  the  standard  deviation  of  the  production  rate. 

The  data  and  output  information  for  phase  1,  cycle  1 are  given  in 
Figures  5.1a  - 5.1-t.  The  cost  analysis  (Figure  5.  If)  indicates  that  the 
factor  feed  was  significant.  Since  no  prior  estimate  of  the  production  rate 
standard  deviation  was  supplied,  the  EVOP  analysis  indicates  none  of  the 
factors  was  significant.  The  response  surface-regression  analysis 
(Figure  8.  l£^_was  in  general  agreement  with  the  EVOP  analysis,  so  phase 

2 was  initiated.  The  feeds  and  speeds  recommended  by  the  EVOP  cost  per 
piece  analysis  (Figure  5.  If)  were  used  for  phase  2,  cycle  1. 

The  MACHOP  output  for  phase  2,  cycle  1 is  given  in  Figures  5.2a  - 
5.21.  Based  on  the  new  data,  a three  variable  regression  equation  was 
fitted.  (See  Section  3.3  for  details.)  Neither  of  the  resulting  EVOP 
analyses  (Figures  5 . 2 f and  5.2h)  indicated  that  any  factor  was  significant. 
Consequently,  another  cycle  of  data  was  taken  for  phase  2. 

The  MACHOP  output  for  phase  2,  cycle  2 appears  in  Figures  5.3a  - 
5.3-t.  The  EVOP  analysis  of  the  cost  per  piece  and  production  rate  results 
(Figures  5.3f  and  5.3h)  again  indicates  that  none  of  the  effects  are  signifi- 
cant. However,  the  predictions  of  cost  per  piece  and  production  rate  in 
response  surface-regression  (Figure  5.3£)  indicate  the  same  direction 
toward  the  optimum.  On  the  basis  of  these  results,  we  could  take  another 
cycle  of  data  for  phase  2 or  begin  phase  3.  We  recommend  beginning  phase 

3 of  the  study  at  either  of  the  following  sets  of  feed-speed  combinations: 


Speed 

Feed 

Speed 

Feed 

166 

0.0210 

166 

0.0240 

166 

0.0240 

or 

166 

0.0270 

192 

0.0210 

192 

0.0240 

192 

0.0240 

192 

0.0270 

However,  before  beginning  this  phase,  a trial  run  should  be  made  at  each 
of  the  new  combinations  to  insure  acceptable  machining  results. 

Comparison  of  the  predicted  costs  in  Figure5.3£  gives  an  indication  of 
the  potential  for  cost  reduction  through  optional  selection  of  feeds  and  speeds. 

*Note  that  an  estimate  of  the  standard  deviation  is  not  required. 
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Table  5 . 1 Summary  of  RIA  Data 


No.  of  No.  of 


Date 

Speed 

Feed 

Parts 

Time 

Tool  Edges 

1)  6/12/73 

192 

0.0187 

18 

410 

18 

2)  6/13/73 

220 

0.0210 

20 

440 

44 

3)  6/27/73 

255 

0.0187 

11 

267 

29 

4)  6/28/73 

192 

0.0210 

14 

377 

18 

5)  6/29/73 

220 

0.0189 

15 

393 

15 

6)  7/5/73 

192 

0.0168 

14 

434 

17 

7)  7/10/73 

220 

0.0168 

13 

396 

23 

8)  7/17/73 

255 

0.0168 

10 

217 

17 

9)  7/18/73 

255 

0.0210 

14 

357 

31 

10)  8/3/73 

192 

0.0187 

14 

374 

16 

11)  8/6/73 

220 

0.0187 

15 

394 

17 

12)  8/7/73 

220 

0.0187 

9 

228 

15 

13)  8/16/73 

192 

0.0187 

8 

214 

16 

14)  8/17/73 

220 

0.0187 

13 

381 

20 

Comments 

None 

Inserts  break  in  half  or  chip.  The  out  of  round  condition 
seems  to  contribute  to  insert  breakage  as  I ran  2 pieces 
on  one  insert  because  they  were  fairly  straight  and  had 
trouble  running  one  piece  with  2 inserts  when  the  piece 
was  extremely  out  of  round. 

None 

Chips  are  still  too  tight.  Tips  of  inserts  are  chipping 
at  times  depending  on  out  of  round  condition  of  piece. 
Good  chips.  Finish  real  good  considering  the  condition 
of  the  pieces. 

None 

Tips  of  inserts  are  chipping  easily  on  irregular  pieces. 
Finish  is  poor  as  insert  breaks  down  at  about  3/4  of  the 
way  down  the  piece.  Chips  are  small  and  tight. 

None 

Inserts  break  down  about  1/4  way  across  cut.  Finish 
is  horrible. 

None 

Inserts  break  down  1/3  of  the  way  through  cut.  Finish 
is  very  rough.  (1  insert  broke  on  end  of  piece  so  was 
able  to  use  replaced  edge  on  next  piece.) 

None 

None 

Finish  was  real  good  even  for  first  time  pieces  done 
with  one  insert  edge,  but  on  third  piece,  insert  had 
a tendency  to  break  down  1/3  of  the  way  across. 
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Figure  5.1  MACHOP  Output  for  Phase  1,  Cycle  1 
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INPUT  VALUES  ANC  COMPUT ED  RESPONSES 


SPEED 

FEED 

PAR)  S 

TIME 

(1) 

(11) 

(2) 
( 12) 

(3) 

(13) 

TOOL-EDGES 
(4)  (5)  (6) 

(14)  (15)  (16) 

( 7) 
( 17) 

(8) 
( 18) 

(9) 
( 19) 

( 10) 
(20) 

COST 
( S/PIECE) 

PROD  RATE 
(PIECES/3IN) 

19  2* 

0. C168 

14. 

434  • 

17. 

9.3  l 

0.0323 

220. 

C. 0187 

15. 

393. 

15. 

9.28 

0.0382 

22  0. 

0.0168 

13. 

396. 
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C. 0 187 

18. 

410. 
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Figure  5.1  (continued) 
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Figure  5.1  (continued) 
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220. 

192. 

192. 
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Figure  5 . 1 (continued) 
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Figure  5 . 1 (continued) 
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Figure  5 . 1 (continued) 
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Figure  5 . 1 (continued) 
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MULTIPLE  REGRESSION  ON 


2 VARIABLES  U I T H 


<i  OBSERVATIONS 
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Figure  5 . 1 (continued) 
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Figure  5 . 1 (continued) 
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Figure  5.2  MACHOP  Output  for  Phase  2,  Cycle  1 
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Figure  5.2  (continued) 
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Figure  5.2  (continued) 
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Figure  5.2  (continued) 
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Figure  5.2  (continued) 


COST  ANALYSIS 


THE  EFFECT  OF  SPtEO  IS 
THE  EFFECT  OF  FEED  IS 

THE  EFFECT  UF  THE  INTERACTION  OF  SPEED  AND  FEEO  IS 

NO  PARAMETERS  ARE  SIGNIFICANT 

THE  PREVICLS  SETTINGS  SHCUIC  BE  RE-EXAVINED 


-0.0018 

NUT 

SIGNIFICANT 

-0.5392 

NOT 

SIGNIFICANT 

-0.2928 

NUT 

SIGNIFICANT 

cn 


RECOMMEND  ADDITIONAL  OBSERVATIONS  BE  TAKEN  AT  POINTS: 


4)  SPEED  = 
FEED  = 


192. 

0.0210 


2)  SPEED  = 
FEED  * 


220. 

0.0210 


1)  SPEED 
FEEO  = 


192, 

C.0187 


3)  SPEED  = 
FE  EC  * 


220. 

0.0187 


f. 


Figure  5 . 2 (continued) 
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Figure  5.2  (continued) 
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Figure  5 . 2 (continued) 
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Figure  5.2  (continued) 
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Figure  5.2  (continued) 
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Figure  5.2  (continued) 
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P.R.  STD.  DEV.  EST.  IS  NOT  SPECIFIED 


NUMBER  OF  TOOLS  IS 

I 

SPEC  I FIEO  TOOL  COST/EDGE 

0.4200 

NUMBER  OF  SPEEDS  IS 

10 

NUMBER  OF  FEEOS  IS 

11 

a . 

Figure  5.3  MACHOP  Output  for  Phase  2,  Cycle  2 
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Figure  5.3  (continued) 
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Figure  5 . 3 (continued) 
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Figure  5.3  (continued) 


EVOLUTIONARY  OPERATION  ANALYSIS 
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Figure  5 . 3 (continued) 
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COST  ANALYSIS 
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Figure  5 . 3 (continued) 
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Figure  5 . 3 (continued) 


PRODUCTION  RATE  ANALYSIS 


THE  EFFECT  OF  SPEED  IS  -0.0000 

THE  EFFECT  OF  FEED  IS  -0.0026 

THE  EFFECT  OF  THE  INTERACTION  OF  SPEEO  AND  FEEO  IS  -0.0013 

NO  PARAMETERS  ARE  SIGNIFICANT 
THE  PREVIOUS  SETTINGS  SHOULD  BE  RE-EXAMINED 


NOT  SIGNIFICANT 
NOT  SIGNIFICANT 
NOT  SIGNIFICANT 


RECOMMEND  A00ITI3NAL  OBSERVATIONS  BE  TAKEN  AT  POINTS: 


4)  SPEED  - 

192. 

2)  SPEEO  * 

220. 

FEED  * 

0.0210 

FEEO  * 

0.0210 

l)  SPEED  » 
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3)  SPEED  - 

220. 
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0.3187 

FEED  * 

0.0187 

h. 

Figure  5 . 3 (continued) 
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Figure  5 . 3 (continued) 


MULTIPLE  REGRESSION  ON 


3 VARIABLES  WITH  12  OBSERVATIONS 


COST  EQUATION 

ANOVA 


SOURCE 

OF 

SS 

MS 

F 

TOTAL 

11 

0.16933578E 

02 

REGRESSION 

3 

0.15972290E 

01 

0.532*0967E  00 

0.27772433E  00 

RES  IDUAL 

8 

0.15336349E 

02 

0 * 191 70437E  01 

COEFFICIENT  OF  MULTIPLE  DETERMINATION  (R**2)  * 0.94 3231 58E-0 I 


COEFFICIENTS 

T VALUES 

STANDARD  DEV. 

BSUBO 

•0.72270166E  03 

BHAT  ( 

1 ) = 

0.  13664<t7<VSE-  03 

T ( 

l)=  0.45348996E 

00 

SD( 

1)«  0 .30 13 1812 E 
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2 ) *- 

0.178687088E  03 

T ( 
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00 
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2)-  0 • 40534497E 

03 
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0.333668671E  02 

T ( 

3>«  0.43842232E 

00 

SD  ( 

3)*  0.76106674E 

02 

THE  PREDICTION  EQUATION  IS 

-722*7017  ♦ 

136.6447*LN( SPEED)  ♦ 

-178 

.687X*LN(FEED) 

♦ 

33. 366 9*LN (SPEED) *LN( FEED) 
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Figure  5.3  (continued) 


MULTIPLE  REGRESSION  ON  3 VARIABLES  WITH 


12  OBSERVATIONS 


PRODUCTION  RATE  EQUATION 


ANOVA 


SOURCE 

OF 

SS 

MS 

F 

TOTAL 

il 

0.101274246-03 

REGRESSION 

3 

0.54460019E-04 

0.181533406-04 

0.114519316  01 

RESIDUAL 

3 

0. 126814226-03 

0.158517656-04 

COEFFICIENT  OF  MULTIPLE  DETERMINATION  (R**2)  » 0.30042887E  OC 
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6SU80 
BHAT  ( 

*-0.139120106  01 
1)*  0.28  304028  C6  00 

T ( 

11-  0.332433826 

00 

SD(  1)« 

0.866459016 

00 

8H  AT  ( 

2)*-0. 3692927366  00 

T ( 

2)— 0.316827606 

00 

SO ( 2)* 

0.116559516 

01 

BHAT  ( 

3)»  0.7433050876-01 

T ( 

3)*  0.339642056 

00 

SD(  3)* 

0.218849546 

00 

THE  PREDICTION  EQUATION  IS 

PI  - -1.3912  ♦ 0.2880*LN{ SPEED)  ♦ -0 . 3693*LN(  FEED)  0 .074  3*LN  (SPEED  > *LN<  FEED  J 

k. 

Figure  5 . 3 (continued) 


PREOI CTEO  COSTS 
(PREDICTED  PRODUCTION  RATES) 
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Figure  5 . 3 (continued) 


6t  Summary 


Optimization  of  machining  parameters,  such  as  cutting  speed  and  feed, 
is  possible,  in  relation  to  time  and  costs,  by  use  of  this  developed  sys- 
tem and  production  machining  data.  This  system,  MACHOP  (Machine  Optimiza- 
tion Program),  is  applicable  to  analysis  and  control  of  single-operation 
as  well  as  multiple-operation  (numerically-controlled)  machining.  Simple 
formats  are  used  to  gather  data  from,  and  return  computerized  results  to, 
production  personnel.  In  application,  two  variables  of  the  machining 
parameters  may  be  adjusted  independently  or  simultaneously,  either  manually 
by  the  machinist  or  by  the  numerical -control  programmer,  using  small-step 
increments  within  the  limits  of  the  machine  tool  or  the  workpiece  tolerances. 
Computation  of  results  from  data  collected  at  the  original  adjusted  machining 
parameters  provides  print-outs  of  four-point  performance  indices.  These 
print-outs  not  only  show  time  and  cost  performance  at  the  machining  para- 
meters from  which  the  data  for  computation  has  just  been  taken,  but  also 
give  production  settings  which  could  improve  performance.  Subsequently, 
machining  at  adjoining  machining  parameters  is  tested  and  analyzed  until 
the  print-outs  repeatedly  indicate  that  the  optimal  conditions  have  been 
reached.  Notably,  the  data  system  and  computer  program  are  simple,  and 
general,  enough  to  readily  allow  optimization  of  other  machining  para- 
meters such  as  cutting  fluids,  work  material  selection,  and  tool  material 
and  geometry,  as  well  as  speeds  and  feeds. 

The  development  of  this  system  and  related  computer  program  has  ad- 
vanced the  state-of-the-art  of  optimization  in  machining.  In  the  past, 
optimization  was  restricted  to  optimizing  an  expanded  Taylor  Tool-Life 
equation,  not  optimization  of  the  actual  production  machining  operation 
itself,  or  required  more  adjustments  of  variables  than  were  practicable 
in  production  operation. 


7.  Recommendations 


It  is  recorrmended  that  this  system  and  computer  prooram  be  applied 
to  analyze  and  control  all  major  machining,  and  similar  process,  opera- 
tions in  which  a large  number  of  parts  are  produced  in  a single  run  or 
repetitive  runs.  It  is  also  recommended  that  it  be  used  to  augment  or 
replace  time-study  methods  presently  used  to  establish  production 
standards;  and,  in  particular,  that  it  be  applied  to  all  new,  major 
machining  operations,  and  operations  where  new  parameters,  such  as  work 
material,  cutting  fluid,  or  tool  material  and  geometry  are  introduced. 
Furthermore,  it  is  strongly  recommended  that  it  be  used  to  provide 
computerized  print-outs  of  cutting  tool -life  and  costs  per  workpiece  for 
procurement  and  use  of  quality  tools  according  to  tool  performance  in- 
stead of  merely  tool  price.  Finally,  it  is  recommended  that  this  optimi- 
zation procedure  be  used  on  other  processes  in  which  operations  parameters 
may  be  adjusted  to  improve  control  of  time  and  costs. 
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Appendix  A.  Verification  of  Regression  Modules  of 
the  PIM  and  MAC  HOP  Programs 


This  appendix  contains  an  output  from  a stepwise  regression  program, 
which  was  used  to  check  the  accuracy  of  the  regression  portion  of  the  PIM 
program.  A similar  check  was  performed  on  the  MACHOP  program. 
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SUP  NUMBER  5 
VARIABLE  ENTERED  I 


MULTIPLE  R 0.6104 

STO.  ERROR  CF  EST.  0.C172 


ANALYSIS  OF  VARIANCE 

CF  SUM  OF  SCUARCS 
REGRESSION  5 0.001 

RESIDUAL  3 0.001 


MEAN  SQUARE  F RAT  IU 

0.000  0.356 

0.000 


VARIABLE 


VAR I AbLES  IN  EQUATION 
COEFFICIENT  STD.  ERROR  F TO  REMOVE 


VARIABLE 


VAP IABLES  NOT  IN  EQUATION 
PARTIAL  CORR.  TOLERANCE 


(CONSTANT 
SPEED  1 
F E ED  2 
SFEDSQ  4 
FEEDSQ  5 
SFEDFD  6 


0.90118  ) 

0.08087 

0.96511 

U .03633 

0.17757 

0.12478 


0.60953 

1.05571 

0.12451 

0.279C3 

0.12969 


0.0176  (2) 
C.7758  (2) 
0.0853  (2) 
C.4050  (9) 
C • 9257  ( 2) 


F-LEVEL  OR 


T CL  FRANCE  INSUFFICIENT  FOR  FURTHER 


COMPUTATION 


SUMMARY  TABLE 


STEP 

NUMBER 

VARIABLE 

ENTERED  REMOVED 

multiple 

R 

RSQ 

INCREASE 
IN  RSO 

F VALUE  TO 
ENTER  CP  REMOVE 

1 

FEEDSQ 

5 

0.1995 

0.0400 

0.0400 

0.2913 

2 

FEED 

2 

0.3584 

0.1284 

0.0885 

0.6091 

3 

SPEOSQ 

4 

0.4038 

0.1630 

0.0346 

0.2067 

4 

SPEDFD 

6 

0.6074 

0.3639 

0.2059 

1.3043 

5 

SPEED 

1 

0.6104 

0.3726 

0.0037 

0.0176 

F TO  ENTER 


NUMBER  OF  INDEPENDENT 
VARIABLES  INCLUDED 


1 

2 

3 

4 

5 


A, I Sample  Stepwise  Regression  Output 


Appendix  B.  Simulation  Program 


This  appendix  contains  a listing  of  the  simulation  program  and  a 
typical  output.  The  instructions  for  using  the  program  are  given  in  the 
program  as  comments.  This  simulation  program  was  developed  for 
preliminary  analysis  of  the  PIM  and  MACHOP  programs. 
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FORTRAN  IV  0 LEVEL  21 


MAIN 


DAT E * 73361 


17/31/li 


OOCi 

IMPLICIT  R E AL  ♦ 8 (A-F,0-Z) 

OOC2 

L>  1 MENS  IC  N P(30),T(3C),1P(30) 

00C3 

DIMENSION  F(3U),  V( 30 ) , PR ( 30) , CU ( 30 ) , P I ( 30) 

ooc4 

CO  1 MON  INT 

C 

c 

NPT  - IS  THE  NUMBER  OF  POINTS  AT  WHICH  THE  l)ATA  IS  TO  RE  GENERATED 

c 

c 

r 

0 - IS  THE  PERFORMANCE  INDEX  COEFFICIENT 

L 

c 

INT  - IS  THE  kANDCM  SEED*  A SEVEN  DIGIT  INTEGER 

00C5 

c 

READ  (5,12)  NPT , 0, 1 NT 

00C6 

12 

FORMAT  ( I 10tFl0.2,I  10) 

c 

c 

RLO  - IS  THE  LABOR  ANC  OVERHEAD  IN  i/MlN 

c 

TLC  - IS  THE  TOOL  COST  IN  R/EDGE 

OOC  7 

c 

READ  (5,20)  RLO, TIC 

OOC8 

2u 

FORMAT  (2F1U.4J 

OOCS 

WRITE  (6,30)  RLO, TLC 

001C 

3u 

FORMAT  ( )0X,»  Si  MUUTI  CN  OF  PRODUCTION  DAT  A*///25X,  'LAPOF  C OVERMEA 
ID  CJST  * **,F5.2,*  /MIN*//3CX,  ' 1 0 JL  CCST  * V»,F5.2,*  /F  DGfc*  //  ) 

OOll 

READ  (5,10)  ( (V(I  ),F(  I)  ),  1*1, NPT) 

0012 

10 

FORMAT  ( 2 F 1 0. 4 ) 

0013 

DC  15  J = l , NPT 

OOM 

CALl.  SIM(J,F,V,NCT,KPT,TPD) 

0015 

P(J)  - KPT 

0016 

T(J)  * TPO 

0017 

TP( J)  « NCT 

0018 

15 

CON  T I NUE 

001S 

CALL  PtRINU(P,T,TP,CJ,PR,PI, 3 , RLO, T LC, NPT ) 

002C 

wRI  rc  (6,24) 

0021 

24 

FORMAT  ( • 1 1 , 10X , ' SPEED  FEED  *,•  PARTS  TOOl 

l TIME  COST/  •,*  PROD.  •,*  PI-  '/ 

213A, • (RP  1) »,5X,  • ( 1PR) \ 17X, • CHANGES  1 , 14X, • PI ECE • , 6X, ' RATE*, 7X , • I NO 
3 E X ' /) 

0022 

wRITE  (6,25)  ( ( V( I) ,F( 1 ) , P(  I ),TP (I ) ,T(I ) ,CU( I ),FR(I ), PI ( I )) ,I«1, 
1NP1  ) 

0023 

25 

FORMAT  (IOX,F10.1  ,F  10  .4,  F ID  . 1 , F 1 0. 1 , F 1 0 . 1 , F l 0 .2, 2F12 . 4/ ) 

00  24 

RETURN 

0025 

END 

FORTRAN  IV  G LEVEL 


21 


RAND 


DATE  - 73361 


17/31 / 1 i 


OOCI 

SUBROUTINE  RANDCINT.IREALI 

00C2 

KEAL^B  DREAl,  oint,  OKEEP 

c 

GUMfcN f NO  2 IN  LIST  IS  A REAL  IN  THE  CALL  STATEMENT,  OUT  AN  INTEGER 

c 

IN  THE  SUBROUTINE 

Oo  C 3 

DINT  « I NT 

OOC  A 

DREAL  - HINT  ♦ 16PC7. 

00C5 

KEtP  « OREAL  / 21A7A836A7. 

OOC  t 

DKEEF  - KEEP 

00C7 

INV  * DREAl  - OKfcEP  * 21A7A836A7. 

ooce 

IRbAt  * (1NT/129)  ♦ 10/37AIE2A 

occs 

RETURN 

001C 

END 

FORTFAN  IV  G LEVEl  21 


* UMNOM 


DATE  » 73361  17/31/11 


03  C 1 
00C2 
00C3 
DOC* 
00C5 
00C6 
UUC7 

occe 


FUNCTION  UN INUM  C t FP  , S IGMA 1 
IMPLICIT  Rb  Al*3  ( A-F,  C-Z ) 

COl^ON  I NT 

CALL  RAMH  INT.REALI 

XNURM  «(*EAL**u.  1926966100-  (1  -REAL  M *0.1 9 2i>'?o5  l DO  1/3.  1306A5600 

UMlNUM  ■ EFR*SIGFA*XNCRM 

RETURN 

CND 


FORTRAN  IV  G LEVEL  21 


PER INO  OATE  « 73361  17/31/11 


U0C1 
00C2 
00  C 3 
OOCA 

00C5 
D0C6 
00  C 7 
OJC  6 
OOC9 


SUBROUTINE  PE RI NO (P»T|TP#CU»PRfPIi3»PL0*TLC»NPTl 

IMPLICIT  REAL*aU-M,0-/ J 

OlMfcNSlLJ  PC30I ,T (3CIVTF(30 I 

DIMENSION  F (301  » VC  301 , PR  C 30  , CU1  331  , PI  C 3 01 

DO  luO  1*1, NPT 
Pk(  I )*P( I I / T t II 

CUC  KI«|kLC«TI  1 1 *T  LC *TP(  II  )/ PC  II 
100  PIC  1 I « C l.-OI^PfvCI  I ♦ CC/CU  ( 1 1 I 
RETURN 


001C 


END 
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FORTRAN  IV  G LEVEL  21 


SIM 


OATE  - 73361 


17/31/11 


OOC1 
0002 
00C3 
OOC  A 


0UC5 
uuC6 
OOC  7 

ooce 

oocs 


00  1C 
0011 


0012 

0013 

00  1 A 

0015 

0016 
0017 
0010 
OOIS 
002C 
0021 

0022 
0023 
00  2A 
0025 
00  26 

0027 

0028 
0J2S 

CO  30 

0031 

0032 

0033 
003A 

0035 

0036 

0037 

0038 


003S 
00  AC 
00  A 1 
OOA  2 
00  A3 
00  A A 


SUUK0U1  ME  SlM(If  F,  V,NTC,  NTPT,TTPO) 

IMPLICIT  RbAL*8  (A-H.0-2) 

OlMtNSILl  F(16),V(16) 

COM  PON  CUT 

c 

c 

c in  this  subroutine  all  the  parameters  of  the  production  ppocfss 

C APE  INITIALISED  rather  than  to  be  read  by  a read  statement 

c 
c 

C 1 AY  LOR • S TOOL  LIFE  CONSTANTS 

c 

ALP  * O.J 
BETA  « C • 2 
CONST  • AOO.O 

C XL  — LENGTH  UF  WORKPIECE 
XL  * A 7 • 5 

C U — CIAMETtR  OF  WCRKPIECE 
0 * fa  • 5 

c tct  --  tuul  changing  time 

TCT  « 1. 

C TM  — MACHINING  TIME  / PIECE 
TM  * XL/ ( V ( I ) ♦F ( I n 
C 1 H — HANDLING  TIME  / PIECE 
TH  * 15.0 

C SHIFT  — UTaL  PkODGCTICN  TIME  / SHIFT 
SHIFT  • A2C.0 

C ERR  - IS  Tml  KANOCM  ERROR 

ERR  - 0*20 

XV  * 3. 1A15^U*V(I) /12. 

SIG  = ( Dl  EG  (CLNST  ) - PETA*  Ui.OG  ( F ( I ) )-DLOG(  XV  I ) / Al  P 
ICULNT  • 0 
TIME  - 0. 

PARTS  - 0. 

WRITE  (6,15)  I • V ( I ) * F ( I » 

15  FORMAT  (' 1' «10X, "FOR  DAY*, 12  /10X,*  SPEED- • , F 6. 1 , 5 X, «F EE 0« ' , F7 . A/ 
l/) 

WRITE  (6,301  TM 

30  FCRM/T  ( 1 MACHMNG  TIME  FOR  THIS  SPEED  C FEEO  • , F 7 .2,  • MIN»//| 
10  ICOUNT  » I COUNT  ♦ 1 
Y « UMN3M(LPR,S1G) 

TLI FE  - OEXPCSIGtY) 

PAkTSl  » TLIFE/TM 

WRITE  (6,20)  ICOUNT, TLIFE,PARTS1 

20  F Ok  MAT  (•  TOOL  * * , I 3, 5X , • L I F E* ' , F5 . 1 , • MIN • , 5 X , • # OF  PARTS  -‘.F5 .2/ 
1/  ) 

1IME1  * PAR7Sl*(TM*TH)fTCT 

PAR  TS  - PARTS+PART SI 

TIME  » T I ME ♦TlMEl 

IF  (TIME  .LC.  SHIFT)  GO  TO  10 

DlFF  - TIME-SHIFT 

PAR VMS  - CIFP/(TM»TH) 

PAR  IS  - PARTS-PARTMS 

WRITE  (6,25)1  .PARTS, ICOUNT,  SHIFT 

25  FORMAT  (•  SO  FOR  DA Y • , I 3, / 1 OX, F 5 . 1 , • P ARTS  ARE  PRODUC ED* /10A, 13, 

1 • TUCLS  ARE  LHANGEC  CU^ ING 1 / 10X , Fo . 1 , 1 MINUTES  OF  PRODUCTION  TIME 
1 • / ) 

PARTS  * PARTS  ♦ 0.25 

NTPT  * PARTS 

TTPD  * SHIFT 

nTC  * ICJGNT 

RETURN 

END 


SPEED*  192.0 


FEED*  0.0187 


4 


PACrilN  ING  TIME  FOR  THIS  SPEED  L 


TOCL  n 1 LIFE*  3.1MIN 


TOCL  * 2 


LIFE*  10.1M1N 


TOCL  * 3 


LIFE*  19.9MIN 


7CCL  n 4 LIFE*  7.7MIN 


TOOL  U 5 LIFE*  3.5  PIN 


TOCL  * 6 LIFE*  2.4MIN 


TOCL  t 7 


LIFE*  27.2MIN 


TOOL  A 8 LIFE*  2.0PIN 


TOCL  * 9 LIFE*  1.7MIN 


TOOL  A 10  LIFE*  10.7MIN 


TOCL  A 11  LIFE*  38.9MIN 


TOOL  A 12  LIFE*  14.1 


TOCL  A 13  LIFE*  9.2MN 


TOOL  A 14 


LIFE*  26.1PIN 


TOOL  A 15 


LIFE*  13.4PIN 


FEED  13 

.23 

OF 

PARTS 

* C 

OF 

PARTS 

* C 

OF 

PARTS 

* 1 

OF 

PARTS 

* c 

OF 

PARTS 

* c 

OF 

PARTS 

* c 

OF 

PARTS 

* 2 

OF 

PARTS 

* c 

OF 

PARTS 

* c 

OF 

PARTS 

* c 

OF 

PAFTS 

* 2 

OF 

PARTS 

* 1 

OF 

PARTS 

* C 

OF 

PARTS 

* 1 

GF 

PARTS 

* 1 

14.3PAKTS  ARE  PRODUCED 
15  T COL S ARE  CHANGED  DURING 
420. C MINUTES  CF  FROCUCTION  TIRE 


MIN 

.23 

.77 

.51 

.58 

.26 

.18 

.05 

.15 

.13 

.81 

.94 

.07 

.69 

.98 

.01 
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Appendix  C.  Data  Collection  Forms 


Two  types  of  data  collection  forms  are  given,  one  for  single  tool 
operations  (Figure  C.  1)  and  the  other  for  multiple  tool  operations 
(Figure  C.2), 
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SINGLE  TOOL  MACHINE  OPTIMIZATION  STUDY  FORM 


Machine 

Date 

Stock  Lot  ID 

Speed 

Operator 

Shift 

Tool  Material  iD 

Feed 

Production  Production  Subtotals 

Startup  Time  Stop  Time  (minutes) 


Total  Production  Time  (minutes) 

Number  of  Pieces  Produced 
Number  of  Insert  Edges  Used 
Comments: 

Figure  C . 1 


Total 
Total 
Cost/edge 
Overhead  and  Labor  Cost 


(Tally) 

(Tally) 
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MULTIPLE  TOOL  MACHINE  OPTIMIZATION  STUDY  FORM 


Sheet  No. 

Machine 

Speed 

Feed 

Operator 

Date 

Shift 

Stock  Lot  ID 

Insert  1 ID 

Insert  2 ID 

Insert  3 if) 

Production  Production  Subtotals 

Startup  Time  Stop  Time  (minutes) 


Total  Production  Time  (minutes) 


Number  of  Pieces  Produced 


Total 


Total 


Total 


Number  of  Insert  3 Edges  Used 


Comments: 


Total 
Cost/edge 
Overhead  and  Labor  Cost 
Figure  C.2 
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Appendix  D.  Analysis  of  the  PIM  Design  Module 


Following  the  determination  of  the  optimum  feed  and  speed  based  on 
a given  set  of  experimental  observations  the  next  set  of  feed  and  speed 
combinations  are  chosen  in  the  following  manner: 

1.  In  the  event  that  the  optimal  point  is  not  on  a boundary,  nine 

or  fewer  points  are  selected  as  follows.  Initially  a 3^  design  is 
selected  with  the  previous  optimal  point  as  the  center  point, 
the  distance  between  speeds  is  L 4 = number  of  usable  speeds/4  and 
the  distance  between  feeds  is  M4  = number  of  usable  feeds/4.  Each 
point  is  checked  with  the  boundary  conditions  to  insure  that  it  is 
feasible.  If  it  is  outside  of  the  boundary  condition,  a replacement 
point  is  selected  by  moving  the  point  in  one  step  horizontally  or 
vertically  depending  on  the  boundary  condition.  The  procedure  is 
repeated  until  a feasible  point  is  selected.  If  this  point  does  not 
duplicate  another  point  already  selected  it  is  retained,  otherwise 
it  is  discarded.  Thus,  fewer  than  nine  points  are  possible. 

2.  In  the  event  that  the  optimal  point  is  on  a boundary,  five  or 
fewer  points  are  selected  along  the  boundary  with  the  previous 
optimal  point  in  the  middle.  As  in  the  previous  case,  each  point 
picked  is  checked  for  feasibility  and  uniqueness.  Ifrihe  initial 
points  selected  are  not  feasible,  they  are  moved  in  until  they  are 
feasible.  In  the  process  if  some  of  the  points  duplicate  other 
points  already  selected,  they  are  discarded  thus  resulting  in  less 
than  five  points. 
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Appendix  E.  Carboloy  Systems  Computerized 
Machinability  Program 

This  appendix  contains  the  results  obtained  from  General  Electric's 
Carboloy  computerized  Machinability  Program  for  the  turning  operation  on  the 
recoil  cylinder.  (See  Table  2.5  for  the  job  description.) 

Access  to  this  program  may  be  made  by  acquisition  of  a G.E.  terminal 
installation  in  the  user's  plant  on  a contractual  basis  or  otherwise  by  th 
sending  the  required  information  to  the  company  for  processing. 

The  program  is  based  on  historical  machinability  data  and  Taylor's 
tool  life  equation.  If  the  tool  material  is  entered  into  the  program,  the 
recommended  feed  and  cutting  speed  for  minimum  cost  and  for  maximum 
production  rate  will  appear  in  the  output  data.  If  the  feed  and  speed  are 
entered , the  recommended  tool  material  will  be  indicated  in  the  output  data . 

This  approach  has  its  limitations,  for  there  are  many  variables  which 
effect  removal  performance  that  cannot  be  accounted  for  in  a generalized 
program.  These  variables  include  the  condition  and  rigidity  of  the  machine 
tool,  type  of  tool  holder,  experience  and  ability  of  the  machine  operator, 
and  the  tool  setting  practice.  Certainly,  the  program  cannot  be  expected 
to  predict  precisely  the  optimum  machining  parameters.  It  may,  however, 
serve  as  a good  starting  point  to  select  the  initial  machine  settings  and/or 
tool  material. 

A copy  of  the  program's  output  is  given  in  Figure  E.  1.  The  Carboloy 
tool  material,  which  was  entered  into  the  program  by  a G.  E.  representative 
and  which  was  considered  to  be  equivalent  to  the  VR/Wesson,  C5W  Titanium 
Coated  Carbide,  was  Carboloy  Grade  350.  The  second  choice  was  Grade  78. 
The  equivalence  of  the  entered  grade  to  the  one  being  used  is  questionable. 
The  program  recommended  the  following  machining  conditions: 

Feed  - 0.014  ipr, 

Speed  for  Minimum  Cost  - 272  rpm,  and 

Speed  for  Maximum  Production  Rate  - 358  rpm. 
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CUTTING  FLUID  0N  C ARB0LOY  MAY  CAUSE  TOOL  FAILURE 
DUE  T0  THERMAL  CRACKS.  CHECK  T0  MAKE  SURE  T00LS  ARE 
WEARING  OUT,  NOT  CHIPPING  0R  BREAKING 

IDENTIFICATION  AISI  4140 


STARTING  SURFACE  C0NDITI0N 
FINAL  SURFACE  FINISH 
BRINELL  HARDNESS 
PART  T OLERANCE 
COOLANT 
MACHINE  T00L 

TYPE  LATHE 

OPERATION  AXIAL 

M0TOR  H0RSEP3WER 
SPINDLE  SPEED  LIMIT 
CUTTING  T00L 

TOOL  MATERIAL 
GRADE 

2ND.  CHOICE  GRADE 
SIDE  CUTTING  EDGE  ANGLE 
NOSE  RADIUS 
OUTPUT  DATA 
SET  UP 

STARTING  DIA. 

FEED 

DEPTH  0F  CUT 
SPEED  MIN.  COST 
SPEED-MAX.  PROD. 

TECHNICAL 

WEAR  LAND 

TOOL  LIFE-MIN.  COST 
SPEED-MI N . COST 
H.P.-MIN.  COST 
TOOL  LIFE-MAX.  PROD. 

SPEED-MAX.  PROD. 

H. P. -MAX.  PROD. 

IF  CHIPPING  OF  THE  CUTTING  EDGE  OCCURS  USE 


IN 


HTR 

250  MU.  IN. 
252 
.015 
EMU 


TURN 

50 

380 


C AR80L0Y 
350 


78 

30  DEG. 
0.047  IN. 


8.500  IN. 

.014  IN. 
0.125  IN. 
272  RPM 
358  RPM 

.030  IN. 


9 

60  6 

1 1 

3 

798 
1 4 

SECOND 


/RE\ 


MIN. 

FPM 

MIN. 

FPM 

CHOICE 


Figure  E.  1 Sample  Carboloy  Computerized 
Machinability  Program  Output 


•AA  ’ 


GRADE. 
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Appendix  F.  MAC  HOP  Listing 


FORTRAN  IV  C LEVEL  21  MAIN  OATE  • 713*2  14/0*207  FACE  0001 


C 

HA 

IN 

10 

C 

MA 

IN 

20 

C 

MACHO P 

MA 

IN 

30 

C 

MA 

IN 

40 

C 

(MACHINE  OPTIMIZATION) 

MA 

IN 

50 

C 

MA 

IN 

60 

C 

MA 

IN 

TO 

C 

THE  MACHOP  PROGRAM  UTILIZES  EVOLUTIONARY  OPERATION  (EvCPJ 

MA 

IN 

80 

c 

AND  regression  ANALYSIS  tc  examine  the  RESPONSE  SURFACE  OF 

MA 

IN 

90 

c 

MACHINING  OPERATIONS  TO  DETERMINE  THE  OPTIMUM  FEEC-SPEFO 

MA 

IN 

100 

c 

CO^RINATIONS.  TWO  PERFORMANCE  INCICIES,  COST/PIECE  ANO 

MA 

IN 

no 

c 

PRODUCTION  RATE  ARE  CONSICFRFO  AS  RESPONSFS.  FIRST  THE 

MA 

IN 

120 

c 

PROGRAM  PERFORMS  an  cvnp  ANALYSIS  AND  RECO^ENOS  Two  SETS 

MA 

IN 

130 

c 

OF  OPERATING  CONDITIONS.  NEXT  THE  REGRESSION  COEFFICIENTS 

MA 

IN 

140 

c 

FOR  THE  FOLLOWING  EOLATION  ( OR  A REDUCED  FORM  OF  IT  ) ARE 

MA 

IN 

150 

c 

ESTIMATED 

MA 

IN 

160 

c 

MA 

IN 

ITO 

c 

P.I.  ■ BO  ♦ Bl*LN( SPEED)  ♦ B2*LN( FEED)  ♦ 63 *LN ( SPEED ) 

MA 

IN 

180 

c 

♦ B 4 • L N ( F E E U ) * • 2 ♦ B5*LN( SPEED* FEED) 

MA 

IN 

190 

c 

MA 

IN 

200 

c 

USING  TH^SE  RESULTS,  THE  PERFORMANCE  INOICIES  ARE  PREDICTED 

MA 

IN 

210 

c 

FOR  f EFO  - SPEED  COMBINATIONS  ADJACENT  TO  THOSE  CBSERVED. 

MA 

IN 

220 

c 

THE  PROGRAM  IS  APPLICABLE  TO  SINGLE-0 PE RATION  ANC 

MA 

IN 

230 

c 

*ULT|PLE-OPERATlGN(NCM£RICALLY  CONTROLLED)  MACHINING 

MA 

IN 

240 

c 

PROCESSES. 

MA 

IN 

250 

c 

MA 

IN 

260 

c 

MA 

IN 

270 

c 

REFERENCES  ARE 

MA 

IN 

260 

c 

MA 

IN 

290 

c 

BOX,  G.  t DRAPER,  N.,  EVOLUTIONARY  OPERATION  A STATISTICAL 

MA 

IN 

300 

c 

METHOD  FOR  INCREASING  PRODUCTIVITY  f JCHN  hILEY  G SONS, 

MA 

IN 

310 

c 

1969. 

MA 

IN 

320 

c 

MA 

IN 

330 

c 

BURCHFIELD, P.B.,  MULTIPLE  LINEAR  REGRESSION  , •JOURNAL  OF 

MA 

IN 

340 

c 

QUALITY  TECHNOLOGY' , VOL.  3,NC.  4, OCTOBER  1971. 

MA 

IN 

350 

c 

MA 

IN 

360 

c 

DRAPER,  N.  C SMITH,  H. , A®PLIEO  REGRESSION  ANALYSIS  , 

MA 

IN 

370 

c 

JOHN  WILEY  t SONS,  NEW  YORK,  1968. 

MA 

IN 

380 

c 

MA 

IN 

390 

c 

MA 

IN 

400 

0001 

C0MM(,n  /Ci0i/SPFE0UCOI,FEFD(lC0),SPND(100),NSr0,NFD,  ISP,  IFO 

MA 

IN 

410 

0002 

COMMON  /C1D2/IPH, I CYC, SCYCLE ,RLO, X ( 24 f 100 ) , Y ( 2 , l CO ) , NOBS , NCB1 

MA 

IN 

420 

* , 10P , BSP • BED, NTLC*  TLC ( 20 ) 

MA 

IN 

430 

0003 

COMMON  /C103/C1FF(2,4|,SLM(2,4) ,AVE 1 2 , 4 ) , SUMS ( 2 ) , A VE S ( 2 , 2 ) 

MA 

IN 

440 

0004 

COMMON  /C107/!NU(4) ,SM AX, SHIN, F MAX, FM IN 

MA 

IN 

450 

0005 

Cl  MENS  ION  BIS), TITLE (20) ,BC0(2) ,BP ( 2, 5) . I RRI 2 ) 

MA 

IN 

460 

c 

MA 

IN 

470 

c 

INPUT  MMN  CONTROL  CARD 

MA 

IN 

480 

c 

MA 

IN 

490 

0004 

WRITE! 6,5999) 

MA 

IN 

500 

0007 

IER-1 

MA 

IN 

510 

0006 

REAC(5,5004, FND*9Dl l ,ERR*9012)TI TLE 

MA 

IN 

520 

0009 

WR  I TE ( 6, 6 COO ) T I TLE 

MA 

IN 

530 

0010 

«EA0(5,5001.€Nn*900l,ERR*9C02) 

MA 

IN 

540 
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MAIN 
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•IPH,1CYC,  ICP,  I MAX  , NTLC,  FLO,  AVES(  l,  1),AVES(2,  1 ),BSP,BFC  HA  IN  550 


0011 

IP(  IPF.GT.OIGO  TO  6C1 

MAIN  560 

0012 

WHITE (6,oCCi)IPH 

MAIN  570 

0013 

ItR*2 

MAIN  500 

ooh 

GO  10  6C2 

*A(N  590 

0015 

601 

xRi T H 6, 6002 ) IP H 

MAIN  600 

OOU 

602 

IF(  ICYC.GT.OJGU  TO  603 

MAIN  610 

0017 

WRITt(6,6C03) ICYC 

MAIN  620 

0016 

Its  *2 

MAIN  630 

0019 

GO  TO  6CA 

MAIN  6A0 

U02C 

603 

*kITE(6,600A) ICYC 

MAIN  650 

00  21 

6CA 

IF ( luP.GT.u.ANU.ICP.LT.AIGO  TO 

605 

MAIN  bdO 

0022 

wRI T E ( Of  6005 ) IOP 

MAIN  670 

002  3 

IfcR-2 

MAIN  680 

002  A 

* 605 

GU  10  (606,607,6001 t ICP 

MAIN  690 

0025 

606 

hRI TE (6,6C06) IuP 

MAIN  TOO 

0026 

GO  1C  bCJ 

MAIN  710 

0021 

007 

WRI  TM6,6CC7)I0P 

MAIN  720 

0028 

GO  TC  6C9 

MAIN  730 

OC  25 

60b 

WRI TE ( 6, 60Jb) IOP 

MAIN  7 AO 

0U3C 

o09 

l F ( IMAX.NE.  UGO  TC  61C 

MAIN  750 

00  31 

wk I Tb ( 6, 6C09) 

MAIN  7 60 

0032 

GO  TO  fcll 

MAIN  770 

0033 

610 

xRI TE (6,6010) 

MAIN  700 

Ow>3a 

611 

wRI  ft (6,oCll )RLC 

MAIN  790 

0035 

lF(AVES(l,l).LC.O.O)GU  TO  612 

MAIN  800 

0036 

Wkl TE (6,6C12)AVES (1  ,1  1 

MAIN  810 

0037 

GO  TO  613 

MAIN  820 

0038 

612 

WRITE (6,6013) 

MAIM  330 

0025 

613 

IF(  AVCo(2,i).LE.0.0)GC  TO  6612 

MAIN  8 AO 

00  AC 

WRI  ft (6,661^)AVtS(2,l) 

MAIN  850 

OOA1 

GO  TC  6o  L A 

RAIN  360 

00  A 2 

,6612 

taRI U (6,6615) 

MAIN  870 

00  A3 

60  1A 

I F ( IGP.NE.3IG0  TO  620 

MAIN  880 

00  A A 

IF( dSP.Lb.O.O. GP.6FE.lt .0.0 )G0 

TO  615 

MAIN  990 

00A5 

wkl TE(6,6ClA)B$P,bFD 

MAIN  900 

OOA6 

GO  TO  62  0 

MAIN  910 

00  A 7 

615 

Wkl 1 E ( 6,oC15) 

MAIM  920 

00A8 

1ER-2 

MAIN  930 

00A5 

628 

IF(MLC.EC.0)NTLC*1 

MAIN  9 AO 

005C 

IF( NTLC.LT.l. Ok. N TLC. GT .20 1 GC  TO  7A7 

MAIN  950 

0051 

wRI T E ( 6, 602  a INTLC 

MAIN  960 

0012 

c EAD( 5,5j03,ENu«9CG9,EKP«9010) (TLC ( I ) , I * l , NTLC ) 

MAIN  970 

0053 

00  2276  I-l,.*TLC 

MAIN  990 

005A 

I F ( T L C ( I ) . LT  .0  )GD  1C  2277 

MAIN  990 

0055 

2276 

CONTINUE 

MA I Ml 000 

0056 

wkl  (L (6,6029) (TLC ( I I, 1*1, NTLC) 

MAIN1010 

0057 

GU  fC  uU 

MAIN  lu20 

00  5e 

2277 

I ER  *2 

MA I NIG 30 

OC55 

WHI Tk (6, 5026) 

maimoao 

006C 

GO  TO  6 1 A 

MAIN105Q 

0061 

7A7 

IER-2 

MAIM060 

0062 

WRI TE (6,6C25)NTLC 

MAIN1070 

0063 

6 1 A 

GO  TU  (616,61 7) , 1 ER 

MAIN1080 
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0064 

617 

Wkl TCI6,6C16) 

MAIN1090 

0066 

CALL  EXIT 

MAINIIOO 

C 

MAI N 1 1 10 

C 

input  speed  environment 

MAIN1120 

C 

MAI Nl l 30 

0066 

616 

kEAOI  5,50C2,END-9CC3,ERR*99G4)N$PD 

MAIM140 

0067 

IFMSPU.LE.O.OR.NSPD.GT.IOOIGO  TO  618 

MAI  Nil  50 

0066 

IE I6,6C18)NSPD 

MAINU60 

006* 

CO  TO  619 

VA I Nl l 70 

007C 

6 Id 

Wkl  f E (6#  6C19) NSPD 

*A  I Ml  00 

00  71 

GO  TO  617 

MA I Nl 190 

0072 

619 

REACI5,5003,ENU»9003,ERF-90C4)  C SPEED (1 1 ,1-1, NSPD  1 

MAIN1200 

0073 

DO  62C  I • l , NSPD 

MAIN  12  ID 

0074 

IFCSmEEUIII.LE.C.CIGO  TO  621 

“AIN1220 

00  75 

620 

CONTINUE 

MA I Nl 2 70 

0076 

GO  TO  622 

MAIN1240 

0077 

e>21 

WR1TFC6, 60201 

MATM250 

0076 

GO  TO  617 

MA  I N 1 2 1>0 

C 

MAIN1270 

C 

INPUT  FEED  ENVIRONMENT 

MA  I Ml 280 

C 

MAIN1290 

007 * 

622 

ktACI5,50C2,ENO«SC05,ERR-90C6)NFD 

MAIM1300 

006C 

I FI NFC.tE  .U.Uk.NFL.GT .50) GJ  TO  623 

m A I hi  3 1 0 

ooei 

WRJ  Triof6022)NFD 

MA IN1320 

0062 

GO  TO  624 

MAIN1330 

0063 

62  J 

Wkl  7 C to, oC2l I NFD 

MA I N 1 340 

0084 

GO  TU  61/ 

MA 1 N 13*0 

0065 

624 

REACI5,5J03,END-9005,ERR-9006)  I FEEOf I ) , I * 1,NFD) 

MAIN 1360 

0066 

Ufl  625  I-UNFO 

MA I Nl 370 

00  6 7 

1 r 1 FECOUI.LE .O.OIGC  TO  626 

MA INI  3 89 

ooee 

625 

CON  1 1 NLE 

MAI Nl 390 

008* 

GO  TU  62/ 

MAIM400 

UO*t 

626 

*Rl T E I 6,6023 ) 

MAIM410 

00  5 1 

GO  10  61/ 

MAIN1420 

C 

MMN1430 

C 

INPUT  LIMITS  FCR  SPEED  AND  FEED 

MAIM440 

C 

MAIN1450 

0092 

627 

IF  I IMAX.EO. 11  RE AD  15,5003, ENO-90C7, E RR -9003) S MAX, SMI  4, FMAX  ,F MI N 

MAIN1460 

00*3 

I FI 6FAX.bC.u.0ISM#X-SPEEUlNSPU) 

MATN1479 

0094 

I FI $MIN.FO.O.OISMIN«SPEEDII 1 

MAIN 14 80 

00*5 

I FI FMAA  .tO.U.UlFMAX-FFECINFC) 

MAIN1490 

00*6 

I FI  FMlN.£g.O,OlFMIN«FEEOI 11 

MAIN1500 

C 

MAIN1510 

C 

INSURE  INPUT  SPEECS  AND  FEEDS  ARE  WITHIN  DEFINED  LIMITS 

MA I Nl 520 

C 

MA INI  530 

00*7 

DO  / J-l , NSPD 

MA I Nl 540 

00*6 

IF  ISPECDIJI.LT •SPIN IGU  TO  7 

MAI  Ml  550 

0099 

1 F I J.F0.1IGU  TU  8 

MA  t Nl 560 

01CC 

NSPU-NSPD-JM 

MA I M 570 

oici 

DO  15  K-l , NSPD 

MAIN1580 

01 C2 

6 PE  ED (K l-SPEEOI J 1 

MAI Nl 590 

01 C3 

15 

J-JU 

MAIN1600 

0104 

GO  TO  8 

MAIN 1610 

01C5 

7 

CONTINUE 

MAIN1620 

94 
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01C6 

8 

00  9 J«1,N$PD 

MA I N1630 

01 C7 

IF(  SPEEDI  J).GT.$MAX)GC  TO  10 

MMN1640 

oice 

9 

CONTINUE 

MAIN  1650 

01C9 

GO  TC  11 

MA  I N 1660 

01  1C 

10 

NSPO*J 

MA  I N 16  70 

0111 

11 

00  14  I»1,NF0 

MA I M 680 

01  12 

IF  (FEtOt  n.LT.FMlNJGC  TO  14 

MMN1690 

01  13 

1 F ( I.EC.llGO  TO  12 

MA I M 700 

01  14 

NFU  *NFU-I ♦! 

MAIM7I0 

0115 

DO  lo  J»1  ,NFD 

MA  INI  720 

01  16 

FEE  C(  Jl-FUOl  I ) 

MAI Nl 730 

0117 

16 

l -I  ♦ l 

MMN1740 

01  18 

GO  TO  12 

M A IN  1 750 

0115 

14 

CONTINUE 

MAIM  760 

01  SC 

12 

UU  13  I-l.NFD 

MA  I M 1 7 7 0 

0121 

IFlFttCm.GT.FMAXlGO  TO  17 

MA I Nl 780 

0122 

13 

CONTINUE 

MAIN  1790 

0123 

GO  TO  64 

MA IN  1800 

0124 

17 

N FO  * 1 

MAIMOIO 

C 

MAI N 1820 

C 

PRtPARk  FOR  E V 0 P ANALYSIS 

MAIN1830 

c 

MAIN184D 

0121 

64 

CALL  NLCYCL 

MA I N 18  50 

C 

MAIN1860 

c 

IDENTIFY  INPUT  SPEED 

MAIN 1070 

c 

m-\  IN  1880 

0126 

DO  22  I-liNSPO 

MAIN1690 

0127 

I M XI 1 f UJBl ) • EJ  *5  FEED (t ))  GO  TO  33 

MAIN  1900 

0128 

32 

CONTINUE 

MAIN 1910 

0125 

33 

ISP-I 

MAIM920 

C 

MA IN  19 30 

c 

IOENTIFY  INPUT  FEED 

MAIN1940 

c 

MA I N 1 9 50 

013C 

DO  J4  I * 1 * NFD 

MAI Nl 960 

0131 

IF( X(2,N JBlI.EO.FEEOt I))  GO  TO  35 

MA I N 1970 

>0122 

34 

CONTINUE 

MA I N 1980 

1133 

35 

I FD  ■ I 

MAIN1990 

124 

IF ( ISP.  LT- IND< 1 ) ) INUC l )»  I SP 

MAIN2000 

0 35 

ISPPl*ISP^l 

MAIN2010 

0 06 

I F { ISPPi.GT.IND(2 )) IN0(2)-ISPP1 

MAIN2020 

0117 

IF<  IF -C.LT.IMJ (3)) INCC3> -I  FD 

MA  1 N 20  30 

one 

1 FDP1-IF0*1 

MAIN2040 

013* 

IFC  XFCP1.GT.IN0(4I) IND(4)»1F0PI 

MAIN2050 

C 

MAIN2060 

C 

PERFJKM  THE  E V 0 P ANALYSIS 

MAIN2070 

C 

MAIN2U80 

014C 

CALL  EVUP 

MAIN2090 

C 

MA  I N'2 100 

C 

PEKFJRM  THE  REGRESSION  ANALYSIS 

MAIN2110 

C 

MAIN2120 

0141 

2 F I INC(2)~IN0(1)-I)41,42i43 

MAIN2130 

0142 

41 

*RI  T E ( 6 v 4 6 ) 

MAIN2140 

0143 

GU  TO  9033 

MAIN2150 

0144 

42 

I REGS-0 

MAIN2160 
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OlAb 

GO  TC  44 

MAIN21 70 

U1  46 

43 

IPEGS-1 

WMN2100 

0147 

44 

IF  ( INC  (4  1- HO (3  1-1)  4 1,45,  46 

PAIN2190 

0146 

45 

IREGF-C 

MA  I N2200 

0144 

GO  TO  47 

MAIN2210 

01  tc 

46 

IRLGF-2 

MMN2220 

0151 

47 

IREGMU3EGSUPEGF 

MA1N2230 

0152 

IRtGl-IPEG 

I N2240 

0153 

00  595  I GC«it 2 

MAIN2250 

0154 

I RE  G* J REG 1 

MAIN2260 

0155 

CmLL  HEGRESIX, Y,NCBS, IkEG, 00, 0» I GO) 

MAIN2270 

C 

MAIN2280 

C 

PK 1 :4T  PREDICT  ION  EQUATION  AND  RESPONSE  SURFACE  PREDICTION 

MAIN2290 

. c 

MAIN2300 

0156 

GO  TO  (1, 2,3,4, 5) , I REG 

MAIN2310 

0157 

1 

WRITE(6,‘>l)hU,bm,e(2) 

NAIN2320 

015e 

GU  ro  6 

MA  IN23 30 

0155 

2 

wRI TE  (6, 52130,0(1), tit  2), B( 31,0(41 

FAIN2340 

one 

GU  TO  6 

MAIN2350 

0161 

3 

wRI  TE  Cc,  331*0,611)  i B ( 2 » ,0(3  1 ,0141 

MAIN23O0 

0162 

GU  TO  6 

KA IN? 3 70 

0162 

h 

wRI  Tk (6,54)00,6(1  1,6(21,5(3)  ,8(4), £(5) 

MMN2380 

0164 

GO  TO  6 

v4 IN2390 

0165 

5 

WRITE (6,35)60,6(1), 8(2), 8(3J 

HMN2400 

0166 

o 

CONTINUE 

MAIN2410 

0167 

BOO  ( I C'J ) *80 

MA IN2420 

0166 

1RF (IGO)*IRLG 

MUN2430 

0165 

00  585  1*1,5 

MA  IN2440 

01  7L 

5o9 

up(  iGO,n«am 

MAIN2450 

0171 

595 

CONTINUE 

VAIN2460 

01  72 

CALL  AEGPLT  (BOO,  BP,  IND,  IRR  1 

MAIN2470 

C 

MAI  F'2480 

C 

ERROR  ANALYSIS 

MAIN2490 

C 

*AlN2S00 

01  73 

9001 

WRI TE(o, 90171 

MAIF2510 

01  74 

9033 

CALL  EXIT 

MAIK2520 

0175 

9002 

wRI TE(6,9U18) 

MAIN 2530 

0176 

GO  TO  9033 

MA  IN25‘*0 

0177 

9003 

wRI rE(6,9Cl9) 

MAIN2550 

01  76 

GO  TO  5033 

MAIN2560 

0179 

9004 

WRI TC ( 6, 9020) 

' MAIN2570 

oi  ec 

GO  TO  9033 

MAIN2580 

oiei 

90C5 

wRI TE(6,9C2l) 

MAIN2590 

01  82 

GO  TO  5033 

VA  I N2600 

0163 

9UU6 

wRI TE ( 6, 9C22 ) 

MA  I N 26 1 0 

0184 

GO  TO  5033 

MAIN2620 

0165 

9007 

WRI TEC6, 70231 

MMN2630 

0166 

GO  TO  5033 

MAIN2640 

0187 

9008 

WRI TE(6,9C24) 

MAIN2650 

0166 

GO  TO  9033 

M7IN2660 

0165 

9009 

WRI  TE ( 6, 6C2  7) 

*A  IN267C 

015C 

GO  TO  9033 

MAIN2680 

0151 

9010 

WRI TEIo, >v23) 

MAIN2690 

0152 

GO  TO  9033 

MAIN2700 
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0193 

0194 

0195 
01.96 


0197 

0198 
0195 


02  OC 


02  C 1 


0202 


02C3 
02  C4 
02  C 5 
02  C t 
02C7 
U2C8 
J?C9 
02  10 
02  11 
0212 
0213 
02  1 4 
0215 
02  it 
0217 
02  18 
u2  19 
022C 
0221 
0222 

0223 
02  24 

0225 

0226 

0227 

0228 
0229 
02  3C 


9011 

9012 

C 

C 

C 

48 

51 

52 

53 

54 

55 


WRI TE ( 6,  6030) 
GO  TU  9033 
mRI TEl6,6C31) 
GU  fO  9033 

FORMATS 


UNABLE  TU  DISTINGUISH  SPEcC/FEEO  ', 

PI  • ' , F 1 0. 4, 
' ,F10.4, 


FOkMATI'I  *♦♦  EPRGR  ♦ ♦♦ 

♦'CCPO  INAT  IONS*  I 

FORMAT  l//'0' , TIO, 'TFE  PREDICTION  EOUaTIJM  IS'/'O 
••  ♦ ' ,FU.4, '*LN(  SPEED)  ♦ ' tFlO.4,  •♦LMFEEOI'  ) 

FOh  Hu  (//'  o',  T10,  'THE  PREDICTION  EQUATION  IS'/'O  PI 
••  ♦ • tFlJ.4,  •*uM(SP6fcC)  » • ,F1J.4,»*IMCFEE0)  ♦ SHO^, 

*•♦1  JISPEEU)*LN<Frt0l'/'0',T9, • f • , F 1 0. 4,  • • LM  SPEED  1 * *2  • I 
F Uk  *AT  (//'O'.TiO,  • TFE  PR  tD  ICT  I ON  EQUATION  IS'/'O  PI  » '.F10.4, 
*'  ♦ • ,Fl  3.*,,  «»lM  SPEEOI  ♦ • ,F10.4,'*LMFEED)  ♦ 'tFlO.4, 

♦•  *L  l(.SPrc0**uN(FCEDI'/'O'»T5t  • ♦ • , F 1C.  4,  •♦LNC  SPEEDI  ♦ *Z  • I 
FOR  (//'o',  TiOt  'THE  PREDICTION  EQJATIJi)  IS'/'O  PI  * ',F10.4, 
♦ • ♦ » , F 13  • <i,  • • lN  I SPEED!  ♦ • , F 10. 4, ' *LN< FE ED ) ♦ '9FIJ.4, 

*•  ♦LN(SPlEUMi.NCFEEU)  '/'O', T9,  • ♦ • , F 10. 4,  • ♦IN  ( SPE  ED  I ♦ *2  ♦ ', 
•F10.4,  ' *LN ( FL tD ) ♦ *2 ' I 

FUR  HAT (// '0»,T 10, 'ThC  PREDICTION  EQUATION  IS'/'O  PI  • ',F10.4, 


5U0 1 
5u02 
5UJ3 
5 JC4 
5959 
6000 
60C  1 
60u2 
ouOi 
6004 
60J5 

OUwO 

60C7 
6003 
o0^9 
60  10 
6011 
6012 

6013 

6014 

t#u  1 3 
6016 


♦'  ♦ ' , F 1 J • i , •♦LMSPEEC) 
♦ •♦L  J 1 SPEED  I *L»^I  FEED)  • 1 
FORMAT  (51  3,T2  1,  5F  10.4) 
FUR  MAT ( I J ) 

F OF  VAT  ( 8F 10.4 ) 

FuR  VAT  (2  )A4) 


♦ • ,F10.4,'*lMFEED!  ♦ ' , F ID.  4, 


FOFHAT (• 1 
FukN/.TC// 
FDR  H*f(// 
FOR  HAT  I// 


Fl)K*A)  ( 
F JRNAYI 
F JR  HAT  I 
FOR  NAT  ( 
F Oh  NATC 
FUh  PF T ( 
FOR  Va  1 ( 
F JR  UH 
F Jh  y A I ( 
F JR  1A7( 
FUR  VAT ( 
FOR  MAI  ( 


forvak 

FUR  HAT ( 


ERROR  *♦*•) 
ERROR  ♦ ♦♦•) 


*•*  • ) 


♦ ♦•  ERR  CF 
SINGLE'  ) 

MULTIPLE*  I 
NUMERIC  CCNTPOL  ' ) 


P A C H U P • ) 

, 20A4 ) 

O' ,T1C, 'PHASE  IS' ,7339110,'  ♦ ♦♦ 

O', TIG,  'PFASfc  IS  ' ,T35,I  10) 

CYCLE  IS  ',135,110,'  ♦ ♦♦ 

CYCLE  IS'  ,T35,I10) 

TYPE  CF  PROCESS  IS ' , T35, I 10, • 

TYPE  CF  PROCESS  IS « , T33 1 I 10 , • 
lYPE  CF  PROCESS  I S ' , T35 , 1 10, ' 

TYPE  CF  PROCESS  1 3 ' , 7 35 , 1 10 , ' 

,>FEIU-Ft£0  LIMITS  SPECIFIED') 

SPFfcD-FCEU  LIMITS  NJT  SPECIFIED') 

LflOCR-DVtRHEAO  < $/  HI N ) • , T 35 , F 10 . 4 ) 

COST  STO.  DfcV.  tSr.  IS', T39 , F 10.4 ) 

CCST  STD.  DEV.  EST.  IS  NOT  SPECIFIED') 

RASE  SPEED  IS' ,T35,F1C.4/'0',TIO, 

♦ •BaSE  FIED  IS',  T35,F1C.4) 

•iASk  SPEED  AND/OR  FEED  INCORRECTLY  S°ECIFIEO') 
♦•*  PRUGRAH  TERMINATICN  DUE  TO  SPECIFIED  ERROR', 
" CCNDI f IONS  *♦•' ) 


,T  6 1 
O'  ,T2 


TIO, 
T 1 0, 
TIO, 
TIO, 
TlO, 
TIO, 
UO, 
TIO, 
T 1 J , 
TIO, 
UO. 
TIO, 
IS'. 
TIO, 
UO, 


MAIN27IO 
MA I N2720 
MAIN2730 
MAIN2740 
MAIN2750 
MAU2760 
PAINZ/70 
MAIN2700 
PAIN2790 
VM  N2800 
VAIN2810 
MAIN2820 
PAIN2830 
"AIN2840 
PAIN2850 
MAIN2860 
PA  I f. 28  70 
MAIN2680 
PMM2890 
PAIN2900 
MAIN2910 
PAIN2920 
MAIN2930 
MA IN2940 
PMN29  50 
PAIN2960 
MAIN2970 
MAIN2980 
MAIN2990 
MA I N3000 
P.AIN3010 
MAIN3020 
MAIN3030 
MAIN3040 
MAIN 30 50 
MA IN  3060 
MA  0.30  70 
PAD  3080 
MAIN3090 
MAIK3100 
MAIN3110 
MAIN3I20 
MAIN3I30 
MAIN3140 
MAIN3150 
MAIN3160 
MAIN3170 
PA  I K 3 1 8 0 


60  18 

FORMAT ( 

•O' 

,710, 

•NCPBER 

OF 

SPEEDS 

I S ' , T35 

,113) 

MAIN3190 

6019 

FOR  HAT  ( 

•O' 

,T  i 0, 

•NUPDtk 

UF 

SPEEDS 

IS' ,735 

,110, • ♦♦♦  ERRCR  •*♦•  ) 

PAIN3200 

6020 

FOF  HAT ( 

»o* 

,110, 

•SPECIFIED 

3 PEfcOS 

ARE  L/E 

ZERO' ) 

MAIN3210 

6021 

FORMAT C 

1 3* 

, TIO, 

•number 

OF 

F t CDS 

IS' • T 35  , 

110,'  ♦♦•  ERRCF  ♦♦♦•) 

MAIN3220 

6022 

FORMAT! 

• J' 

, T 1 0, 

• NUMBER 

CF 

FEEDS 

IS' , T 35, 

110) 

MAIN3230 

6w23 

FOR  HAT! 

•O' 

, T 10, 

•3PECIFI  ED 

FEEDS 

ARE  L/E 

ZERO'  ) 

MAIN3240 
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0231 

6026 

FOR  *A  M • J* 

0222 

6023 

FUR  HAT  I' O' 

02  33 

6026 

FuR  vat  (•  J' 

02  36 

6027 

FOP  M A 1 t ' 1 
♦,  /'  O'  ,Tl >, 

0225 

6028 

FORMAT!  ' l 
* * T l 9 f 1 Cl 

02  3< 

6029 

FUE  MAT ( « J ' 

0237 

oO  JO 

FOE  *aT  I • i 
*,  /'  O'  * T i ), 

0238 

6031 

FUR  HAT ( • 1 
• , T 1 ) f ' u J 

0225 

6ol3 

FOR  vat  (•  J' 

026C 

6615 

F )P  VAT(  ' J' 

0261 

9 011 

furv  'A1  I ' l 
*,/'  J*  ,U  ) f 

0262 

9018 

FUf  MAT  I ' 1 

* » T i )»'  u 

0263 

9019 

F0>»  VAT  ( ' 1 

*,/' O' . f 19, 

0266 

9u2U 

F Q*  HAT  ('  l 

• ,ti  't  4 i 

0265 

90  21 

FOE  vat  (•  l 
»t/'C',TUt 

0266 

9u22 

RUE  -AT  1 ' 1 
fl  if'CN  I 

0267 

9023 

FORMAT  I ' 1 
*,/'  J'  , 119, 

0268 

9L26 

FUK  HAT  ( ' 1 
* , T 1 9 , 'UN  I 

0265 

END 

EHisCft 
CC51S 
ERRLR 
CLOT  If. 
'SPECIE 
fcRUCK 
b CCRC 
ERROR 


*♦  * 


)LS  IS'  » T35f  I 10»  • ***  ERROR  ***•) 

1L  COST  NEGATIVE  ERROR  ♦ **') 

unexpected  termination  of  cat  a 

5'  ) 

INVALID  CHARACTER  ENCOUNTERED  '/'O' 


UNEXPECTED  TERMINATION  CF  CATA 

)'  ) 

INVALID  CHARACTER  ENCOUNTERED  '/'O' 


'P.R.  STO.  DEV.  ESf • I3',T35,F10.6) 

'P.^.  SfU.  DEV.  b ST . I*  NUT  SPECIFIED') 
t RP  L’K  UNEXPECTED  TERMINATION  OF  DATA 


ERROR 
PFCDLEM 
h RR  LK 
r SPCEC 
fPKCR 
SPEED  P 
E RRO 
T FELD 
f RRLR 
F E EC  P* 
E R R C W 
r SPfcEi 
LRRCH 


INVALID  CHARACTER  ENCOUNTERED  '/'O' 

) 

UNEXPECTED  TER  v I NATION  OF  CAT* 

:TE  K S MISSING' I 

INVALID  CHAPACTEP  ENCOUNTERED  '/'O' 
:R  S'  I 

UNEXPECTED  TERMINATION  OT  C*TA 
rERS  HISSING') 

IIVALlO  CHARACTER  ENCOUNTERED  '/'O' 
:S»  ) 

UNEXPECTED  TER  FI  NAT  I ON  OF  CATA 
.PUTS  MISSI  IG'  ) 

INVALID  CHARACTER  ENCOUNTEKEC  '/'O' 


MAIN3250 
MAIN3260 
MA I N 32  7u 
MAIN328D 
MA I N 3790 
HA  I M3  ADO 
MAIN3310 
MAIN3 j20 
MAU3330 
PI  I N3360 
M \ I N J 3 5 0 
VUN3360 
H A I N 3 3 7 0 
MAIN338U 
MAIN33D0 
MAIN  3600 
•MA  IN3610 
MAIN  3620 
MAIN3A30 
VMNJ660 
,h\INJ650 
*4lN3660 
MA  IN  H67U 
M A IN  36 80 
VA INJ690 
MAINJ5D0 
MAIN3510 
v A I N35?0 
M A IN  3530 
*AIN3560 
MAIM3550 
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NUCYCL 


0*TE  • 74053 


10/32/02 


uOCl 

SUBROUTINE  NUCYCL 

NUCY 

10 

C 

NHCY 

20 

C 

NUCY 

30 

c 

SUBRUUTINE  NUCYCL  15  USEO  FOR  PROCESSING  ALL  PHASE/CYCLE 

nucy 

40 

c 

COMdlKAlIOUS.  NUCYCl  PREPARES  THE  INPUT  DATA  FOR  BOTH 

NUCY 

50 

c 

T Ht  EVJP  ANU  REGRESSION  PROCESSES. 

NUCY 

60 

c 

NUCY 

70 

c 

NJCY 

80 

00C2 

DIMENSION  TEMPI 4 « 24 ) 

NUr  Y 

90 

00C3 

C 01  MON  /Cl  02/  If’N,  ICYCtSCYCLEfRLCtM  24*  loO)«  YI2,  100)  , NUfi  5,  NDBl 

NUC  Y 

loo 

»t  IOPtRSPvaFi)tNll(  « T 1 C (20) 

NUCY 

110 

00C4 

COMMCN  /C103/0IFF  (2f4)fSUM(  2, 4 ) f AVE  i 2f  4)  # SUMS 1 2 ) t AVES  (2  f 2 ) 

NUCY 

120 

00C5 

CU1VUN  /C  1C7/III0I4)  fSMAXfSMlNfFMAXtFMlN 

NUC  Y 

no 

OOC6 

702 

SCYCt  E*  IC  YC 

NUCY 

140 

00  C 7 

NTl  CF«N7teC^4 

NUCY 

150 

ooce 

DO  9 1*1,2 

NUCY 

160 

0Ct9 

AVES (1,2) *0.0 

NJCY 

170 

OOIC 

9 

SUMS! I )»0.0 

NUC  Y 

100 

C 

NUCY 

190 

C 

INPUT  Nt A OBSERVATIONS 

NUCY 

200 

C 

NUCY 

210 

00  11 

NUBS-1 

NUCY 

220 

0012 

703 

NUB 1-6CB5  *1 

NUCY 

230 

00  12 

NOB  4*f.0BS  ♦ 3 

NUCY 

240 

0014, 

F EAC(5,53C2»EN0«9C1 7,ERP-90 16)  1 XU  ,NOB$)  , J- 1,NTLCP) 

NUC  Y 

250 

0015 

l F 1 X( IfNORS) .EG .0.12345)00  TO  701 

NUcY 

260 

0016 

UU  66  J *NCR 1 , NUB4 

N"CY 

270 

0017 

66 

PEA  Cl  5, 5302 9 END-901 5, ERF-90  16)  ( X (I , J ) , 1 ■ l,  NTl  CP ) 

NUCY 

230 

OG  IB 

NOtt  S* NLbS  ^4 

NUCY 

290 

0019 

CO  TO  703 

NUCY 

3 JO 

C 

NUC  Y 

J10 

C 

INPUT  THE  U«1A  FROM  THE  PREVIOUS  CYCLE 

NUCY 

320 

C 

NUCY 

330 

0020 

701 

K EAlS  1 5t!»3Ci«  END-901  lvERR»90 12) 

NUCY 

340 

MISJMUtl  )f  1-1,4) , SUMSU)  9 < AVEIJ,! ) 9l*lf4)9AVESU92)9  >1»2)  v 

NUCY 

3>0 

-UNCI  I),  1*1,4) 

NUC  Y 

3 SO 

00  21 

I F | AVES(1«2)*NE  .U.OAVESI  1,  1)-AVES(  1,2) 

NUC  Y 

370 

0022 

lFUVbS(2  92).NE.0.C)AVES(2,  D- AVES  1 2,2) 

NUCY 

3B0 

C 

NUCY 

390 

C 

SCRT  PU1NT5  INTO  CURT  EL  T SkOUkuCt 

NUCY 

4 jO 

C 

N»  ICY 

410 

0023 

9017 

NUh  i-NJ'IS-l 

NJCY 

420 

0024 

NJBi-NLJS-3 

NUCY 

430 

0025 

UU  45  ! •! ,4 

NUCY 

440 

00  26 

00  45  J«1,NTLCP 

NUCY 

45o 

0027 

45 

TEMPI  1,J)«XU, NOH l ♦ 1 — 1 ) 

NUCY 

460 

00  2 € 

SPMI-TEM^Cl,l ) 

NUCY 

470 

0029 

SPM4«TE«-M  l,  l ) 

NUCY 

4 dO 

003C 

F DM  I -Tfc  MP  1 1,2) 

NUC  Y 

490 

OQjI 

F DM  A* TEMPI  1,2  ) 

NUCY 

5 JO 

0032 

DO  10  I * c,4 

NUC  Y 

410 

0033 

IMTIMPU  ,1 ) .LE.SFKT  )SPM1 -TEMPI!  ,1) 

NUC  Y 

520 

0034 

IF!  f L MF ( I ,l).Gt.SPM*) SPHa -TEMPCI,!) 

NUCY 

5?0 

0o35 

lF(Tt*PU  ,2).Lt.FCMI  IF0M1 -TEMP (1,2) 

NJCY 

540 
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0036 

I F(TEKP(I ,2) .GE.FCMA)FOMA*TEMP( I f 2) 

NUCY  550 

0037 

10 

CONTINUE 

NUCY  560 

0036 

00  11  I * 1 f A 

N’JCY  570 

0025 

J*0 

NUCY  580 

UUSO 

IF  (1  EMP(  I,  l)  .6  0.SPM  .AND.  TEMP  (I  # 2)  .EO.FOMI ) 

J-l 

NUCY  590 

00  A 1 

IF  (IE  AIM  I til  • fcQ.SP*'A.AN0.TEMP(I,2)  .EU.FUMa) 

J *2 

NUCY  600 

00A2 

IF  (TEMP  (I  ,1)  .f.0.SP6A.ANU.TEMP(l,2)  •CO.FOMI  ) 

J*3 

NUCY  610 

00  A3 

IF  (TEMPI  1,1)  • F0.SPM.AHU.TEMP(I,2I .EO.FOMA) 

3*A 

NUCY  620 

00  A A 

I F(  J*  EC. J ) Gu  TO  15 

NUCY  630 

C 

NUCY  6A0 

C 

CALCULATE  THE  RESPONSE  VARIABLES 

NUCY  650 

C 

NUCY  660 

00A5 

DO  il  K*i, NTLCP 

N.JCY  670 

00A6 

X(  K ,NCBU  J-U*TEKF(  I,K) 

NUCY  68C 

00A7 

11 

CONTINUE 

NUCY  690 

uo  Ad 

00  A6  1*1, NOBS 

NUCY  700 

00  AS 

Y(2,I)*A(3,II/X(A,I  ) 

NUCY  710 

005C 

Y ( L , I )*K.O*X(  A,  I ) 

NUCY  720 

0051 

00  A/  3*5 , NTLCP 

NUCY  730 

0052 

A 7 

Y( 1 , A )*Y(  1 , I) ♦Xf J,I  ) ♦TLC ( J-A ) 

NUCY  7 AO 

0053 

A6 

Y(  1,1  ) » Y ( 1,  1)/X(3,I ) 

NUCY  750 

C 

NUCY  760 

C 

PLOT  1HC  INPUT  FOINTS 

NJCY  77C 

c 

NUCY  780 

00  5A 

CALL  0PG°T3(SPMI ,SP6A,FDMI,  FCMA,  3) 

NUCY  790 

c 

NJCY  800 

c 

PRIi^T  THE  INPUT  DATA  ANC  RESPONSES 

NUCY  810 

c 

NUCY  820 

0055 

WRITE (6,6C02) 

NUCY  830 

0056 

1)0  60CA  J * l , NOBS 

NUCY  8A0 

0057 

w hi  Tfc ( o, jC  03) ( X ( I , J), 1*1, NTLCP I 

NUCY  850 

0056 

600A 

WKI Tf (6, j005) ( Y(K,J ),K*l,2) 

NUCY  860 

t 

NUCY  to 70 

C 

TA&JLAfE  THE  FCSPCNSES 

NUCY  8to0 

c 

NUCY  890 

005S 

CALL  FLCTPMINJBl) 

NUCY  900 

0060 

I F ( ICYC.EC.DGO  TO  12 

NUCY  910 

006  1 

RETURN 

NUCY  920 

c 

NUCY  930 

c 

PREPARE  THE  OATA  FOR  E V 0 P 

NUCY  9 AO 

c 

NUCY  950 

0062 

12 

00  IA  1*1, A 

NUCY  960 

0063 

UO  IA  J*1 , 2 

NUCY  970 

006  A 

SUMCJ, I)*0.0 

NUCY  980 

0065 

IA 

AVfc(J,  11*0.0 

NUCY  990 

0066 

RETURN 

NUCY 1000 

0067 

9011 

rlRL  TE  (6,9027 ) 

NUCY1010 

006  £ 

GU  TO  9033 

NUCY 1020 

0069 

9012 

WKI TE( 6, 9C28) 

NUCY 1030 

007C 

GO  TO  9033 

NUC  Y 10 AO 

0071 

9C13 

WRI  (6(6,  )C29) 

NUCY 1050 

0072 

GO  TO  5033 

NUCY 1060 

0072 

90  l A 

WRI TE(6, 90301 

NUCY1070 

007A 

GO  TO  9C3  3 

NUCY 1080 
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0075 

0076 

0077 
00  76 
GO  7S 
ooet 
OOtil 
0062 
0063 
0065 

0065 

0066 
0087 
0086 


0089 

OUSC 

0051 


0052 

0053 
0054. 

0055 


00.56 

0057 

0056 

0055 

01CC 


0LC1 


01C2 


0 1C  3 


9015 

9016 
9033 
15 


17 

18 

19 

C 

C 

C 

53ul 

5 3v2 

6 U02 


WRI TE ( 6,9031 1 
GO  TO  5033 
wRI  Hit,  >032 ) 

Call  exit 

JK«l 

I FI  TEMPI  I ,1)  .NL.SFM1.ANU.TEMF1I,  1 ) .NC.SPMA ) JK»3 
IM  TEMP(  4 ,2)  .Nt.FCMI  .AND.  TEMPU  ,2)  .NE.FDMA)  JK»JK*l 
GO  TL  <1 7,17,  18,19), JK 
WRI  T f ( C , 9 C 35  ) 

GO  TU  9033 
WRI Tt (6, 9C35) 

GU  TU  903  3 
WRI  TUG,  7036) 

GO  TO  5033 


HJRMA1S 


6003 

oOC5 

9027 

9U2H 

9029 

9030 

9031 

9032 
9035 

9035 

9036 


FORM*  T(G6  13.B/6tl3.fl/tC 13. a 
EOR*A(  I OF  10.5  ) 

FUkMaI  (*1*,T20,31X,  • INPUI 
♦*0*,*  SPtfcO  FEED 

♦ *T  J 3L-E0G ES * , T 1 35 , * C J 

* • I 

♦ • 19)  110)' ,1105,  • (S/P 

* /•  * , T 5 l,  * III) 

120)' ) 

( * 7* ,FIO.O,F IC.9,5X,2 
(•  ♦ *,1  405, F 10. 2, F 10.5 
(•l  ♦ ♦♦  tRRCR  ♦ ♦♦ 

T 19, ‘PUNCHED  LUTPUT  F 
( • 1 ♦ ♦♦  tRRCR  ♦♦♦ 

IM  PUNCH  UUTPLT  FPGM 
(•1  ♦ ♦♦  ERRLk  ♦♦♦ 

T19, ‘PUNCHED  LUTPUT  F 
<‘l  ♦ ♦•  EFRCK  ♦ ♦♦ 

IN  PUNCH  CUT  PUT  FKJ* 

I 1 L ♦ ♦♦  ERPCR  ♦ ♦♦ 
r 17, »NEw  CbStkVATIONS 
|‘l  ♦ ♦♦  FFFCR  ♦ ♦♦ 

IN  NEw  ObSLRVATIuNS  - 
l‘J  ♦ *♦  FRRCR  **♦ 
IN  THE  NEW  C6SERVATI 


/2E  13.8, 512) 

VALUES  AND  COMPUTED  RESPONSES  •/// 

parts  time* ,t*i,2ox, 

ST  PRTO  RATE*/*  *,T5l, 

1)  (21  (J)  (5)  (*)  (6)  (7)  CHI* 

IECE)  (PI  fcCES/f«IN)  • 

(12)  (13)  (19)  (15)  (16)  (17)  ( 18 ) • , 


♦*  (19) 

FOR  4AT 
FOR  M 4 T 
FU*'  MAT 
♦,  /•  J 
FORMAT 
* , T l 9,  • 
FOE 

*»/•  O'  , 
FORMaT 
♦tTl  7,  ‘ 
FORMAT 
*,/'  J*  , 

FOr.MAl 

*•  T 1 5,  ‘ 

FOR  **T 

♦ ‘FuUNO 

•'  ) 

F OR  M AT 

♦•FOUND 

♦ , • T ER  M 
FORMAT 

♦•more 

♦ * • 

ENO 


F 10. 0 , ( * • ,T5l,10F5.0) ) 

) 

UNEXPECTED  TERMINATION  OF  DATA 
ROM  PREVIutlS  EV3P  CYCLE  MISSING1) 

INVALID  CHARACTER  ENCOUNTERED  ‘/‘0‘ 
PREVICJS  E V J°  CYCLE* ) 

UNEXPECTED  TERMINATION  Of  DATA  ..• 
RCM  PHE-VIOUS  REGRESSION  STEP  MISSING*) 
INVALID  CHARACTER  ENCOUNTERED  ‘/‘0‘ 
PREVIOUS  R EuRESS I CN  STEP*) 

UNEXPECTED  TERMINATION  OF  CATA  ..* 
MISSING* ) 

invalid  character  lncountfpec  */*o* 

- CHECK  YOUR  KEYPUNCHING* ) 

-MUR  E THAN  Two  FEED  SETTINGS  SERF  *, 
QNS  * / • 0* ,T 10 , ' PPGGKAM  TERMINATING 


( • U ♦♦♦  ERRCR  ♦♦♦  MORE 

1 X TML  NEK  CBSERVATIONS 

INU  ING*  ) 

(')  ♦♦♦  ERRCR  *•♦  MCRE  THAN  TWO 

Th*N  TWU  Ft  EC  SETTINGS  WERE  FOUND  IN 
/ * 0 * ,T1 0, ' FRCGPAM  TERMINATING*) 


than  Two  SPEED  SETTINGS  WERE 

* / ' 0* , T 10 , * PROGRAM  • 


SPEED  SFTT  IKGS  AND  *, 
THE  NEW  03SERVATIUNS* , 


NUC Y 1090 
NUCY1100 
NUCYll 10 
NUCYU20 
NUCYll  30 
NUC Y l 150 
NUCYll 50 
NUCYll  60 
NUCY1170 
NUCYll 80 
NUC  Y 1 190 
NUC Y1 200 
NUC Y 12 l 0 
NJCY1220 
NUC Y 12 30 
NUC  Y 1250 
NUC V 1 2 50 
NUC Y 1260 
NUC Y 1270 
N’JCY  1280 
NUCY1290 
NUCY 1300 
i NUC Y 131 0 
NUCY (320 
NUCY 1330 
NUC Y 1350 
NUCY l 3 50 
NUC  Y l 360 
NUCY  1 370 
NUCY 1380 
NUC  Yl 390 
NUC Y 1500 
NUC  Y 1 5 1 0 
NUCY  1*20 
NUCY 1530 
NUC Y 1550 
NUCY 1550 
NUC  Y I960 
NUCY  15  70 
NUCY1580 
MIC  V 1590 
NUCY 1500 
NUCY 15 1 0 
NUCY l 520 
NUCY 1530 
NUC Yl 550 
NUCY 1550 
NUCY 1 560 
NUC  Y 1 5 70 
NUCY 1580 
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OOC1 

SUBROUTINE  EVOP 

EVOP 

10 

C 

EVOP 

20 

C 

EVOP 

30 

C 

SUBROUTINE  EVCP  PEP  FORMS  THE  EVOLUTIONARY  OPEPATICN  ANALYSIS 

E V IP 

40 

C 

fOR  THE  LrtPUl  OBSERVATIONS.  THE  PROCEDURE  USER  ANC  THE  RESULTSEVOP 

50 

C 

PEPU/UtU  ARE  PC0IFICA1ICNS  PE  THPSt  GIVEN  BY  BOX  ANO  OPAPtK. 

EV7P 

60 

c 

FVOP 

JO 

c 

EVOP 

80 

uuo  2 

COiVCN  /ClOl/SFEECllOC)  ,FEEC(  100  ),$PND(  100) , NSP0,r FO,  ISP,  IFO 

EVOP 

90 

OOC3 

COR*LN  /C  1 02/ I PH,  1 CYC,  SC  YLLF,RL0,X(2<.,  LOO),  Y(2, 100)  , NCPS,  NOBl 

EVOP 

100 

*,  IP  P , PGP , liFli,  LC  ,TLC  ( 20) 

EVOP 

110 

OOC  A 

CG^CN  /C  103/01  EF  (2  ,4  ) , SUM  ( 2 , 4 ) , AVE  ( 2 , 4 ) , S'J  4S  ( 2 ) , AVES  (2,2) 

EVOP 

120 

OOL5 

LG*  <UN  /C104/F4IM  10) 

EVOP 

130 

00  tc 

CO  1 MC  N /C  107/  I Nil (4  ) ,5 MAX,  SM  IN,  FMAX,  FMIN 

EVOP 

140 

OOC  7 

DIMERS  iC  4 KA.Gt  (2  ) , SLAVES  C 2 I , SNA  VES  C 2 ) , SNSUM  ( 2 ) , SNE  W ( 2 ) , S EMFANl  2 ) 

EVnp 

150 

c 

CV'P 

1 60 

c 

BEGIN  POINTING  EVOP  TABLE 

CVOP 

170 

c 

EV3P 

180 

OOC8 

i F.s  r- i 

CVUP 

190 

oocs 

IF(  ICYC. ;E.2IIES)*2 

FVOP 

200 

OOIC 

WRI TE (6,6600) 1PH, ICYC 

EVOP 

210 

00 11 

N2*  NUB i ♦ L 

EV-IP 

220 

u01< 

N3  - »\2  ♦ 1 

CVOP 

2 30 

0013 

rt(6,o553)  (X(  l,NCBl),X(l,N2)  ,X(l,N3)  ,X(1,NCBS)  ,1  *1  ,2  ), 

FVOP 

240 

*(X(2,MC0l)  ,x(2,N2),X(2,N3),X(2,N0r)S),I«l,2) 

EVOP 

250 

001** 

JLYC*SCYCLE 

FVMP 

260 

v/015 

I F ( JCYC.LT.2) JLYC»2 

EVOP 

2 70 

OJU 

IM  JCVC.if  .1J)  JCYC-10 

E V TP 

280 

0017 

wkl  U (6,6601)  ( ( SUM  J,  I),  I-l  ,4)  ,J*l,2) 

EV  JP 

290 

oo  ie 

WRI TE (6,6602)  ( ( AV E ( J , I),I-l,4),J-l,2) 

EVOP 

300 

0C1S 

SCYCL  e* ic  yc 

FVOP 

310 

c 

EVOP 

320 

c 

CCPPJfE  THE  DIFFERENCES 

EVOP 

330 

c 

tv  up 

340 

002C 

DO  10  1-1,4 

FVOP 

350 

0021 

UIFF( 1 ,1)«AVE (1,1  )-Y( I , NOB  1 ♦ I- 1 ) 

EVUP 

360 

0022 

10 

D IFF (2,1) -AVE (2, 1 )- Y ( 2, NOB  1 ♦ I-l ) 

EVOP 

J70 

c 

F VHP 

380 

c 

CALCJLATE  THE  RANGE 

EVOP 

390 

c 

FVOP 

400 

0023 

DO  12C  1*1,2 

EV  JP 

410 

0024 

YMA  x*AMAXl(l>IFF(I,l),CIFF(  I , 2 ) , 0 IFF  ( I , 3 ) , DI F F ( I , 4 ) ) 

EVOP 

420 

0025 

YMIN«AM  YllUlFF  (1  ,1  ) , C I T F ( I , 2 ) , U I FF  ( I ,3),0IFF(  I,  4)) 

EVOP 

430 

0026 

120 

RANGE ( l ) -YMAX-YMIN 

EVOP 

440 

0027 

I F ( IC  YC  .c  T *2)GJ  TC  44 

EVOP 

450 

c 

EVOP 

460 

c 

CALCJLATE  NEW  S,  NEW  SUP  S,  NEW  AVFRAGE  S 

EV  IP 

4 70 

c 

EVHP 

480 

0028 

00  466  T-1,2 

Evnp 

490 

0029 

SHE Wl  I ) * <ANGL ( I ) *F4N( JCYC) 

EVOP 

500 

003C 

$ NS  UP ( I ) - ( SNfc* ( I ) ♦ SUMS ( I ) ) 

EVUP 

510 

0031 

1 F 1 SUMSI I ) .EO.C.O)SNSLM(I )-SNEW(I )*2. 

EVOP 

620 

0032 

SHAVES ( 1 ) - (SNSUP(I) )/ 2 • 

evnp 

530 

0023 

466 

CONTINUE 

E VHP 

640 

102 


FORTRAN  IV  G LEVEL  21 


EVOP 


OATE  - 7 3j6  1 


15/21/16 


0036 

GO  TO  65 

EVOP  550 

0035 

66 

DU  667  1-1,2 

EVOP  560 

0036 

SNEM1  )«J.O 

EVOP  570 

0037 

5NSUM!I)*SUrtSm 

EVOP  580 

0038 

667 

SNAVES(l)-AV(S(I,l) 

EVOP  590 

C 

EVUP  600 

C 

CALCULATE  2 STANDARD  ERROR  ESTIMATES 

EVOP  610 

c 

EVOP  620 

002S 

65 

00  668  1-1,2 

EVUP  630 

U06C 

5E*  VES ( I ) -2 • /SUKl (SCYCLfc)*SNAV£S(I  ) 

EVUP  660 

0061 

SfcMEAM  I)-l.7  3/SORT(SCYCLE)*SNAVESC  I ) 

EVOP  650 

0062 

IF  1 SNA VES ( I 1 .EO.AVESl I, l) )S NAVES! I ) *A VES! I ,2) 

EVJP  660 

0063 

668 

CONTINUE 

EVOP  6 7u 

C 

EVOP  680 

C 

COMMUTE  SUMS  AN C AVERAGES 

EV  IP  690 

C 

EVUP  700 

0066 

67 

00  11  1-1,6 

EVOP  710 

0065 

UO  11  J -1 ,2 

F.ViJP  72  C 

0066 

SUMJ  , 1)  «$Uil  J,  I ) ♦ Y ( J , NOR  1 I — 1 ) 

EVOP  730 

0067 

11 

AVE(J,i)« SUM ! J , I ) /SCYCLE 

F V IP  760 

C 

EVUP  750 

C 

FINISH  PRINTING  THE  EVOP  TABLE 

EVOP  760 

C 

EVUP  770 

0068 

WRI  TE ! o,6603) ( ( Y! J, I ) , I -NOft 1 , N TBS ) , J *1 , 2 ) 

EVOP  700 

0065 

WRITE  !o, 6606 ) ( ID1FF (J, 1 ), I- 1,6) , J«l, 2) 

EVUP  790 

005C 

WRI  rb(6f66C5)  ! ! SIP!  J,  I ),I  -l  ,6)  , J-l,2) 

EVGP  800 

00  51 

aRI TL 16, o6Uul ! IAVE! J, I ) , I«l,6) , J*l,2) 

EV  )P  810 

0052 

WRITE !o ,6620) 

EV  TP  820 

C 

aRI TE!6, 6621) SUMS ( l ), SUMS !2) 

EVnP  830 

0053 

*R) Tt !o,o622 )AVEU 1, IEST ) , A VES! 2, 1EST) 

EVOP  06 0 

0056 

WRI  rt!o,6623)SNEW(i),SNFWl2) 

EVOP  850 

0055 

WRITE  l 6, 6626 ) RaNG  E ! 1),RaNGE(2) 

EVOP  060 

0056 

WRITE  (o,6  625  )^NSUM(  I)  fSNSU'4C2) 

EVOP  87C 

0057 

WRI  fC  (6,6626)  SHAVES (l  ),  SNA  VES  (2) 

FVOP  880 

OU56 

WRITE (6,6607)  LEAVES 

EVOP  890 

C 

EVOP  900 

c 

SAVE  COMPUTED  VALLES  FOR  FUTURE  CYCLES 

EVOP  910 

c 

EVOP  920 

0056 

WRITE(7,7COO) 

EVOP  930 

*(  ISLMt  J,  I ) , 1-1,6)  , SNSUM  ( J ) , (AVE(J,1  ), 1-1,6),  SNAVESU  ) ,J«l  ,2  ) , 

FV OP  960 

*(INCI  I), 1-1,6) 

EVOP  950 

c 

EVOP  960 

c 

ANALYZE  THE  EFFECTS 

FVOP  970 

c 

EVOP  900 

006C 

DO  166  1-1,6 

EVOP  990 

0061 

166 

AVE(2,I)—  AVfc(2,l) 

FVOP 1000 

0062 

I GO  - 1 

EV0P1010 

0063 

ISPi-ISP 

EV0P1020 

0066 

IFUi-IFD 

E VO PI  030 

0065 

165 

GO  TO  ( 166, 167) , IGO 

fc V TP  1060 

0066 

166 

wRITE(6, 168) 

EV0P1Q50 

0067 

GO  TO  169 

EVOP 1060 

0068 

167 

WRI TE (6, 150) 

EVOP 1070 

0065 

1*9 

I F(  SEAVESUGO)  .EG.O.O)GO  TO  999 

EVOP 1 080 
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c 

EV0P109U 

c 

calcjlatc  the  effects  of  the  variables 

Evnptioo 

c 

EV3P1110 

007C 

EFFA*tAVEI  IG0f2M/VH  IGO,  3 I-AVEI  IGO,  l )-AVE(  IGC,  4)  1/2. 

FV0P1120 

00  71 

fcFF  B*UVfcl  IGu  ,2  J»AVE(  IGC  , 4 )-A  Vl  I IG  J,n-AVE(  I GC  i 3 1 )/2  • 

EV0P1130 

0072 

EFF  Ab«(AVE( IGU, 1 J ♦/VfcUGU,  2 )-AVE  (IGO,  2I-AVE  i IGO,  4 ) ) /2  . 

EVOP1I40 

OC73 

I A*  0 

EVOP 11 50 

0074 

I F| AbS I EFFAI-SEmVES I IGO 1)30,21,31 

CV0P1160 

c 

EV0P1170 

c 

THE  EFFECT  OF  SPEEO  IS  NOT  SIGNIFICANT 

tviPiiao 

c 

EV3P1190 

0075 

30 

WFI  TE I o ,66l G ) EF  FA 

EV  OP  1200 

0076 

GO  TO  32 

E V JP12 10 

c 

EVOP 1220 

c 

THE  EFFECT  CF  SPEED  IS  SIGNIFICANT 

EVOP 1230 

c 

EV0P124U 

0077 

il 

I A«  I A ♦ 1 

EVTP1250 

u076 

WRI  TE I c, 66 1 1 ) EFFA 

FV0P1260 

0075 

32 

1 FI  AES (EF FBI- SLAVES  1 1 GO)) 33 » 34 • 34 

EVTP1270 

C 

EV^P  1280 

c 

THE  EFFECT  UF  FEEO  IS  NOT  SIGNIFICANT 

EVQP1290 

c 

EV  JP1300 

0C6Q 

33 

*R1  TE(6,6oU)EFFB 

EV0P1310 

0061 

GO  fO  35 

EVJP1320 

C 

EVOP  13  30 

C 

THE  EFFECT  OF  FEED  IS  SIGNIFICANT 

EVOP 1 340 

C 

EV0P135O 

0062 

34 

I A*  IA*2 

EVOP 1360 

0063 

*R1 TI  I 6»66  L3) EFFB 

EVJP1370 

00  6 A 

35 

IF!  ABS(EFFAi5)-SEAVE$UG0)  13  6,37,37 

EV0P1380 

C 

EV0P1 390 

C 

TFE  EFFECT  OF  THE  INTERACTION  IS  NOT  SIGNIFICANT 

EVOP l 400 

C 

EVOP 1410 

0065 

36 

WRI  TE  I 6, 66  14 ) EFFAB 

EV JP1420 

one 

GO  TO  38 

EVOP I 43 C 

C 

EVOP 1440 

C 

TFE  EFFECT  OF  TFE  INTERACTION  IS  SIGNIFICANT 

EV0°1450 

C 

EV0P146U 

0087 

37 

1 A*  IA ♦ 3 

EV0P1470 

0066 

W RITEI6, 6615 )EF FAB 

EVOP 1480 

0U6S 

38 

I A*  I A ♦ 1 

EV  OP  1490 

C 

EV  IP1500 

C 

RECO^MEnD  THE  NEW  OPERATING  CONDITIONS 

EV0P1510 

C 

EV0P1520 

C 

1FE  :IASE  PuINT  FOR  SHIFTING  OPERATIONS  IS  TAKEN  TC  BE  THE 

EV0P1530 

f. 

LC«cK  LEFT  PCINT.  ALL  CTFEF  POINTS  ARE  DETERHl NEC  AS 

EVOP 1540 

C 

UI SPLACE^EN 1 S 7 F OF  THIS  PCINT.  SHIFTS  TO  A NEW  SET  OF 

EV  OP  1550 

C 

OBSERVATIONS  ARE  CETLRN INOU  BY  SHIFTING  THE  BASE  PGINT 

EVJP1560 

C 

AND  RECOMPUTING  THE  OTHER  THREE  POINTS  FRCM  THE  NEW  BASE. 

EV^P  1570 

C 

EVOPl 580 

OOSC 

GO  FC  (99c’,lo0l«lC02,lUl2, 1013,1023, 1123), IA 

F VHP  1 590 

C 

EVOP 1600 

C 

NONE  UF  THE  EFFECTS  ARE  SIGNIFICANT 

EV0P1610 

C 

EVOP 1620 
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0051 

995 

WRI U (6,6617) 

EVOP  16  30 

OOS2 

GO  TO  lOoO 

EVOP 1640 

C 

EV0P1650 

c 

SPIED  IS  significant 

EVOP1660 

c 

EV3P1670 

0053 

1001 

1 H EHFA)20CC,2001,20Q1 

FV0P1680 

0054 

2000 

1SP«ISPH 

EVOP 1690 

C 

EVUP1700 

C 

DETERMINE  IF  THE  SELECTED 

POINTS  AKE  WITHIN  THE  0 EF IN  ED 

EV3P1710 

c 

ENVlKCHNtNT. 

EV0P1720 

c 

EV  OP 1 7 30 

005  5 

loco 

I F ( ISPU-NSPD)2033, 2033,2034 

EVOP 1 7 40 

0056 

2034 

I$P*ISP-1 

EV3P1750 

0057 

. 

WRI TE(b,6018) 

EV1P1760 

JJ5e 

20  3 j 

I F ( Ui;a-KFD)  2035,2C35I2036 

EVOP1770 

0055 

2036 

I FU-lFC-i 

EV0P1730 

01CC 

Wkl  It ( 6,6019) 

EVTPL790 

01  Cl 

2035 

1 FI  ISP)2J37, 2037,  2038 

EVOP l 800 

01C2 

2037 

ISP«ISPM 

EVOP 1810 

01C3 

wP  I Tfc  ( o , «jC20) 

EVOP 1 820 

01C4 

203d 

I F ( 1FO20 39, 2039, 2040 

EVOP1330 

0 1 C 5 

2CJ9 

1 FU* IFD*1 

EV  JP 1840 

01C6 

WRITE ( 6, o021 ) 

EVOP1850 

U1C7 

2040 

WRI  TL  ( 6 , 66  16)SPCFC(  1 SP 1 , S PE  EO  (I  S P ♦ 1 ) , FEED  ( I F0» 1 1 , FEE Ct  IF  OH)  , 

E V )P I860 

ISPECDt  ) SP)  ,SPEEU(ISPU)  ,F EECCIFD) , H EEC(  IFO) 

EV  TP  1 870 

C 

EVOP 1880 

C 

PLul  RECJPMENUEC  CPEFATING 

CONDITIONS 

EVOP1890 

C 

EVJP1900 

0 1 C 6 

CALL  OlPTS  (SPEED  (I  SP  ),SPEEO(  ISP  ♦!  ) , F£ED(  IFO)  ,FEECl  IFO  1)  , 1 GO  ) 

EVOP 19 1 0 

01C5 

ISP-ISPl 

EVOP 1920 

01  1C 

I FL)  * I FC 1 

EVCP1930 

0111 

l GO* I CC* 1 

EV0P1940 

01  12 

I r ( | GO • E 3 • 2 )G0  TO  145 

EVOP 1950 

01  13 

RETURN 

EV0P1960 

01  14 

20C1 

I SP -ISP-1 

EVOP 19  70 

01  15 

GO  Tu  1000 

EVOP 1980 

C 

EVOP 1 990 

C 

FEEL)  IS  SIGNIFICANT 

EV0P2000 

c 

EV0P2010 

01  16 

1002 

IK  EFFC)2CJ3,2004,2C04 

EVTP2020 

01  17 

2003 

IFO»IFOH 

EV0P2030 

oi  ie 

GO  ro  IGOO 

EVOP2040 

01  15 

2 GC4 

I FDMFO-l 

EV1P2050 

012C 

GO  TO  1000 

EVTP2060 

C 

EV0P2070 

C 

SPEED  C FEED  OP  THE  INTERACTION  ARE  SIGNIFICANT 

EVOP2080 

C 

EV0P2090 

0121 

1012 

1F( AHS(EF FA) -SELVES ( IG J) )23C5, 

20C6,  200 o 

EV0P2100 

0122 

2005 

YH1N-AMDK  Y(  IU0,NUE1),Y(  1GC,N0DIM)  , Y ( 1 GO,  NCE  1 *2  ) , V ( 1G0,NUP$  )) 

EVOP2110 

0123 

UO  2U07  1*1,4 

EV0P212O 

0124 

l F(  r(  lGu,NJbl*I-l  ).EC.YM1N)GC 

TO  2008 

EVOP  21 30 

0125 

2007 

CON  riKUE 

EVOP2140 

01  26 

2006 

GO  TU  (2JC9, 2010,  2011, 2012), I 

EV0P2150 

0127 

2009 

ISP*ISP-1 

EV0P2160 
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0128 

IFD«IF0-1 

EV3P2170 

ul  25 

GU  TO  10J0 

EVOP2100 

0I3C 

2010 

I SP*I  SP*1 

EVnP2l90 

0131 

I FU*lFO»l 

EV0P2200 

0122 

GO  TU  1000 

EV3P2210 

0123 

2011 

ISP»13P*1 

FV3P2220 

0135 

! Fl>*  I FD— 1 

EV0P2230 

012* 

GO  ro  1000 

EV1P2250 

0136 

2012 

ISP-ISP-I 

EV0P2250 

0137 

IFUMFDU 

EVOP2260 

0136 

GO  TO  1C0C 

EVDP2270 

0135 

2006 

i M E*FFA)2013,  2015, 205 

EVUP22BO 

0 1 AC 

2012 

ISPMSFH 

EV3P2290 

0151 

GO  TO  2015 

EV0P2300 

0152 

2015 

ISP*  ISP- 1 

E VnP23 1 0 

01  52 

2015 

1 F ( cFF:))2JU>,2Gl7,2Cl7 

EVGP2320 

0155 

2016 

I FU» 1 FU» l 

EV3P2330 

0155 

GO  TC  lUoO 

E V0P2  350 

0156 

2017 

lFO*IFn-l 

FVJP2350 

0157 

GO  TO  1330 

EV0P2360 

C 

EV0P2370 

C 

S P E 1 3 ANl>  1 ME  INTERACTION  ARE  SIGNIFICANT 

EV0P2330 

C 

EVOP2390 

0158 

1013 

I F ( EFFAJ201U, 2GIS,2C15 

EV IP 25 00 

0155 

2018 

ISP«I$P*l 

EV0P2510 

01  5C 

GO  TO  2C2C 

EVOP  2520 

0151 

2015 

ISPMSP-l 

FVOP2530 

0152 

I F( Y( ICC,  NOBS >-Yl  IGC,  NObl  1)  2021,2021,  202  2 

EV0P2550 

0152 

2021 

l Fl)*IF0#  1 

EVOP2550 

0155 

GO  TO  ICOC 

EV0P2560 

0155 

2022 

i FU  * I FC-1 

EV0P2570 

0156 

GO  TO  ICOC 

EV0P2580 

0157 

2020 

IF(  ii  ICU.NJdi  ♦ 1 1— Y ( IGC,  NOU  l ♦ 2 ) 12021, 2022,  202  2 

E V 1P2590 

C 

EV3P2500 

C 

FEEO  AND  THE  INTERACTION  ARE  SIGNIFICANT 

EVOP  25 10 

C 

EV0P2520 

0156 

1023 

I F( 6FFB)2023,2o25,2C25 

EV0P2530 

0155 

2023 

1FLWFIUI 

CV0F2550 

016C 

GO  TO  2025 

EV0P2550 

0161 

2025 

1F0-IFC-1 

EV0P2560 

0162 

i PC YC IGC,NJbL  *Z 1-Y ( IGC, NOB  1 ) ) 202  7, 2026,2026 

EV3P2570 

0163 

2G26 

I SP»ISP-l 

EVOP  2580 

0165 

GU  TO  UJO 

EV0P2590 

0165 

2027 

ISP-ISPU 

EV3P2600 

0166 

GO  TO  IC3C 

EV0P26 10 

0167 

2025 

i FI Y I IGC,  H Obi  *1 1 — Y 1 IGC, NOBS )) 2027,  2026, 2026 

EVJP2620 

C 

EV0P2630 

C 

ALL  1 ERMS  APE  SIGNIFICANT 

EVQP2650 

C 

EV0P2650 

0168 

1123 

I FC  EFFA)2028,2029,2C25 

EV0P2660 

0165 

2028 

ISP*ISP»l 

EVOP2670 

01  7C 

GU  TO  2C30 

EV0P2680 

0171 

2029 

ISP*ISP-l 

EV0P2690 

0172 

2030 

I F ( EFFd)2031,2032,2C32 

EV0P2700 
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0173 

0174 

0175 

0176 


0177 

0178 
0175 


oiec 


0 1 E 1 


0 1 C 2 


0163 

0164 


Oitil 
01  kt 
0167 
0166 
0185 
0190 
0151 

0192 

0192 

0154 

0155 

0156 
0197 


0156 


2031 

2032 

C 

C 

C 

148 

150 

6018 


6019 

6020 

6021 

0553 

6600 

6601 
66  02 
66C3 
OtfC4 
6oC5 
6606 
66  C7 

6010 

6611 

6612 

6613 

66U 

06  15 

6616 


IFD*IFDU 
GO  TO  ICOO 
I F0«IF0-1 
GO  TO  1000 

FORMATS 

FORMAT ( • 1 • , T46, 'CCST  aN/LYSIS* ) 

F Oh  mat  ! * 1 ' *T46,  ' PRliCUCT  ION  PATE  ANALYSIS  * I 

FUK  UT(  • JThE  RECCFmCNDL'D  SPEED  HAS  EXCEEDED  THE  MAXlMuN  DEF  INCD  * 
l'SPELU  LIMIT* /'O' • lOx, *TmE  PkEVIJUS  SCTTINGS  SHOULD  0 E PE-fcXAMiN* 
2* ED* I 

FuRM-l < • 7THL  RlC(JKM£NCED  FEED  HAS  EXCEL  DEO  THE  MAXI  MUM  DIF  INED  ' 
1 'FELL)  LIMIV/'O*  fl0X,*7HE  PREVIOUS  SCTTINGS  SHOULD  3E  RE-bXAMIN* 
2'EO*  ) 

F Uk  T ( 1 DTHE  kECUMM  ENDED  SPEED  HAS  FALLEN  BELCVi  THE  MINIMUM  *, 

♦ •UEPIMJ  ', 

I'SPLED  LIMIT* /*0*  , 1GX,*1HE  PREVIOUS  SETTINGS  SHOULD  6E  FE-EXAMN* 
2' ED*  ) 

FORMAT! 'OTHL  F.ECURM  ENCEO  FEED  HAS  FALLEN  OELCW  TEE  MINIMUM  *f 
♦'DEFINED  *, 

1 'FEED  LIMIT1 /» C'  ,10X,* THE  PREVIOUS  SETTINGS  SHCULD  BE  PE-fX/VIN* 
2'  ED*  I 

FORVkK*  j*  f TiOt  •SPEED  ( PPM)  • , t 25,  4F  1 0. 0,  10X,  4F  1 0 .0/ * O'  , T 10  , *FCEO 

♦ , • ( 1PR ) ' ,Tc5,4F  1C. 4,  10X, 4F lu.4) 

FUP  M#  T C • 1'  fT42t  * EVOLUTIONARY  OPERATION  ANALYS  IS  * , /T  4 t 9 
♦'PHASE  ' » I J» 

1'  CYCLE  '» 13//'C' ,T46f 'CALCULATION  OF  AVERAGES  ' /' 0 • , T45* • COST ' » 
2T9Ut  * PkuJtCT  ILN  RATE'/  'JTPERAr  INC  COnO  IT  IONS  ' ) 

FOR  1A  T l / ' OPnE  V I OUS  CYCLE  SUP  * , T25, '.F  1 0.  2 t 10X,  4F  10  .4  | 

FUn  M/.T  ('  JPkcV  i JUS  CYCIC  AV  fc  R AGE  * , T 2 5 * 4F  1 0 . 2 » iCX  , 4F  l D . 4 J 
FUR  MAT!*  3NtW  UtlSERVATIUNS*  , T25 , 4 F 1 0 . 2 f lOX , 4F  10.4  1 
f Ok  1A1  ( ' JGlFr  CS ENCE5 ' , T 25 , 4F 10 . 2, 10X , 4F 10 . 4 ) 

I- On  MAT  ( ' JNfcrt  SUMS  * fl  25,*»F  ID.  2,  l OX,  4F  10.4  J 
F UR  m AT  ( ' J.Mt  n AvErAGES  : Y ! I ) ' , T 25  * 4F  1 o.  2,  l OX f 4F  1 0 . 4 I 
FORm/,1  !//«  j*fTA6f  *CALCJ|/TICN  OF  2 S.C.  L I M I T S : ' / / ' 0 • , 

♦ •FUR  Nt  *»  EFFCCTS'  ,T36,F10.4  ,T85,F10.4) 

FURIMl'OlHt  EFFECT  OF  SPEED  IS  ' , T 5 3,  F 10 . 4 , • 


*• SI  GN IF IC  ANT 
FORMAT! • JtHE 
♦, ' — ♦*♦' ) 
FUn  VA 1! • UTHE 

♦ f * •) 

FORMAT if UTHE 
* , * — • } 
FORMAT ! ' J f HE 

1F1U.4,1  

FURmaT! ' UTHt 
IF  10  • 4 f • *♦♦ 


. > 

EFFECT 

EFFECT 

EFFECT 

effect 


SPEED  IS  ' , roatFio.4f • 
FEED  IS  * »T53»  FiO • 4* ' 
FEED  IS  ' ,T88, F10.4, ' 


NOT 

♦•♦  — S I GN! I 


UF 
CF 
OF 

OF  THE  INTERACTION  OF  SPEED  AND  FEED  IS  ' 


FI  CANT 

NTT  SIGNIFICANT 

♦♦♦  — SI GN IF ICANT* 


EFFECT  UF  THE  INTERACTION  JF  SPEED  AND  F EE  0 IS*,T88, 
— SIGN! FICANT  — *♦♦' I 
PUP  MAT ! // / 1 0*  1 1 1 5 1 * PECOMME  MC  ADDITIONAL  OBSERVATIONS  BE  T*kfcN  AT 
1, 'PCINT^:  '///'  ' , T 1 6t  * 4 J SPEED  * ' , F 1 0. 0,  T4  5 , * 2 ) SPEEC  » *,F10.0/ 

1'  ' , T 1 5 * * FEED  » * ,Flu.4,T45, * FEED  - ' , F 10 .4// ft  *0  * , 

1115, '1)  SPEED  » * , F 10.0,1 45, • 3)  SPEED  * *,F10 .0/'  ', 
l*  ' ,T  1 5, • FEED  * •, FIO*  4, T45 , ' FEED*  '.F10.4I 


EV0P2710 
EVDP2720 
EV0P2730 
EV1P2740 
EV0P2  750 
EV0P2760 
EVOP2770 
EV0P27  80 
CV3P2790 
, EV  )P2800 
, £ VOP28 l 0 
EV'*1P2820 
, FVOP2830 
,EVUP2340 
EV0P2850 
EV0P2860 
EV. TP 2 8 70 
, tV  UP2080 
C V"'P2890 
FV0P2500 
EV7P2910 
, EV3P2920 
EV0P2930 
' E V0P2940 
EVUP2950 
EVJP2960 
EV0P2970 
EV0P2980 
EV0P2990 
EVUP30G0 
EVJP3010 
EVHP3020 
EVUP3030 
EVUP3040 
EV JP3050 
EV0P3G60 
EVOP3070 
E VOP3080 
EVJP3090 
•FV0P31O0 
EV0P3110 
•EVUP3120 
EV0P31 30 
EVUP3140 
EV0P31 50 
EV0P316G 
EVOP  31 70 
FV0P3 l 80 
FV0P31 90 
•EV0P3200 
EV0P3210 
EV0P3220 
EV0P3230 
EVOP  3240 


107 


FORTRAN  IV  G LEVEL  21 
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OIS9 


O2C0 

02C 1 
02C2 
02L3 
0204 
02  C 5 
02C6 
02C7 


6617  F0hvAT('3NJ  PARAMETERS  AR  E S IGN ! F IC ANT • / • 0* , *THC  PREVIOUS  SFT1 I NG • E VOP 3250 


It  * S SHUULO  BL  KE-CXAVINCO* I 

6620  FuR-Al  (//•;>•,  T44, 'CALCULATION  OF  STANDARD  DEVIATIONS*//) 

C 6621  FOrMAU  •0*,'PREV1CL$  SUM  S * ,T3 5,  F 10 . 4,  T 35,  F 10 . 4 ) 

6622  FOP  VAT ( * JPkEV  I'JUS  aVEFaGF  S * , T35, F 1 0. 4 , T85,  F l C. 4 ) 

6623  ruM*AT  PJ.NF*  S « RANGE  * F4  , N»  , T3:>,  F 10. 4,  T85,  F 10.4) 

6624  FOR  I • ORANGE  * ,T35 ,F 10.4 , X 85, F 10.4 ) 

6625  FURM/T  (*ONbw  SUM  S*  ,T35,F 10.4,Td5,F 10.4) 

6626  FUPMaTI*  JNth  AVERAGE  S = NER  SUM  S/2  * ,T 35, F10 .4 , T85, F 10 .4  ) 

7000  FORMAT l * J. i2345*/6E13.8/6El2.fl/6ELJ.0/2El3.8,4I2l 

END 


FVJP3260 
EVDP3270 
EV3P3200 
FVJP3290 
EVJP3700 
EV0P3310 
EVUP3320 
EVJP3330 
EV0P3340 
E V0P3350 
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00C1 

sueRcurnE  orgpis  (sfi,$P2,fli, ro2,  igo? 

OGP 

10 

C 

ORGP 

20 

C 

SU3RLUT  1 NE  UKGPIS  FLCTS  EACH  OCFINED  FEED-SPEED  CO '*8  IF.  AT  ION 

ORGP 

30 

c 

AS  **  LN  A GFAPF  AND  INC!  CAT  E3  THE  FOUR  PL  I NT  6 OF  PARTICULAR 

ORGP 

AO 

c 

INTEREST  WITH  / N •*' 

OFGP 

50 

c 

CRGP 

60 

00C2 

CO  V-VCN  /:  101/5PEEC!  100)  , r EEC ( l 00 1 , oPND C 100)  ,NSPD,NFD,  ISP,  I FO 

CRGP 

70 

0003 

CO  1KCN  /Cl J5/PT1, PI  2, PT3, NAME C 501, CHART (ICO, 501 

GRGP 

80 

OOCA 

01  HE.mS  1C  i FND150) 

JR  ZP 

90 

0uC5 

00  6 1*1, dC 

ORGP 

100 

00  C 6 

6 

FNU< l > »0.0 

URGE 

110 

UJC7 

NS*  10C/NSPI) 

ORGP 

120 

ooct 

NF*  50/NF 0 

CRGP 

130 

OUCS 

NSl -NS^NSPO 

CRGP 

1AU 

00  I C 

NF  1 a,NF*,\FD 

ORGP 

160 

0011 

ISP  1*1 

C«GP 

160 

GOU 

1SP2*1 

09  GP 

170 

0U13 

N FD  1*  1 

ORGP 

180 

OOIA 

N FL)  2*  1 

CRGP 

190 

00!t 

UO  l li*l,NSPO 

ORGP 

200 

00  It 

1 F 1 SPl  .t  J.SPf  C LU  S)  IISPl»l  1S-1  )*NS*1 

C»GP 

210 

0017 

IH(  SE2.E  J.SK ted  SIIISP2*  C IS-l  )*N5*l 

CRGP 

220 

ooie 

1 

u on  r i i\ue 

ORGP 

230 

0019 

00  2 JS*  L , NFO 

ORGP 

2 AO 

002C 

1 M HLl.E  J.FkEUC  JS)  INFCl-l  JS-l  MNFfrl 

ORGP 

2 GO 

0021 

I f ( FU2  • t ?.FtEO( JS) )NF02-( JS-l)*NFf 1 

C°GP 

260 

0022 

2 

CONTINUE 

CROP 

270 

0023 

K * 0 

ORGP 

280 

0024 

OU  3 l*l,NFl,NF 

CRGP 

290 

0025 

K *K  ♦ 1 

•OR  GP 

300 

0026 

3 

F NU ( I 1 -FEEU(K) 

ORGP 

310 

C 

UU  J J*l,NSl,NS 

ORGP 

320 

C 3 

LHA  ITCJtl )*PTl 

CRGP 

330 

0027 

CM/  KT(  IS  P l , NFD1  )■  PT  2 

ORGP 

3 AO 

0026 

CHART ( I SP  1 , NFU2 ) * PT2 

ORGP 

350 

0029 

CHAKTC 1SP2,NH01 ) * PT2 

ORGP 

360 

00  3V. 

CHART! 4SP2,Hr02l*FT2 

ORGP 

370 

0031 

GO  TO  (31  , 32,^3), IGC 

ORGP 

380 

0032 

31 

WRITE (o,6A) 

CRGP 

390 

0033 

GO  TO  3A 

CRGP 

AOO 

003  A 

32 

r»R I Tt  ( 6,65) 

OKGP 

A10 

0035 

GO  TO  3A 

CRGP 

A20 

0036 

33 

wRI  T E ( 0 , v*) 6 ) 

ORGP 

A30 

0037 

3 A 

I *50 

ORGP 

A AO 

0036 

A 

IHl FNC ( I ) • EO • 0 • 0 ) GO  TO  5 

ORGP 

A 50 

0035 

wRITE(o,o7|NAPk(l  ),FNU  (I)  , ( CM  j (j  f ||  , J- l , l OC  I 

ORGP 

A60 

00  AC 

GO  TJ  8 

ORGP 

A 70 

OOA  1 

5 

WRI  rt  (6,69)NAHC ( I ),  (CHART ( J , I ) , J* l , 1 00 ) 

CRGP 

AOO 

OOA2 

8 

1*1-1 

OP  OP 

A90 

00  A3 

I F ( I.GT.OJGO  TO  A 

CRGP 

500 

00  A A 

K»0 

ORGP 

510 

00A5 

00  11  1-1,10 

0«JC,p 

520 

OOA6 

DU  LI  J«I,NSPO,iO 

ORGP 

530 

OOA? 

K-K  ♦ L 

OP  GP 

5 AO 
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OUC  • 74045 


14/G4A4  7 


0048 

11 

SPNCl  JJ*SPkEDU) 

OPGP  550 

0C4S 

0*176(6.70 

ORGP  560 

00*0 

*KI 1 T (6. 7 i ) ( SPNU ( I ) . 1*1. NSPD) 

ORGP  570 

00*1 

*kl7E(6. 721 

ORGP  580 

00*2 

CH6»TIISP1.NFD1  J-PT3 

ORGP  590 

00*3 

C HA  *T ( ISP i.NF  D2 )■ PT 3 

ORGP  600 

i*0  54 

CHAkT( ISP2.NF01 J-FT3 

ORGP  610 

005* 

CH*e)(I$P2.NFD2)«PT3 

CRGP  620 

0056 

RETURN 

ORGP  630 

C 

ORGP  640 

C 

FCRMAlS 

CRGP  650 

C 

ORGP  660 

0057 

64 

ECKVAT  ( • l • »T54.  * N E W CCST  ENVIRONMENT  • ) 

ORGP  670 

00*8 

65 

PukM/T ( • 1 • .T54, *NtW  PRODUCTION  RATE  ENV I RCNMENT • ) 

OPGP  680 

00  59 

06 

FOR  PAT ( ' l • .7*4. • S Pt C I FlfcU  ENVIRCN1ENT •/•  ') 

ORGP  690 

006C 

67 

FOR  T ( • •.  1A1.FIC.4,  2X,  • I • .2X,  100AU 

CRGP  700 

0061 

69 

F OF  VAT ( • '.1A1.12X. >I'»2X.100A1) 

ORGP  710 

0062 

70 

FORMAT! • • . 13X, • J •/ 

ORGP  720 

♦ ' • « 13X. • I 

— • . ORGP  730 

•) 

ORGP  740 

0063 

71 

FORMAT ( • ••7X.10F10.C) 

ORGP  750 

0064 

72 

FORmat(»  • . T60. • S P E E 0 ( PPM  J • ) 

ORGP  760 

0065 

END 

ORGP  770 
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UOCi 

SUBROUTINE  PLUTPM(Il) 

PLOT 

10 

c 

PLOT 

20 

c 

SlbR3UTINE  PLC1PM  TABULATES  THE  RfcSPfJNSFS  VERSUS  FEEDS 

PLOT 

30 

c 

AND  S*>EbDS. 

PCUT 

AO 

c 

PLOT 

50 

OOC? 

CU4MUK  /Cl02/IPH»lCYC»SCYCLE»RL0,X(2Af 100lf Y(2» 100) 

i NOBS  * NVB1 

PLOT 

60 

•t  lUP,BSP»aF0,HlLC,TLC(20) 

Pl.OT 

70 

OOC3 

I 2*  1 1 ♦ 1 

PLOT 

BO 

UOCA 

I 3*  C l ♦ 2 

PLOT 

90 

uOC5 

i A*  1 1 ♦.> 

PLOT 

100 

occt 

wM  TfcU»  72) 

PluT 

110 

OOC  7 

whlJtiof71)X(2fI<i)*Y(ltIA)fYClfI2)fYI2fI^)»Y(2fI2)fX<2»li)f 

PLOT 

120 

• Y(lfn),YU»tJ)»Y(2tm»Y(2,l3)»X(l,m»X(l»l3l 

PLOT 

130 

ooca 

KE1  JR  N 

PLOT 

1 AO 

c 

PI  uT 

150 

c 

F C R N A*i  S 

PLOT 

160 

c 

PLOT 

1 7 J 

oocs 

71 

FOKNATCO't  JOX,FlC.Af  • COST  * • » F 10 . 2 t 1 CX  , • 

COST  c ' 

'#  PLOT 

1 Bo 

♦F10.2 

PLUT 

190 

•/'U'r^jA, • PfcUi).  FI.  « • »F  10.A» 10A» • PRCQ.  R1  . » 

• * F10. A//// 

PLOT 

200 

*• 0*  »3QxtElu.At 1 CCST  * ' « F 10. 2 1 10X  * • CCST 

* •tFlO.2 

Pi  OT 

210 

* / * 0 * t ACX  f • PKJD.  n.  ■ • 9F  lJ.AtlOX,*  PPCD.  FI.  • 

• i F10.A/ 

PLOT 

220 

♦•w*  »5AXfF  U.0»2AXtF  lU.O/^O*  *69Xt  »S  PEED  (RPPP) 

PLOT 

230 

00  1C 

72 

FUF  ^TPISjJX,'  TABLE  OF  RESPONDS  ,///»0*  »33X»* 

FEED  IIPPJM 

) PL  OT 

2 AD 

0011 

END 

Fl.QT 

250 
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OOCl 

SUb*UUT  ME  RtGPLT  CbC,  e,  INC,  I PEG,  IGO) 

PEGP 

10 

C 

RCGP 

20 

c 

SLHRJUTIME  REGPLT  PREDICTS  THE  PERFORMANCE  INC1CIES  FOF 

PFGP 

30 

c 

FELDS  AND  SPEEDS  ADJACENT  TO  THE  ObSE R VA 11 ONS • 

PEGP 

50 

c 

FEGP 

50 

00 1 2 

CO 4 MON  /C  LOl/SPkECilOOl  ,F  EECUOO ) , SPhC(  100)  , NSPO , NFO,  ISP,  1F0 

RCGP 

60 

00C3 

D I MENS  It-  1 I.NU(4],IJ(2,5),SPDUI2,lOO),bO(21 ,IRGG(2) 

REGP 

70 

00C4 

DAY  A PTTl , PT'I 2/ *(*,•)  ’/ 

RFGP 

80 

00C5 

wRI TE(6, 71) 

REGP 

90 

c 

PFGP 

100 

c 

CCMPUTE  PHtOI CTfcO  VALUES 

FEGP 

110 

c 

REGP 

120 

OOCfc 

81 

ISP  1*1  NLA  1 1 1 — t 

REGP 

130 

00  C 7 

1 F ( ISPl.'.E.OJ  ISPl«i 

REGP 

160 

OuCtt 

1 SP  2* I M)  ( 2 ) ♦ l 

FEGP 

150 

U0C5 

i K<  ISP2.GT  .NSPO  11  5P2*NS PL) 

FEGP 

160 

00  1C 

IFOUIM  (3)-l 

REGF 

170 

00  ll 

l F ( IFC1.LE.0)IFDI*1 

RFGP 

160 

oO  1 1 

IFU2«IND(4)U 

FEGP 

190 

0013 

I F ( IF02  • iT  •NFO) I F C2*NFD 

FEGP 

200 

001* 

I «I FD2 

REGP 

210 

00  15 

505 

c*a '.uG ( feed (i ) ) 

FEGP 

220 

0016 

UQ  506  J*  ISPl * I SP2 

REGP 

230 

00  17 

A«ALCG( i »LcO( J) ) 

PEGP 

250 

ooie 

DO  sOo  <*1,2 

REGP 

250 

OOIS 

IRt  «IREGU) 

RFGP 

260 

0O2G 

SPOJ(K,J)  « bO(K)  ♦ H C K,  1 ) * A ♦ 0 ( K, 2 ) *C 

REGP 

270 

0021 

bu  TO  (500,511,512,513,515),  IRE 

REGP 

280 

Oo  22 

511 

SPJC(K,J)  « *»PUD(K,J)  ♦ b ( K , 3 ) * A*C  ♦ e(K,5)*AAA 

PEGP 

290 

o 023 

GU  fU  5J& 

REGP 

300 

00  24 

512 

SP)C(<,J)  * SPOO(K,J)  ♦ B ( K , 3 ) *A*C  f B(Kf6)*C*C 

FEGP 

310 

00  25 

GU  TO  506 

REGP 

320 

0026 

513 

SP)C(K,J)  » S POO  ( K , J ) ♦ (HK,J)*A*C  * B (K  | 5)  * A*A  ♦ B(K,5)*C*C 

FEGP 

330 

0027 

GU  1C  506 

REGP 

350 

0028 

515 

SPOl)(K,J)  - S POO  ( K , J ) f P(K,3)*A*C 

RFGP 

350 

00  25 

506 

COfrfTlMjc 

PfcGP 

360 

C 

REGP 

370 

C 

PLUT  PREDICTED  VALUES 

PEGP 

380 

C 

REGP 

390 

003C 

ID« ISP2-ISPlf 1 

REGP 

500 

0031 

1 F ( IC.GT.IOIGO  TO  11 

RFGP 

510 

0032 

wPI  rt(o,  72)  Ft  EDO  ),  (SPJDU,  J) , JMSP1,  ISP2) 

PEGP 

520 

0033 

WRI  TElo, 78)(PTi  l, SPC0(2, J), PTT2, J»I SP1, IS P2) 

REGP 

630 

0035 

GO  TC  77 

PC  GP 

660 

0035 

11 

1ST-ISPI 

REGP 

650 

0036 

33 

I fchC*  1ST  »9 

PFGP 

460 

00  2 7 

IfcNC**I\J(ICNU, ISF2  ) 

RFGP 

470 

0038 

I F ( 1ST  .L0.1SP1  ) WRITE  16,  72)  F EEl)(I  ),  ( 5PDU(  l , J)  r J»!ST,  IE  NO) 

REGP 

480 

0035 

1 F(  1ST  .NE  . IS  PI  IaHT  t (6,  62  ) ( SPDD(  l , J ) , J«ISY,  IE  ND  ) 

REGP 

490 

004C 

wP  I TCCo,  78  M PTTl  , SPCC(2,J  ) , PTT2,  J-IST  • I END) 

REGP 

500 

0051 

IFC  IST.Cfc.  ISP2JG0  TC  77 

REGP 

510 

0052 

I ST  *1  S I ♦ 1 0 

REGP 

520 

0053 

GO  TO  33 

REGP 

530 

0055 

n 

I-I-i 

PEGP 

540 
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00  A 5 

1F(  1.GE.IF0UG0  TC  505 

PEGP 

550 

0046 

-WR!TE(6t  70 

REOP 

560 

0047 

WRITE  16, 73)($PEEUU)»  I«lSPlfISP2l 

REGP 

570 

0048 

I F ( I U • GT . 5 )G0  TC  98 

PFGP 

580 

004$ 

WR1  TE ( 6f 74) 

REGP 

590 

u05C 

99 

CALL  EXIT 

PFGP 

600 

0061 

96 

WRI  rC(6»75) 

»C3P 

610 

0062 

GO  TO  99 

PFGP 

620 

C 

kCGP 

630 

C 

FORMATS 

PEGP 

640 

C 

PFGP 

650 

C 

KCGP 

660 

0053 

70 

FORMAT!'  1 « 12X9 ' 1 '/ 

PEGP 

670 

• ' *912X9*1 • , 

KCGP 

680 

CEOP 

690 

0054 

71 

F0RMAU*l*tT5J,  *PKtCICTEO  COSTS  •/•  • fT44 ,»< PREDICTED  PRODUCTION* 

PEGP 

700 

• •*  WJT  tS I * / / • 0 FEED  1 1 PR  ) ' / • 0 * ) 

kFCP 

710 

0055 

72 

F UK  4AT ( * * , 12 X, • I */•  • , F 1 0 . 4 , 2X , • l • 9 2 X, 1 OF 1 0 . 2 1 

PEGP 

720 

00  56 

73 

FOts"*T  (*J*,15X,10FIC.C) 

PEGP 

7 10 

0057 

74 

FOR  -AT(*J*9T3Jt*S  P fc  E 0 <RPH)*I 

RFGP 

740 

Uu58 

75 

FORMAT!* J*rT609 *S  P t E 0 <RPM)») 

PEGP 

750 

005$ 

7 B 

FURRAT!  • *,12X9  'I*93x9lO(2A,Al  ,F6.4fAl  )) 

PtGF 

760 

0060 

82 

FUk  VA  T ( * ' 9 12X v ' I •/•  12 X, *1* ,2X,10F10.2) 

FCGP 

770 

0061 

END 

PFGP 

780 
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ouci 

SUbRuUT I 4E  KfcGRES  ( X, Y,N, IREG, 8SU30, HHAT, IDEP  ) 

REGR 

10 

C 

REGR 

20 

C 

SLbNJLTIuE  P.EGKES  CCMPJTF5  t 

(V  MULTIPLE  LINEAR  REGRESSION  ON* 

REGF 

30 

C 

M VAXIAULES  HlTh  N CjBSLRVATIOIS,  EQUATIONS  UStO  ARE: 

REGR 

40 

C 

REGR 

50 

c 

PI  • BC  ♦ bl  LNCSPE'fO)  ♦ 02 

LN( FEED) 

REGR 

60 

c 

REGR 

70 

c 

PI  « DO  ♦ Bl  LN (SPEED)  ♦ 02 

LN(  FEED)  4 03  IN(SPfcEC) 

LN( FEED) 

REGR 

80 

c 

PF.GP 

90 

c 

PI  * 00  * Ul  LN (SPEED)  ♦ B2 

LNI FLED)  ♦ B 3 LN(SPEEC) 

IN(FCED) 

REGR 

100 

c 

♦ B4  LM  $ PE  ED  ) * *2 

REGR 

110 

c 

REGR 

120 

c 

PI  * 130  ♦ Ul  LMSPErDJ  ♦ 02 

LNCFCEU)  ♦ 03  LN(SPFEC) 

LNCFEFU) 

RFGR 

130 

c 

♦ H4  LMFEEU)**  2 

REGR 

140 

c 

REGR 

150 

c 

PI  * uC  ♦ bl  IN (S PEED ) ♦ 02 

LN( FEED)  ♦ HJ  LN(SPEEC) 

LN( FEED) 

REGR 

16C 

c 

♦ 64  IN ( SPELD ) * *2  ♦ 

ti5  LN(FEE0)**2 

REGR 

170 

c 

REGF 

100 

c 

FEGR 

190 

00  C2 

I) I MENS  16  J l*H/iH5)  ,Y8AR(  2) 

REOP 

200 

00C3 

01  MENS  IE  1 X ( 2 , 10L)  , > (2  , 1 00  1 , 0 JM  l ( 5 , 5 ) , XP  X ( 5 » 6 ) 

RFGR 

210 

OOC4 

uIMENSIl.  1 XBAk  ( 5 ) f 

RAC(5) ,ST0(5),T(5)  ,YMAT(100), 

REGR 

220 

•Oil'  F(  lb 01  , XPYC  2*  5 1 1 SST  12) 

REGR 

2 30 

ooc  5 

0 I MENS  1 L 1 XPXC(5,6) ,XFYC<2,5) 

REGF 

240 

0 JC6 

NP  • N 

REGR 

250 

00C7 

EPS  ■ 1.0E-05 

REGR 

260 

ooce 

50 

H * 5 

REGP 

270 

00C5 

XN  * N 

PC  G*> 

280 

001C 

A?  * P*l 

RFGR 

290 

0011 

IFC lOEP.EC.ilGU  TC  30 

REGR 

300 

0012 

DO  10  1*1,5 

REGR 

310 

0012 

DO  10  J* 1 , 5 

REGP 

320 

00  1 A 

lu 

AHXCI,J)*UUM1(1,  J) 

REGR 

330 

0015 

DO  11  1*1,5 

PEGR 

340 

00  16 

11 

xPxd  ,6)*XPY(2,I  1 

REGR 

350 

0017 

GO  TO  13u 

REGR 

360 

C 

REGF 

370 

c 

COMPUTE  THE  MEANS 

Pf  GP 

380 

c 

REGR 

3 90 

0016 

30 

DO  40  J * l ,NP 

REGF 

400 

OUiS 

X ( 1 , J)  * ALOG(X( l,J|) 

REGR 

410 

002C 

X ( 2 ,J)  * ALuG (X ( 2 , J ) ) 

RFGR 

420 

0021 

XC J,j)  * X( i, J)+X(2,J) 

REGR 

430 

o022 

X(<t,J)  * Xt l,J)*X(l,J) 

REGF 

4<»0 

0023 

X( 5 , J)  * X(2, J) *X(2, Jl 

FEGF 

450 

0024 

40 

CONTINUE 

REGP 

460 

0025 

DO  60  J* l , 2 

REGP 

470 

0026 

YBA^CJ  )-0.0 

RFGF 

480 

0027 

Du  60  1*1, N 

RFGR 

490 

0026 

bU 

YdA-UJ  )**0AMJ  )*YCJ  ,I)/XN 

PCGR 

500 

0025 

DO  JO  J*l,M 

REGF 

510 

003C 

XBARtJ)  • 0.0 

RFGR 

520 

0031 

DO  7C  1*1, N 

REGR 

5 30 

0032 

70 

X0AH  (J)*XBAk(J)+X(J,  11/ xn 

REGF 

540 
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0033 

80 

CONTINUE 

REGP  550 

C 

REGR  560 

C 

COMPUTE  THE  SLMS  GF  SQUARES 

RCOP  570 

C 

PFGR  580 

0024 

OU  90  1-1,2 

KCGP  590 

0035 

S ST (1  )»0.0 

PCGP  600 

0036 

00  90  J«l,N 

PFGP  610 

0037 

90 

SST  ( I )«SS  f 11)  ♦ ( YC  I,  J)-YBAR(  I MMY!  I , JI-YfMftC  1)1 

REGR  620 

0038 

DO  110  J - l , H 

«M.GR  630 

0035 

00  110  I - l » M 

REGF  640 

004C 

apa u,n  « o,o 

RFGP  650 

00  41 

00  IJO  K-1,N 

RFQR  660 

0042 

100 

aPA  ( J , i )aXPAl  J,1  ) MXI  I,K)-X8AR(I  ))MX(J,K)-X3AR(J)) 

PEGF  67  o 

0043 

APXCl  J,  I ) * A PX ( J , I ) 

REGR  680 

0044 

110 

CONTINUE 

F EOF  690 

0045 

00  13C  X* l, 2 

? EGF  7v/0 

0046 

00  130  J-1,M 

REGP  710 

0047 

APYCK  ,J)  - O.C 

RFGT  720 

0048 

00  120  I • l * N 

RE  GE  730 

004$ 

1 2 C 

X«>Y  (x  , J)»APY(K  , J)  * ( X ( J , I )-A8AR  ( J ) ) * ( Y ( K , IJ-YK/MK 

)) 

c FGF  740 

00  5c 

APYCi A,  J) - X PY ( A , J ) 

RE GF  750 

0051 

130 

CONTINUE 

REGF  760 

0052 

DO  1 4 C 1*1, M 

PEGR  770 

0053 

140 

aPX  1 I ,6  )*XPY( 10EP,  I I 

RF GF  760 

0054 

03  135  1-1,3 

REGP  790 

0055 

00  135  J - l , 3 

REGF  800 

0056 

O'JM  1 ( I , J)  - A PX  f I , J ) 

k fGP  810 

00  5/ 

135 

COf.  TINGE 

RE OR  820 

C 

REGR  830 

c 

DETERMINE  THE  FCRP  OF  THE  REGRESSION  EQUATION  ANO 

REGR  840 

c 

INVERT  i H E X»X  PATRIX 

PFGP  850 

c 

FEGR  860 

0056 

136 

GU  TC  (15  1,143,  142,  146),  IREG 

REGP  870 

0055 

142 

00  144  1*1,3 

REGP  880 

00  6c 

XPa  (4  , i ) *XPX( 5,1) 

REGP  890 

00  61 

144 

XPX ( I ,4) * X P X ( I ,5) 

REGR  900 

0062 

XPX  (<*,4)»XPX(  5, 5) 

REGP  910 

0063 

143 

rf-4 

PEGP  970 

0064 

00  141  I * 1 » M 

REGP  930 

0065 

141 

a PX  ( I , 5 ) - XP  X ( 1,6) 

PEGR  940 

0066 

C Ac  L ClMCcN(XPX,P,CFS,D,  1) 

PEGP  950 

0067 

IF  1/85(3)  • Gc  • EPS)  GO  TO  145 

PEGF  960 

006  £ 

GO  TC  143 

REGP  970 

00  £ S 

151 

M « 2 

REGP  980 

00  7C 

XPX  ( 1,3)  » aPX( 1,6 ) 

REGP  990 

0071 

X Pa ( 2 , 3)  « XP  a ( 2 , 6 ) 

RFGP 1000 

00  72 

Gu  TO  154 

PEGR1010 

00  73 

146 

CA*_L  CIMcGN(APX,M,EPS,U,l) 

RFGP 1020 

0074 

IF  (Atfi(  ))  .GE.  EPS)  GO  TO  145 

REGP 1030 

0075 

148 

IREG  » 5 

RFGP1040 

0076 

00  149  1-1,3 

RFGP i 050 

0077 

00  149  J-i, 3 

REGR 1060 

00  7 8 

149 

XPX ( I , J)-DUM1 ( I , J) 

FEGR 1070 

00  7$ 

M * 3 

REGP 1080 
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ooec 

00  15/  1*1,3 

PEGP 1090 

UU  8 1 

157 

XPX  l l , 5 ) 3 X PX ( 1,6) 

RtGRllOO 

U0«2 

155 

Call  DlPECMXPx»Mt£PS,0,  l ) 

P6GPI110 

C 

PFGPH20 

C 

COMPJTE  THE  REGRESSION  PARAMETERS 

REGP1130 

c 

hEGR!150 

0083 

165 

00  15C  l * l , M 

°FGC 1150 

0U65 

150 

OHAT  m»XPX(K,M»l» 

PEGF1160 

0085 

W • 0.0 

REGEll 70 

UU  66 

SSR  * 0.0 

PEGP11P0 

0067 

00  16C  t 3 1 , M 

PEG® 11 50 

00b8 

W * w ♦ 3HAU  I ) *XBAH ( 1 1 

REGP1200 

0065 

USUCU  « YHAR(  I0EP1  - W 

P.FGP1210 

005C 

1 6 w 

jSK*SSR  + eHAMIl  • XPY  (IOfcP,I) 

f CGF 1220 

0051 

XM  3 R 

PFGC 1230 

00  52 

XMSR  * SSR/XM 

REG*1  )250 

005  J 

>SE  * SSfUOcP)  - SSR 

FEGF1250 

0055 

X E 3 N-M-l 

KECP1260 

u05  5 

K * XC 

FEGP1270 

0056 

XMSE  3 SSE/XE 

HFGR1230 

0057 

h * X V$R/XMSC 

&FGH290 

Ouse 

oo  i/c  1*1, y 

PEGP 1300 

0C55 

KAim  ) * xpx(  i,n ♦xpse 

PEGP l J l 0 

OICC 

sruti  )-i)RrcRAo(i  n 

PEGP 1320 

01C1 

T ( I ) 3 HMATC I )/STO( 1 ) 

K TOP  13  30 

0 IC2 

i ru 

CONTINUE 

FF.GM350 

Ol  C 2 

OU  190  I 3 1 , N 

PEGP1350 

01  C5 

YMArm  3 osubu 

REGF1360 

01  C 5 

UU  100  J3l,M 

REGP 13  70 

01C6 

iso 

YHA  T ( 1 )«YFaT( I ItBHATCJ )*X( J*l ) 

REGM300 

01  C 7 

DlFHl)  3 Y(IDtP,l)-  YMAT(I) 

PEGP 1390 

U1C6 

190 

CON  riNUE 

PEGP 1500 

01C9 

Wkl re  C6, 503)  M,N 

REGF 1510 

c 

PEGP 1520 

c 

OUTPUT  THE  PEGPESSION  RESULTS 

wkGR 1530 

c 

PEGP.1550 

c 

RCGP1550 

c 

-1 

RFGP1560 

c 

IE  QJTPU1  OF  THE  X»X,  X'Y,  (X*X) 

C RESICUALS  IS  DES1PE0, 

KtGP 1570 

c 

THE  FCLLUwlNG  CCMRENT  CAROS  SHOULD 

BE  ACTIVATED  BY 

PEGP  1530 

c 

KEMJVING  T'Hfc  C’S 

PF  GP 1590 

c 

PEGP1500 

01  10 

Gti  TO  (52l,  522),ICEP 

PFGf  1510 

01  11 

521 

wr.l  ffcU,52J) 

P t GK 1 520 

01  12 

GU  TO  525 

PEGP 15  30 

0113 

522 

A RI TE(6,525) 

RCGP1550 

0115 

525 

K SO  3 SSR/SST ( I OEP ) 

PEGP 1550 

C 

WRITE  ( 6 , 5C5 ) 

PLGF15O0 

C 

OU  200  1*1, M 

REGP15/0 

C 200 

WRITE  ( o , 5 06  ) 1,  (XPXCCI  , J),  J*1,M) 

KEGP1530 

C 

WRITE  (o,515) 

KEGP  1590 

L 

DO  210  l«l,M 

REGF 1600 

C 2 10 

WR I r l ( u , 5 1 o ) 1 , XPYC ( 10EP , I ) 

PEGP 1610 

C 

WRITE  (6,507) 

REGR1620 
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% 


V 


c 

C22U 

0115 
01  16 
0117 

oi  le 

0115  230 

c 

c 

C240 

012C 

C 

c 

c • 

0121  503 

C 506 
C5l‘> 
L5uo 
L 5 07 
L^ea 

0122  505 


0123 

510 

U W4 

51  1 

C 5 12 

C 

C 51  3 

L5  In 

C 5 15 

L 5 16 

0125 

623 

0126 

524 

0127 

00  220  I * It  M 

wRI  TE  ( 6 « 500)  It (XPXIlfJ) » J«l»M) 
hi  * N — 1 

WP.I  a (6,5C9)  Nl,S57notP),.1fSSR,XM5Rf  F,K,SSEfXMSefPSO 
WRITE  (6,510)  OSUBC 
DO  230  I * 1 ,M 

*RI  Tt  (ov511)  I , BFAT  ( I ) , I , T ( I ) ,t,$TD<I) 

WRITE  (6,512) 

L>U  240  I * 1 • N 

WRITE  (6,5li)I,Y(  IGEP,l),YHAT  m,l>IFF(I) 

RETURN 

FORMATS 

FO*vaT  ( l «( l ,2  OX  , 22F"LL7(PlE  REGKF.S  $ 1 UN  ON,  13,  1 5H  VAF  I AET 
3! 5,  13H  U16EKVA7ILN5  ///) 

F'1C  4aT  (IH  , iX.AHRCw,  13, 3X,  6f  l 3 . 8/ ( 12X,6E10.0)  ) 

FUk | ( 1 (CfluA, •CCSFECTED  X**X  MATRIX1///) 

F vii  W A T (IM  ,lx,OltF  JR,  13,3X,6Eld.8/(  !2X,oE  18.3)  ) 

FJF  «»AT(  /////  1 1 X,  •CLHktCl  LO  X'*X  INVFOAF  MATRIX*  ///) 

FUR  (AT  ( in  9lx,<.MFO*,lJ,3X,oEl<I..O/(  12xtoEld.a)  ) 

FUr  (AT  (l MO,  / ,**0X,5FAf,jVA///l  W,tMS3URCt,9X,2H0F,  10X,2HSS, 

* loA  ,2HJ*a,loA,  IMF  , 1 7 X/ //  1 1 X , 5HT*)T  AL  , 9 A , 1 3 , 2X  , 6 1 fl  . B // U X, 
♦lUMRbr.RtjG  IuM,4X,  I 3 ,2  X,  if  1 J.0.//1  lX,nH=FSIO')AL,oX,I  i,  2X, 
WEld.d/////  11X,*C  IT  EFFICIENT  OF  MJlTIPIE  UETEFMINVT  I.JN  (R**2) 
*,  C l ).  d ///  ) 

FJF  (ATC22X,l2HCUtFF!ClFliTS,  15X,6HT  VALUES  , 20X  13MSTANUAF  C CFV. 

♦ VHdSLHO  * ,tl5.0) 


PEGFL630 
PEGF 1640 
REGFlo50 
«HGF 1660 
PHGF1670 
RFGPL680 
PEGF1690 
PCGP  1700 
FEGR1710 
RE  Gr  1720 
FTGR1730 
FEGF1 740 
RL-GP1750 
REGP1760 
w l I Mw  EGR  1 7 70 
RFGC 1780 
REGF1790 
REGP 1800 
RfcGF  IB  10 
R E G c 1 820 
"EGF1830 
REGF 1040 
L.  EGK  18  50 
RE  SF 1060 
* • R EOF 1 3 70 
KEGF  1 9 HO 
//IIXPEGFIBVO 
F FGC  1900 


F jf  M/.T(11X,  jhOHAT(  , I2,2M)  «,  E 16 .9 , 1 OX,  2H  T ( , 1 2 , 2b ) * , F 1 5 .0  , l 0X,3HSD(  , FCGF  19  10 

♦ 1 2 , 2H  )«, c 15. d)  ^T.GP  192  0 

FQk  *aT  ( IM1, 10X, * CbSEFYATI JW*  ,12X, • Y-CBCERVEO* ,10X, • Y-E ST  IM AT  ED*  REGP 1930 

• .n,  * DIFFERENCE  • , ////////)  FEGF1940 

F Of'  (AT  (IM  ,lUX,I3,19X,U5.8,5X,Ei5.8,5X,ri5.0)  REGF  1950 

FORMAT  ( I 2, 2X, I 3, 2X, 13)  REGR1960 

F Jr  wat  ( ///////////////,  l IX, 20HCOFRECTEO  X*Y  MATRIX////)  RcGF 1970 

FUF.vit  ( 2X,4HNHw  ,I3,JX,6E10.0/(llX,6E18.8))  R«=GF  1900 

F OF  '(AT  (*J*  , T60,  *CCS1  EQUATION*)  RCGFL990 

FOaMC  J*  , r50,*PRCCUCTICjN  FATE  EQUATION*)  3 EGR  2000 

END  RCGP2010 
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FORTRAN  IV  G L£  VFL  21  OIMFQN  DATE  ■ 73361  15/21/16 


OOC1 

SUbKCUTME  UlMfcON  l A , N , EPS , OCT , NSOLUT  ) 

OTMF 

10 

C 

OIMF 

20 

C SUBROUTINE  DiMCUN  I « A SINGLE  PRECISION  M/TRIX  INVERSION  ROUTINE 

Cl  ME 

30 

C ThaT  SfEKS  MAXIMUM  FlVI  f T L LM ENT 3 A >10  INVERTS  IN  PLACE. 

01  ME 

40 

C 

\«  ARK  AY  CUnTaIUING  M^TkIX  OF  C lEFTICIEHTS  TO  BE  INVERTED. 

OHE 

50 

c 

N*  CKJkR  UF  MATRIX  A. 

01  ME 

60 

c 

ESP3  MI  ALI  CWAUIE  VALUE  OF  MATRIX  PIVIT  BEFORE  MATRIX  IS 

Cl  ME 

70 

c 

TERMED  SINGULAR. 

DIME 

90 

c 

OET*  VmLUL  OF  DETERMINANT  CF  MATRIX  A. 

01  ME 

90 

c n 

jOLUT*  NO.  SJLL'i  ION  VECTORS  IN  THE  AUGMENTED  MATRIX. 

OIME 

100 

c 

01  Mf 

IIO 

JOC2 

DIHELSIL4  I PI  V I 2 » 2C),ICJ  I20),Y|  23I,A(5,6I 

CMC 

120 

OUC3 

IF  IN.^T.l) RETURN 

0 1 MF 

130 

0GC4 

ir  CN.  A.  1 )uO  TO  250 

O’ ME 

140 

UUC5 

DEI  3 4 1 It  1 ) 

CiME 

150 

Ou  c6 

A(  L,  l)3l./A|  l, n 

DIME 

160 

OUC7 

M 1.21  * A 11. 1)*A(192) 

01  U 

170 

UOC8 

RbTUK  1 

OIME 

180 

uocs 

250 

CD  NT  I NOE 

DIME 

190 

00  1C 

M 3 N ♦ NSOLUT 

Cl  *E 

200 

0011 

JET  3 1.2 

ni^t 

210 

oU  12 

ASSIGN  2J5  TO  IZEPC 

OIME 

220 

OU  I J 

00  iuO  K * 1 » N 

dime 

230 

JO  is 

KM1  = K-I 

DIME 

240 

00  1 5 

ircKri.Gr.oi  assign  95  to  izero 

OPE 

250 

0016 

d 1 G A 3 0.0 

CIME 

2 60 

0017 

on  loi  i»i,n 

0 1 .MF 

270 

OJ  18 

on  101  j * 1 . n 

0 1 MF 

280 

0019 

r*G  Tu  I/tRU  , 195,  205) 

01  ME 

290 

0G2C 

95 

NO  1C2  1 1 3 1 , KMl 

OIME 

300 

JU21 

1FI  L .CU.IPIVU,  II  ) .CR.J  .EO.  I P I V C 2 , 1 I))  GO  TO  101 

OIME 

310 

00  22 

102 

CONTINUE 

dime 

320 

0 0<J 

205 

CONTINUE 

CIME 

330 

0024 

IF  ( AilSCAl  i,J)).LT.R!GA)  GC  TO  101 

DI  ME 

340 

0025 

bIGA  3 AjSCAI  1.  J)  ) 

OIME 

350 

0026 

IPI VI l . K ) 3 I 

DIMF 

360 

0027 

I PI  VI2.K)  3 J 

OIME 

370 

0028 

101 

CON  fINUE 

OIME 

380 

OG29 

I F| EICA.GE.EPS)  GC  TO  201 

DIME 

390 

003C 

OET  3 C.J 

DHE 

400 

U031 

GO  TG  20 J 

oimf 

410 

0032 

201 

IR  3 IFIV11.K) 

DIME 

420 

0022 

JL  3 IFIVI2.K) 

niM£ 

430 

0034 

BIGA  3 A 1 IR, JC) 

OIMF 

440 

0025 

JET  3 LET*ttIGA 

OIMF 

450 

0026 

00  103  LL*  l,M 

Cl  MF 

460 

0037 

103 

All  l L ) 3 A<  IR,LU/BIGA 

OIMF 

470 

00  2 E 

A I I R , JC  ) 3 l.u/BIGA 

OIME 

480 

0035 

DO  IuO  LLL3  1 . n 

OIME 

490 

004C 

AJCK  « AILLL.JC) 

DIME 

500 

0041 

IF! LLL.EO.IR)  GC  TO  100 

DIME 

510 

0042 

All.LLfOC)  3 -AJCK  /B  IGA 

OIME 

520 

0043 

00  1G4  L4  3 l.M 

OIMF 

530 

0044 

IFIL4.tC.JC)  GO  TC  104 

CIME 

540 
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FORTRAN  IV  G LEVEL  21  OIMfcQN  DATE  » 73361  13/21/16 


OOA3 

A (L  LL  »L* 1 • A(LLL,l*>  -A JCK  *A ( 1 k »L* ) 

DIME 

550 

00*6 

10* 

coot  iNue 

DIMC 

560 

00  *7 

ICO 

Cun  f 1NLE 

01  Mr 

570 

00*6 

DO  IU5  I * 1 1 N 

0!  MF 

5 BO 

00** 

IK  « IPlVdtl  1 

0 I H t 

590 

0050 

1CJ  IIP)  « JC 

CIME 

600 

0031 

105 

CONTINUE 

DIVfc 

610 

OUii 

ICT  ■ 0 

01  ME 

620 

0053 

NMi  » N-i 

DIME 

630 

00  5* 

UO  106  IMfNMl 

DIME 

6*0 

0055 

1 PI  * I ♦ 1 

0 1 ME 

650 

0056 

on  10  6 J*  I P l » N 

01  MF 

660 

UO  57 

I F C U J(J) .GE. 1CJC 111  GO  TO  106 

m af 

670 

0058 

ITfcMP  * ICJ(J) 

01  MF 

630 

00  5 * 

ICJ(J)  * 1C J ( 1) 

0 1 MF 

690 

006C 

1CJ11)  - I TEMP 

01  <E 

700 

0061 

ICT  ■ ICf  ♦ 1 

DIME 

710 

00  62 

1C6 

CGr«mtE 

OME 

720 

0063 

IFC  (ICT/2)*2.NE.ICTJ  DET  - -DET 

n i mf 

730 

006* 

UU  10  7 J * 1 » M 

UI  *F 

7*0 

0065 

DO  10  H I * 1 » N 

DIME 

750 

0066 

JC  * IP1V(2«1I 

01  ME 

760 

0067 

Ik  « iPIVIltl) 

DIME 

770 

0066 

100 

Y( JC)  * AC  Ik# J1 

oi  *.r 

78u 

006* 

DO  107  K*  1 1 N 

0 I Mf 

790 

007C 

1C7 

M K , J J * Y ( K I 

dime 

300 

0071 

00  110  1*  1VN 

01  ME 

310 

0072 

UU  111  J * 1»N 

01  ME 

62  0 

00  73 

IR  - iPlVCltJJ 

DIME 

330 

007* 

JC  ■ 1P1V(2«J) 

01  ME 

8*0 

0075 

111 

YC1R)  * A ( 1 1 JC  ) 

DIME 

850 

volt 

UU  lie  K * If  N 

OIME 

860 

uO  77 

110 

ACI fk)  * YCK) 

UI  me 

870 

0078 

2 CO 

RETLkN 

01  MF 

880 

007* 

END 

DIME 

090 
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FUR1FAN  IV  G LEVEL  21 


BLK  OATA 


DATE  ■ 73361 


15/21/16 


OOC1 

BLOCK  DATA 

BLCK 

10 

C 

DICK 

20 

C 

DEE  IN  IT  ION  UF  CCNSTANTS 

QIC  K 

30 

c 

BLCK 

',0 

00  C 2 

C01VLN  /C10A/F4NC  1C) 

61CK 

50 

00C3 

CO-I^CJN  /C106/PT1  f FT2.PT3,  NAME  ( 50  ) , CHART  ( 100,  50) 

QLCK 

60 

00  C* 

LU4KLN  /C  IC7/INUK)  , SMAX , SM  IN , FMAX , FMl N 

QLCK 

70 

00C5 

*)Ai  ^ F^N/Ot  i . .40,. 42,. 43,. 44, .45,  .45, .46,  .46/ 

BLCK 

80 

oott 

DATA  FTl,FT2,P73/,.,,,*,t,  •/ 

BLCK 

90 

00C7 

DAI  A NAM£ / 50* • • / 

BICK 

100 

OOCE 

DATA  NAME  (I  6)  , NAME  ( 1 7 )f  NAME  ( 18  ),  NAME  ( IV  ),  NA  *E  ( 20  )/•-•,' M,' P • , 

BLCK 

110 

*•!• ,*-*/ 

BLCK 

120 

OOCS 

DATA  NAME  ( H I f MA  K E ( 24  ) , NA  M E (26), NAME!  2H)/,D',,£,,,L,,,F'/ 

BLCK 

130 

001C 

DATA  ChA<l/5J00*»  •/ 

BLCK 

no 

0011 

DAT  A lM)/lUO.Of  ICC,  C/ 

BLCK 

150 

00  u 

DATA  S MAX ,SMIN,FMAX ,T Ml N/O.C, 0.0,0. 0,0.0/ 

hLCK 

160 

JO  12 

END 

BLCK 

170 

120 


Appendix  G.  MACHOP  Program  Documentation 


G.l  Program  Description 
G.1.1  Synopsis 
G.1.2  MACHOP  Block  Diagram 
G.l. 3 Program  Sequence 

G.1.4  Program  Requirements  and  Restrictions 

G.2  Input-Output  Descriptions 
G.2.1  Program  Set-ups 
G.2. 2 Input  Formats 
G.2. 3 Input  Requirements 
G.2. 4 Output  Description 

G.3  Systems  Material 
G. 3 . 1 Flow  Chart 

G.3.2  Glossary  of  Important  Variables 


G.l  Program  Description 
G.1.1  Synopsis 

The  MACHOP  program  is  written  in  FORTRAN  IV  for  an  IBM  System 
360/65  Data  Processing  System  using  the  standard  IBM  FORTRAN(G) 
compiler.  The  following  documentation  outlines  all  aspects  of  this 
application  unique  to  this  program. 
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G . 1 . 2 MAC  HOP  BLOCK  DIAGRAM 
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G . 1 ♦ 3 Program  Sequence 


MACHOP  consists  of  one  main  routine  and  eight  subroutines.  The 
normal  sequence  is: 

MAIN 

NUCYCL 

EVOP 

ORGPTS 

PLOT  PM 

REGPLT 

REGRES 

DIMEQN 

BLOCK  DATA 

Functions  which  are  required,  but  are  not  supplied  with  this  program, 
are  available  in  most  IBM  FORTRAN  IV  compilers.  These  are  as  follows: 

Square  Root  (SQRT) 

Absolute  Value  (ABS) 

Maximum  (AMAX1) 

Minimum  (AMIN1) 

Natural  Logarithm  (ALOG) 

G.1.4  Program  Requirements  and  Restrictions 

The  MACHOP  routine  uses  a building  block  approach  to  examine  the 
response  surface.  The  output  from  a cycle  must  be  included  in  the  input 
for  the  next  cycle  in  order  to  evaluate  correctly  the  new  observations.  If 
at  any  point  in  this  process,  output  from  the  previous  cycle  is  omitted 
from  the  input  of  a subsequent  cycle,  the  results  will  be  invalid. 

The  following  restrictions  apply  to  the  MACHOP  program: 

1.  The  number  of  spedds  defined  for  this  operation  must  be  less 
than  or  equal  to  100. 

2.  The  number  of  feeds  defined  for  this  operation  must  be  less 
than  or  equal  to  50. 

3.  Observations  must  be  inputted  in  groups  of  four. 

4.  No  constraints  other  than  the  limits  on  the  speeds  and  feeds  are 
taken  into  account  by  this  program. 
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G.2  Input-Output  Descriptions 


G . 2 , 1 Program  Set-ups 


(7* 

[PUNCHED  OUTPUT  FROM  PREVIOUS  CYCLE 
NEW  OBSERVATIONS  (4) 


OLD  OBSERVATIONS 


f LIMITS  FOR  SPEEDS  & FEEDS 
[DEFINED  FEEDS 
fffOF  FEEDS 
("DEFINED  SPEEDS 
f#  OF  SPEEDS 
(^TOOLING  COSTS 
('MAIN  problem  card 


i 


V 


f TITLE  CARD 
('//GO.SYSIN  DP  * 


ILOCK  DATA 


fblMEQN 


(REGRES 


fREGPLT 

["PLOTPM 


Control 

Cards 


Program 

Deck 


User 

Supplied 

Data 


(ORGPTS 

(EVOP 

('nucycl 

f MAIN 

(V/FORT. SYSIN  DD~* 
| '//  EXEC  FORTRAN 

07  JOB 


I Control 
j Cards 


MAC  HOP  Deck  Set-up 
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G.2.2  Input  Formats 


The  input  formats  for  the  MAC  HOP  program  are  given  in  Figure  G.l. 
These  include: 

1.  Title  Card  (20A4)  (Required) 

The  Title  Card  consists  of  user  supplied  identification 
information  punched  in  columns  1-80. 

2.  Main  Problem  Card  (513,  T21,  5F10.4)  (Required) 


cc 

1-3 

Variable 

IPH 

Description 

Current  Phase  Number  (Right  justified) 

4-6 

ICYC 

Current  Cycle  Number  (Right  justified) 

7-9 

IOP 

Operation  Type:  IOP=  1:  single  operation 
(Right  justified)  2;  multiple  opera- 

tion 

3:  numerically 
controlled 
operation  (N/C) 

10-12 

I MAX 

1:  If  limits  for  speed  and  feed  are 
specified  (Right  justified) 

0:  Otherwise 

13-15 

NTLC 

Number  of  different  tools  (Right  justified) 

21-30 

RLO 

Rate  for  labor  and  overhead 

31-40 

AVES(COST)  Preliminary  estimate  of  the  standard 
deviation  of  the  cost  observations 

41-50 

AVES(PA) 

Preliminary  estimate  of  the  standard 
deviation  of  the  production  rate  observa- 
tions 

51-60 

BSP 

Base  speed  for  N/C  operations  (Must 
be  present  when  IOP=3,  cc  7-9) 

61-70 

BFD 

Base  feed  for  N/C  operations  (Must  be 

present  when  IOP=3,  cc  7-9) 
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3 .  Tool  Cost  Card  (Required) 


cc 

Variable 

Description 

1-10 

TLC(l) 

Cost  of  a tool 

edge  for  tool  no.  1 

11-20 

TLC  (2) 

Cost  of  a tool 

edge  for  tool  no.  2 

21-30 

TLC(3) 

• • « 

Note:  there  must  be  exactly  NTLC  (cc.  13-15,  M.P.C.)  tool 
edge  costs  reported. 

4.  Speed  Environment  (I3/8F10.4)  (Required) 

Card  1:  cc  1-3,  NSPD,  Number  of  speeds  defined  in  the  speed  v 

environment  (NSPD^  100;  Right  justified) 

♦Cards  2,  3, . . Defined  speeds  (rpm)  for  the  environment 

in  ascending  order.  (There  must  be  exactly  NSPD 
speeds  inputted.) 

5.  Feed  Environment  (I3/8F10.4)  (Required) 

Card  1:  cc  1-3,  NFD,  Number  of  feeds  defined  in  the  feed 

environment  (NFD^  50;  Right  justified) 

♦Cards  2,  3, . . .:  Defined  feeds  (ipr)  for  the  environment  in 

ascending  order.  (There  must  be  exactly  NFD 
feeds  inputted.) 

6.  Feed-Speed  Limits  (4F 10. 4)  (Optional) 

This  card  must  be  present  if  IMAX  = 1,  cc.  10-12  M.P.C. 


1-10 

11-20 

21-30 

31-40 


Variable  Description 

SMAX  Maximum  speed  allowed  for  this  operation 

SMIN  Minimum  speed  allowed  for  this  operation 

FMAX  Maximum  feed  allowed  for  this  operation 

FMIN  Minimum  feed  allowed  for  this  operation 


♦For  N/C  operations  the  speeds  and  feeds  are  the  percentage  overrides 
to  be  considered  in  this  analysis  (not  fractional  equivalents) . 


126 


7.  Observation  Cards  (8F10. 4)  (Required) 


cc 

1-10 

Variable 

X(l,i) 

Description 

Observed  speed  setting  (rpm) 

11-20 

X(2,i) 

Observed  feed  setting  (ipr) 

21-30 

X(3,i) 

Number  of  parts  produced 

> 

31-40 

X(4,i) 

Time  (in  minutes)  expended  on  this 
operation 

41-50 

X(5 , i) 

Number  of  tool  edges  for  tool  #1 

* 

51-60 

X(6,i) 

Number  of  tool  edges  for  tool  #2 

61-70 

X(7,i) 

• • « 

NOTE: 

Tool  edges  must 

be  reported  for  exactly  NTLC  (cc. 

M.P.C.)  tools  for  each  observation. 


G.2.3  Input  Requirements 

1.  General  input  requirements  for  the  use  of  the  MACHOP  routine  are: 

a.  All  variables  indicated  to  be  right  justified  are  integer  values 
and  must  not  have  decimal  points  punched. 

b.  All  other  values  are  real  numbers  and  should  have  their  decimal 
points  punched. 

c.  Observations  must  be  read  in  groups  of  four. 

d.  Phase  and  cycle  numbers  must  be  accurate.  (Phase  indicates 
which  set  of  four  points  is  being  observed  and  cycle  indicates 
the  number  of  observations  which  have  been  taken  at  each  point.) 

e.  Punched  output  from  one  cycle  must  be  included  in  the  input 
of  the  next  cycle. 

2.  Specific  input  requirements  for  the  MACHOP  routine  differ  for  each 

category  of  input.  For  clarity,  sample  deck  set-ups  are  defined 

for  each  major  situation. 
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a.  Case  I:  Initial  phase  and  cycle  with  an  estimate  of  the 
standard  deviation. 


OBS.  #4 


f 


f 


OBS.  #3 


OBS.  #2 


OBS.  #1 


S/F  LIMITS 


FEEDS  DEFINED 


# FEEDS 


SPEEDS  DEFINED 


# SPEEDS 


f 


TOOL  COSTS 


MAIN  PROBLEM  CARD 


TITLE  CARD 


1 .  Title  Card  (Required) 


2.  Main  Problem  Card  (Required) 


(cc  1-3) 

(cc  4-6) 

(cc  7-9) 

(optional) 

(cc  13-15) 

(cc  21-30) 

AVES  (COST)  = Standard  Deviation  Estimate  (cc  31-40) 
AVES  (P/R)  = Standard  Deviation  Estimate  (cc  41-50) 


IPH  = 1 
ICYC  = 1 
IOP  = 1,2,3 
IMAX 
NTLC 

RLO  = Value 


3.  Tool  Costs  (Required) 

4 . Speed  Definitions  (Required) 

5 . Feed  Definitions  (Required) 


6.  Speed  and  Feed  Limits  (Optional) 

7.  Observations  (Exactly  four  (4)  required) 
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b. 


Case  II:  Initial  phase  and  cycle  without  an  estimate  of 
the  standard  deviation.* 


OBS.  #4 


Jf 


OBS.  #3 


OBS.  #2 


OBS.  #1 


S/F  LIMITS 


FEEDS  DEFINED 


# FEEDS 


SPEEDS  DEFINED 


# SPEEDS 


TOOL  COSTS 


MAIN  PROBLEM  CARD 


TITLE  CARD 


1.  Title  Card  (Required) 

2.  Main  Problem  Card  (Required) 

IPH  = 1 (cc  1-3) 

ICYC  = 1 (cc  4-6) 

IOP=  1,2,3  (cc  7-9) 

I MAX  (optional) 

NTLC  (cc  13-15) 

RLO  = value  (cc  21-30) 

3.  Tool  Costs  (Required) 

4.  Speed  Definitions  (Required) 

5 . Feed  Definitions  (Required) 

6.  Speed  and  Feed  Limits  (Optional) 

7.  Observations  (Exactly  four  (4)  required) 


*A  limited  output  will  result  since  no  estimates  of  the  standard  deviation 
are  available. 
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c.  Case  III:  Non-initial  cycle. 


f 


P/O  FROM  PREVIOUS  CYCLE 


NEW  OBSERVATIONS  (4) 


OLD  OBSERVATIONS 


S/F  LIMITS 


FEEDS  DEFINED 


# FEEDS 


(l 


SPEEDS  DEFINED 


# SPEEDS 


TOOL  COSTS 


MAIN  PROBLEM  CARD 


TITLE  CARD 


1.  Title  Card  (Required) 


2 .  Main  Problem  Card  (Required) 

IPH  = Phase  Number  (cc  1-3) 

ICYC  = Cycle  Number  (cc  4-6) 

(NOTE:  IPH  = 1 and  ICYC  = 1 not  used) 
IOP  = 1,2,3  (cc  7-9) 

I MAX  (optional) 

NTLC  (cc  13-15) 

RLO  = value  (cc  21-30) 


3.  Tool  Costs  (Required) 

4.  Speed  Definitions  (Required) 

5.  Feed  Definitions  (Required) 

6.  Speed  and  Feed  Limits  (Optional) 

7.  Old  Observations  (Required) 

8.  New  Observations  (Exactly  four  (4)  required) 

9.  Punched  Output  from  Previous  Cycle  (Required)* 

♦Note:  After  each  cycle,  the  old  punched  output  is  discarded  and  is 
replaced  by  the  new  observations  and  new  punched  output. 
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Figure  G.  1 MACHOP  Input  Formats 
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G.2.4  Output  Description 


Output  from  the  MACHOP  program  generally  consists  of  less  than  2000 
printed  lines  and  five  punched  cards.  The  output  is  as  follows: 

a.  Input  information  such  as  feed-speed  limits,  labor  and  overhead 
rates,  tool  costs,  and  prior  estimates  of  standard  deviations 

of  the  response  variables  (if  any) . 

b.  The  specified  feed- speed  environment  represented  "graphically.  " 

c.  The  cumulative  input  data  and  the  computed  responses  (cost 
per  piece  and  production  rate) . 

d.  A table  of  the  computed  responses  for  this  phase  and  cycle. 

e.  Evolutionary  operation  calculations  for  both  the  cost  per  piece 
and  production  rate  responses. 

f.  A recommended  set  of  operating  conditions  for  the  next  experiment 
(which  may  be  identical  to  the  current  set),  based  on  the  cost  per 
piece  analysis. 

g.  A graphical  representation  of  the  operating  conditions  recommended 
in  f. 

h.  The  recommended  set  of  operating  conditions  for  the  next 
experiment  (which  may  be  identical  to  the  current  set),  based 
on  the  production  rate  analysis. 

i.  A graphical  representation  of  the  operating  conditions  recommended 
in  h. 

J.  The  results  of  the  regression  analysis  and  the  accompanying 
analysis  of  variance  table  for  the  cost  per  piece  analysis. 

k.  The  results  of  the  regression  analysis  and  the  accompanying 
analysis  of  variance  table  for  the  production  rate  analysis. 

t . The  predicted  response  surfaces  for  cost/piece  and  production 
rate  calculated  from  the  respective  prediction  equations. 

Optional  output  includes  a listing  of  observations,  the  predicted  values, 
the  residuals,  the  matrices  X'X  and  (X'X)-*,  and  the  vector  X'^;. 
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G«3  Systems  Material 


G,3«l  Flowchart 


MAIN  PROGRAM 


SUBROUTINE  NUCYCL 


(nu^T) 


•“(  ORGPTS  j 


. - -(^PLOTPM^) 


SUBROUTINE  EVOP 


i 
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<J1 
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SUBROUTINE  PLOTPM 


SUBROUTINE  ORGPTS 


U> 

"NJ 


0 


SUBROUTINE  REGPLT 


SUBROUTINE  REGRES 


G.3.2  Glossary  of  Important  Variables 


Variable  Name 

Description 

ANS 

Vector  of  regression  results  returned  from  REGRES 
subroutine 

AVE 

Vector  of  the  averages  in  the  EVOP  routine 

AVES 

Estimate  of  the  standard  error  in  the  EVOP  table 

B 

Vector  of  estimates  of  the  b-coefficients  in  the 
regression  equation 

CHART 

100  x 50  matrix  containing  the  plot  output  in  the 
ORGPTS  routine 

DIFF 

Vector  of  the  differences  between  the  average  cost 
at  a given  point  and  the  newly  observed  cost 

FEED 

Vector  containing  all  defined  feed  settings  in 
ascending  sequence 

F4N 

Vector  of  the  values  of  f^  n as  defined  in  Box 
and  Draper  [ 3 ] 

ICYC 

Cycle  number 

IFD 

Index  containing  the  current  feed  setting 

IPH 

Phase  number 

ISP 

Index  containing  the  current  speed  setting 

NFD 

Number  of  feed  settings  defined  in  the  environment 
of  this  operation 

NOBS 

Number  of  observations  in  the  current  cycle 

NS  PD 

Number  of  speed  settings  defined  in  the  environment 
of  this  operation 

RLO 

Labor  and  overhead  rate  as  reported  on  the  main 
problem  card 
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Variable  Name 


Description 


SPEED 

Vector  of  all  defined  speed  settings  in  ascending 
sequence 

SUM 

Vector  of  the  sums  of  the  observations 

SUMS 

Sum  of  the  estimates  of  the  standard  error  of  the 
observations 

TLC 

Tool  edge  cost 

TLCH 

Observed  number  of  tool  edges  used 

X 

Matrix  (5xn)  of  the  input  observations,  both  old  and 
new.  For  observation  n, 

X(l,n)  = speed  (rpm) , 

X(2 , n)  = feed  (ipr) , 

X(3,n)  = number  of  parts  produced, 

X(4,n)  = production  time  in  minutes,  and 
X(5,n)  = number  of  tool  edges  used. 

Y 

Matrix  (2xn)  of  the  performance  indices,  both  old 
and  new.  For  observation  n, 

Y(l,n)  = cost  (dollars/piece),  and 
Y(2,n)  = production  rate  (piece s/minute) . 
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